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This paper describes new equipment and its applications in the precise determination 
of the spectral radiant energy emission of light sources having low intrinsic brightness, in 


particular fluorescent lamps. 


light sources having superimposed bright lines is discussed in detail. 
number of commercially available fluorescent sources. 
for example absolute spectral emissivity measurements on tungsten, spectral 


problems, 


The evaluation of the spectral radiant energy distribution of 


Data are given on a 
Applications are noted in other 


reflectivities, spectral transmittance measurements on dark shade welding glass, total ozone 


studies, 


I. Introduction 


The production of light (or other radiant energy) 
through the frequency changing process within a 
fluorescent lamp tube, wherein the generated radiant 

ergy (principally) of wavelength 2537 A is emitted 
through a band of frequencies (wavelengths) extend- 

over an appreciable spectral range, results in 
sources of high luminous (or radiant) efficiency but 
of low intrinsic surface brightness. The low radiant 
ensity per unit area makes it difficult to determine 
irately the spectral radiant intensity with ordi- 
nary spectrophotometric apparatus. Precise meas- 
ements are further complicated by the presence of 
mercury emission lines superimposed upon the con- 
tinuous fluorescent spectrum. These lines contribute 
ppreciably to the color of the lamps and hence to 
their relative spectral quality and must therefore be 
irately evaluated in any precise spectral energy 
surements 

lt is a fortunate circumstance that many of the 
ivailable phosphors have a maximum spectral ex- 

tion sensitivity near the wavelength 2537 A in 
mercury emission spectrum, especially since ap- 
oximately one-half of the wattage dissipated 
thin the low pressure mercury are discharge is 
emitted at this wavelength. Otherwise the surface 
brightness of fluorescent lamps would be even lower. 

The introduction of fluorescent lighting into the 
many industries and pursuits of man has aroused 
much mterest in its relationship to colors, whether 

be that of a paint, a dyed fabric, a chemical reac- 

or even an individual’s facial appearance [{1, 2]. 

is recognized that the simple matching of a 

rescent lamp color to a particular tungsten fila- 

color temperature, or to average daylight, does 
nsure equal color appearance of the illuminated 

ct. A food, for example an egg, an orange, or 
ipple may appear unappetizing under a fluorescent 
amp, Which matches perfectly average daylight in 
all color. A vegetable, flower, or garment may 
wise appear of different hue. The spectral qual- 

of the two sources also must be considered since 


ires in brackets indicate the literature references at the end of this paper 


or measurements in other situations where low intensities may be involved. 


the colors of the materials, whether they be food, 
flowers, clothing, household furnishings, or a human 
face depend upon both the optical spectral reflective 
characteristics of the material and upon the relative 
spectral radiant energy distribution of the source of 
illumination. 

The recent widespread use and diversity of colors 
of fluorescent lamps, as improvements and changes 
are made in the composition or treatment of the 
phosphors, have resulted in an increased demand for 
more precise information regarding their spectral 
emission characteristics. It is the primary purpose 
of this report to describe a new instrument and its 
application to this problem. 


II. Instruments and Methods 
1. Spectroradiometer 


A double quartz prism spectrometer, manufactured 
by the Farrand Optical Co., has been modified 
slightly for use in this investigation. (See fig. 1 
This instrument employs two sets of paired 30- 
deg crystalline quartz prisms. The are of 
30 cm focal length, have 75-mm free aperture, and 
are doublets, constructed of quartz-lithium fluoride 
and achromatized at two positions within the 
ultraviolet spectrum (about 220 and 400 may). 
The slits are each manually adjustable. The 
entrance and exit slits are straight, the center slit 
being curved. Changes in wavelength are 
complished through a tandem drive that rotates 
the four prisms arranged in Thollon type mountings 
All other optical components of the instrument are 
fixed in position. 

For recording purposes a three-speed reversible- 
gear drive from a small synchronous motor was 
attached to the center wavelength adjusting knob 
(shaft These three speeds allow scanning the 
spectrum from 200 to 780 my in about 4, 8, or 12 
min in either direction. the calibrating 
factors are slightly different for the two directions 
all measurements, in the present investigation, 
were made with increasing wavelength. Cams and 
switches are provided to automatically stop the 
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2uc- 


Since 


437 





and doult 
oyed in fl 


e monochro 


emp iorescent 


motor drive at either end of the spectrum and to 
give index marks on the recorder chart, by means of 
an auxiliary pen, at 18 positions (equally spaced in 
wave number) through the spectrum. For manual 
operation the motor drive is simply disengaged 

A sector disk having 17 equally spaced openings 
and operated at 1,800 rpm by a synchronous motor 
modulates the light beam at 510 c/s before it enters 
the spectrometer. This frequency was chosen as 
one of the possible ones between 500 and 1,000 
cycles that could be easily obtained by means of the 
disk and motor setup and that was incommensurate 
with any of the harmonics of the light and power 
frequencies The use of a modulated light beam 
makes possible the employment of this instrument 
in a fully lighted room having normal light fluctua- 
tions, without noticeable effect upon the recorded 
data, as only changes within the light beam path 
are registered and then only if at a frequency near 
that of the light modulation. Disk openings of 
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circular shape give approximately sinusoidal | 
modulation. They require a disk of larger diamy 
than do angular openings which result in the gen 
tion of an approximately square wave. Eit 
type gives satisfactory performance 
The photoelectric employed in 

instrument in the present work are a standard | 
sensitive (RCA 935) and a red sensitive (GE PJ-1 
to cover the range of about 250 to 780?my. A qui 
tube having an S-4 (blue surface 
available for use with the short wavelengths. Ot 


detectors 


sensitive 


types of standard or special tubes are of cours: 
adaptable for use in this instrument 


2. Tuned Amplifier 


A tuned amplifier (see fig. 2) peaked at 510 
has certain characteristics that may be noted 
The input stage employs a specially selected tube 
high signal to noise ratio, the type 1620. This st; ore 
circuit is designed for high signal gain with litth 
attention being given to frequency shaping. 

The second stage, employing a 65J7 tube, contains 
a tuned secondary consisting of a high @ variabl 
inductor shunted by a fixed mica condenser. This 
combination results in a band pass of about 100 
cycles. The amplifier is stable at all gain settings 
yet has a sufficiently narrow pass band to reduc 
noise and power line frequencies to a relatively iow 
level. The pass band peak is sufficiently broad 
however, when properly adjusted to 510 cycles, that 
normal fluctuations in the power line frequency ar 
unnoticeable. 
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(mplifier gain is controlled by a tandem duogrid 
riable resistance linkage in the second and _ third 
The gain is fixed in the first and fourth 


ages 


ves, except that provision is made for substituting 
proper value of grid resistor in the first stage for 


proble m at hand. Ample shielding of the first 
ve and of the second and third stage grid connec- 
ns reduces pick-up noises to a minimum. Suf- 
ent space between amplifier stages (which are 
ranged in line) and individual plate-supply filter- 
tworks reduces positive amplifier coupling to a 
egligible value. The use of an isolated power- 
pply chassis eliminates transformer magnetic fields 
thin the amplifier, thus reducing 60- and 120-cycle 
hum 
The output of the amplifier is connected with a 
vacuum-tube voltmeter (Ballantine model 300), or 
after rectification by a 6H6 twin diode the resulting 
direct current is fed through a filter and into a 
standard commercially available Leeds & Northrup 
Speedomax recorder. The alternating-current out- 
put is linear for all values below about 20 v.  Re- 
corded direct-current values are slightly nonlinear as 
the result of contact potential within the 6H6 diode 
rectifier. Increased linearity is obtained with in- 
creased load resistance. Various circuits have been 
employed for eliminating nonlinearity, among these 
being the use of a balanced direct-current amplifier 
1] or balanced cathode follower. Usually, however, 
the resulting zero drift counterbalances any gains 
obtained over the use of the standard diode rectifier 
circuit, which is illustrated in figure 2. With this 
circuit the zero is stable, and a linearity correction is 
easily made. However, it is often preferable in 
precision work to also read the intensity values 
directly on the alternating-current vacuum-tube 
voltmeter and manually plot the resulting data. 
Most of the data displayed in this report were 
obtained by the latter method, the recorder curves 
being employed primarily as supplementary aids. 


3. Standards of Radiant Energy 


The precise evaluation of radiant energy requires 
the use of a standard on which the measurements are 
based. At this Bureau all work of this type is 
necessarily referred to a fundamental black body. 
In practice, however, thermal standards [5] that have 
been prepared in the form of carbon filament lamps 
are employed. In the present work secondary 
standards in the form of a special tungsten-in-quartz 
lamp [6] and a quartz mercury-are lamp [7] are 
required to obtain a precise calibration of the spec- 
rometer for use with a continuous source having 
superimposed bright lines. 


(a) Tungsten-Filament in Quartz Lamp 


lhe tungsten-filament-in-quartz lamp employed 
this work (see fig. 3) is of special construction by 

the Westinghouse Electric Corp. It is based upon a 
oratory model described in a previous report [6]. 
constructed of drawn tungsten wire arranged in 
ane formed by a series of “hair-pin’” loops. The 


Figure 3. Tungsten-filame nt-in-quartz standard lamp. 


filament wires are evenly spaced, uncoiled, and 
covered by shields at the looped ends in order to 
insure approximately uniform temperatures through- 
out the filament area. The use of a high quality 
fused silica bulb permits a high transmittance of 
radiant energy throughout the region from below 250 
to above 2,000 my. The high quality of the fused 
silica envelope allows the employment of a projection 
system where the use of a single filament wire is 
desired. However, in the present work no projection 
system is needed or desired. The light from the 
complete lamp is simply allowed to fall directly upon 
the entrance slit of the spectrometer. 

The spectral distribution of the radiant energy 
from a lamp of this type is a function of the filament 
temperature, the emissivity of the material compos- 
ing the filament, and the transmittance of the lamp 
envelope. For this purpose, therefore, black-body 
radiant energy, the emissivity of tungsten, and the 
transmittance of the fused-quartz envelope must be 
considered. The relative spectral radiant intensity 
from a black body depends only upon the tempera- 
ture. Data for any desired temperature may be 
caleuiated through the use of the Planck radiation 
law. However, data may be found in the literature 
(24 to 26]. 

The emissivity of tungsten has been determined 
by a number of workers [8 to 11] for different tem- 
peratures. Until better data are available we are 
necessarily making use of a weighted mean curve of 
the published values for a filament temperature of 
about 2,900° K [6]. These values are probably 
adequate in radiometric work for a temperature 
range of 200 to 300 deg. The values employed in 
the present work are tabulated in table 3 for selected 
wavelengths. 

As the transmittance of the fused silica envelope 
in the spectal region of 250 to 2,000 my is nearly 
uniform, no correction is necessary in most practical 
applications of the lamp. Any deviation from uni- 
formity is no doubt exceeded by lack of accurate 
knowledge on relative spectral emissivity within the 
spectral range employed in the present work. There- 
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fore, in the present case no note is taken of the (b) Mercury Arc Line Source 


transmission characteristics of the lamp-envelope. 
The method for the absolute calibration of this lamp 


is displayed in tables 1, 2, and 3. 


The quartz mercury are lamp employed as a st 
ard line source has been previously described 
where [7]. It is a commercially available mer 
pool 110-v direct-current Uviare having about 
TaBie 1. Calibration of quartz mercury arc lamp for radi- | third of each end of the lamp tube shielded so 

OM CREED EF OS mp & SEES effects of end-blackening and arc-unsteadiness 
the two ends of the tube are eliminated from 

Output of thermal standard of radiation, («w/em *) 7.9 observations. 

Galvanometer deflection, cm 9. 5 

uw/em? per cm deflection 3. 956 The intensity of the principal spectral lines { 


Galvanometer deflection for He arc lamp (with G5860 filter, 


t 50.7 en 365.5 to 578 mu are evaluated by means of a photo- 


transmittance 0.324 at 365.5 mw) at 50.7 cn 


Radiant energ (mainly 65. 5 mg)* from Hg are through filter, 3 electric radiation meter consisting of a 935 photot hy 


pw i/cm at 50.7 em 


Equivalent radiant energy of 365.5 ma from lamp, sw/em? (at 2 (and associated electronic amplifier), a surface ther- 

Equivalent radiant energy of 365.5 ma from lamp, »wiem*at | 2 mopile and associated galvanometer, and a grou; 

aL. _____.. oe 2.5 glass filters. First, the absolute intensity of wa 
length 365.5 muy is determined by the thermopile 

a ry a ey ee ee ... | galvanometer through the use of the Corning filt: 
eS ee ee or lee . 5860, and a standard of thermal radiation [5]. 

, table 1). 

The relative intensities at the lines 405, 435.8 
546.1, and 578 mu are next obtained simply through 
sets of measurements employing the photoelect 
meter, or the thermopile and galvanometer, and a s 
of filters suitable for isolating the various mercury 
emission lines in groups of one, two, or even three an 
through a calculation of the resulting simultaneous 
equations. The intensities of the mercury Il nes 
tabulated in column 6 of table 4 were thus obta 
at the distance of 72.6 cm from the lamp tube. 


TaRLe 2. Calibration of 935 phototube radiation meter (pho- 
totube covered with G5860 filter 


TABLE 3 aNen @ Hancere Cungen filament in qua (c) Evaluation of the Spectral Line Intensities in Terms of th 


ctrometer * Continuous Spectrum Parameters 
Since the spectral energy distribution from 
fluorescent lamp consists of a continuous spect! 
upon which is superimposed a series of bright lines 
a special method must be developed for precise eva 
uation of the respective components concurrenth 
present. Such an evaluation has been calculated by 
Mary P. Lord [12] and employed in revised form 
Charles W. Jerome [13]. In the present work, hov 
ever, the evaluating factor has been experimenta 
determined through the use of two standards of 
diation—a tungsten-filament-in-quartz lamp an 
mercury arc lamp. It was felt that the mechan 
operation of a double monochromator could not 
faithfully relied upon to follow a calculated relatior 
ship between line and continuous spectrum intens!- 
ties for use in precise work. Indeed, experience iF 
shown that neither of the instruments in this lab 
tory follows precisely a calculated intensity relat 
ship. 
In order to obtain a relationship between 
response for spectral lines relative to a contim 
_ | spectrum, definite radiant energy intensities wv 
: Sion of 49,600 was 7.49 pa, with the lamp 17,1 em from the meter." admitted to the entrance slit of the spectrom 
7 ann part of this reading resulting trom the Sap band centered at 85m s | for the two types of radiant energy. For this | 
"Therefore, the part of the microammeter reading due to the 365-ms band is | pose the mercury are and the tungsten filan 
Os hexciore, from table 2. the absolute value of this 365-my band as measured is | lamps were alternately placed at the calibr: 
2.13 swjcm, a8 shown in column 8. | (Calculated values for the mens- | distances from the spectrometer slit, which 
» absolute energy value for the ngsten filament-in-quartz lamp for 365-my covered with a plate of finely rround fused S 


band accordingly is 2.13/0.32 w/em?. Other values for the other . . . . 
Bn Nan oy Be ey ety Cy ls Phe nm glass to insure uniform and complete illumina 





er 4. Calibration of spectrometer for integrated response of 
ctral lines, using standard tungsten-in-quartz and quartz- 
cury arc lamps 


Ratio, 
line to 
continu 


Spectrom- 
ur eter 
lamp * recorder ous 
inten- integrated spectrum 
sities response 


Spectrom He 
eter 
recorder 
integrated 
response 
for 5-my 
band 


re Spon se 

tor columns 
0.25-mm 7 
4x7 


slits 
x 


tions 


om of the instrument aperture. The respective responses 
ctrun were obtained as area integrations of the plotted 
lines or recorded) data for 5-my bands for the tungsten 
bute lan p and for the complete triangle in the case of each of 
rently the spectral lines for the mercury arclamp. The values 
ed by ius obtained for the four intense mercury emission 
rm 1) nes (ler 0.25-mmi slits) within the visible spectrum are 
. how- valuated in table 4. From this it appears that (for 
ntaiy the spectromeier slit widths employed in most of the 
of ra- this work; namely 0.25 mm) the spectral lines give 
ANG responses on the average about 0.896 times that for 
anica equal amount of continuous spectral radiant 
Lot be rgv when integrated in this manner and based on 
ation- )-myu spectral interval. For wider spectrometer 
itensl- slits this factor is slightly larger, and for narrower 
ce as slits slightly smaller. In the evaluation of spectral 


nes relative to continuous spectra the integrated 
ies for the lines are therefore divided by 0.896 
multiplied by 1.12) in the present investigation 
en 0.25-mm slits are employed. 
4. Method Employed in Determining Spectroradi- 
ometer Response Factors 


lhe practical use of a spectrometer for the de- 
mination of the relative spectral energy distri- 


thin the slit width range of 0.15 to 0.50 mm a change of 
chanical slit widths resulted in a change of about 10 percent 
between the line and cont 


100 percent in 
in the intensity 


NUOUS Sper tra 


bution of sources of radiant energy is dependent 
upon an accurate evaluation of all factors that 
affect the instrument output in relationship to the 
energy incident upon the entrance slit. Since in 
the apparatus herein described no condensing lens 
or mirror is employed between the lamp and the 
spectrometer slit, the radiant energy from the 
entire lamp, or from a representative portion of it 
in the case of a long fluorescent tube, enters the 
spectrometer without having undergone spectral 
distortion. However, within the spectrometer the 
lenses and prisms, and even the surface reflections, 
produce selective wavelength attenuation of the 
light beam. Furthermore, selective attenuation 
results from uncorrected errors in prism angles in 
their relationship to unavoidable mechanical defects 
in prism rotation and to slit alinement (especially 
in a three-slit instrument). The refractive char- 
acteristics of the prism material itself introduces a 
variable slit width factor. Differential loss of light 
at the margins of the optical apertures of the instru- 
ment as it is driven through the spectrum produces 
a selective spectral response effect. Hence, it is 
important that the light beams of the standard and 
unknown sources not only have closely the same 
geometrical characteristics, but that at all times the 
illuminated area remain well within the optical 
limits of the instrumental apertures. The detector, 
especially when it is a conventional phototube, has 
a response characteristic that is selective with 
wavelength. 

The various factors that determine the response 
characteristic of a spectroradiometer may be deter- 
mined by approaches from a number of directions 
and combined in the reduction of the data. Some- 
times, however, it is difficult to obtain the required 
data for individual factor determinations. Partly 
because of this difficulty, but mainly because it was 
considered a more accurate method, all factors having 
to do with the spectral output response of the instru- 
ment were combined and determined experimentally 
by means of the standard tungsten-filament-in- 
quartz lamp (supplemented by a color temperature 
standard [27 to 31]) and designated as the spectral 
response factor of the spectroradiometer. The spec- 
tral response factor is of course different for different 
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Ficure 4. Relative spectral factors for 
photoelectric spectroradiometer for selected phototubes 
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mechanical slit widths or for different phototubes 
see fig. 4 A slightly different curve is obtained if 
a color temperature standard is employed since the 
approximate spectral energy distribution associated 
with a standard of color temperature is somewhat high 
in both red and blue relative to the green [8]. In 
the present work the temperature standard (tungsten- 
filament-in-quartz lamp) is employed for all wave- 
lengths shorter than 510 my. Above this 
wavelength a mean between the temperature and 
color temperature standard is emploved since much 
of the interest in the present work is in the field of 


about 


color. 


5. Method Employed in Evaluating the Spectral Data 


The chart of the emission spectrum of a fluorescent 
lamp as observed re presents not only the character- 
istics of the lamp emission but also such modifica- 
tions as result because of the mechanical slit width 
and the relative spectral response factors for the 
instrument and phototube See fig. 4 

The mechanical slit width affects not only the 
relative spectral response factors, but also the ratio 
of the integrated areas for the bright lines relative 
to the continuous spectra. Furthermore, the wider 
the slits the wider the triangle bases associated with 
the mercury lines and the pure the 
continuous spectrum becomes In the present work 
the mechanical slit widths were set at values such 
that there was a distinct separation between the 546- 
and 578-mu lines. If the slits are set too narrow the 
transmitted energy reaches a level so low that signal- 


less 


emission 


to-noise ratio suffers. In the present work the 
mechanical slit widths chosen gave spectral-slit 
widths of about 5 to 10 my within the visible 


spectrum, 

The first step in the evaluation of the observed 
spectral data is to plot the smooth curve representing 
the continuous spectrum. (See fig. 5.) In this 
illustration the observed data in volts as read on 
the vacuum tube voltmeter are plotted (circles) as 
a function of wavelength. The solid curves represent 
the envelope of the energy distribution as drawn 
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Figure 5. Spectroradiometer spectral re sponse for a “‘daylight”’ 
fluorescent lan p 

Circle how voltage readings on VTVM for the indicated wavelengths for 
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data for 


by the recorder or as would be obtained by 0 
a very large number of voltage readings. Nexi {| 
integrated areas of the triangles representing {| 
mercury emission lines are determined and eval 

in terms of continuous spectrum areas per 5 
centered at the wavelength of the mercury emis 
The continuous spectrum curve is next reduc: 
relative spectral energy through the use of the 
trometer response factors (see fig. 4) and the 
plotted (see figs. 6 to 9) Finally, energy v: 
corresponding to the mercury emission lin 
determined from the relative integrated areas 
correction by the bright line to continuous spect 
ratio are plotted on the curves. They are show 
triangles whose bases are 10 muy, as would be t! 
case with fixed wavelength slits of 5-my spe 
width. 


s 


The use of two phototubes, the one more sens 
to the wavelengths and the other mor 
sensitive to the longer wavelengths, has advantag 
comparable to that of a similar usage in the Beckmar 
spectrophotometer [35, 36]. The apparent shift 
wavelength (see fig. 5) resulting in spectral regions 
wherein the photocell has an increasing (or decreas 
ing) spectral response may cause difficulties if th 
shape of the energy distribution curve of the lam 
under investigation differs too greatly from tl! 
of the standard. The mean between measurement 
with cells of the two types of spectral responses 
should give better data than obtainable with 
single-cell type In the present work 
the shapes of the curves in the spectral region b 
tween about 520 and 640 my is a good example of 
what may be expected when the two types of photo 
cells are used. 

If it is desired to obtain an absolute evaluatio 
of the total spectral radiant energy from the lamp 
at a specified distance and direction from the lam 
it may be accomplished through the use of a filter 
(or combination) having a transmission band withu 
the range of principal spectral emission of the lamp 
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(see fig 
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and the application of the method illustrated in 
table 3 for the tungsten-in-quartz lamp. The 
mercury emission lines are placed within the proper 
wavelength intervals or prorated between intervals 
in those cases where the wavelength of mercury 
emission is not closely coincident with the center 
of the wavelength interval chosen. 


Il. Fluorescent Lamp Emission Spectra 


Fluorescent lamps employed for illuminating pur- 
poses appear to be made of varying quantities of two 
phosphors, one having a peak near 480 my and the 
other near 590 my. Changes in composition or treat- 

: 40 S00 ; “oP ment of the phosphors result in varying emittance in 
ee ee the blue relative to that in the yellow-orange spectral 
spectral energy distribution « region {16 to 20]. Hence, the lamps are thus made 

white’? fluorescent lamp ? o, 99 . >? . 
: - ; ‘““white’’, ‘cool white’’, or “davlight , ete., in general 
———e color character. Special colors may be made by 
further control of the manufacturing process or 
through the use of dyes or colored tubing in the lamp 

manufacture. 

As previously mentioned and illustrated in figures 
5 to 9, fluorescent lamps have an appreciable amount 
of radiant energy within the four visible mercury 
emission lines, 405, 435.8, 546.1, and 578 mp. Hence 
their color is a combination of the two spectra. The 
phosphors employed control the continuous spectrum 


— but have little effect upon the output of the mercury 
ith a | | emission lines. In the illuminating lamps examined 
S | | an appreciable amount of energy was present at 
n be- s | | 365.5 mu. A small amount of radiant energy, often 
ple of ' - ne aes > © | approximating 10 percent of the value at 365.5 was 
hoto- : <a also observ able at 313.2 my The spectral energy dis- 
ative spectral energy distribution of a fluorescent | tributioncurvesdisplayed were obtained on individual 
lalior ee lamps available to the author and are reproduced 
lamp Balthasar ttc sc tii primarily to illustrate the use of the methods and 
lamp equipment employed. 
filter As previously noted, a knowledge of the precise 
vithir ee — spectral energy distribution of a fluorescent lamp is 


lamp | | important for use in the evaluation of its color 
characteristics. Using presently published data on 
the standard observer and evaluating the spectral 
data according to the z, yy, 2 coordinate system for 
colorimetry [31] a preliminary calculation on one 
“cool white’ fluorescent lamp gave values of + and y 
of 0.372 and 0.357 as compared to 0.376 and 0.368 
visually observed by members of the Colorimetry 
Laboratory of this Bureau. While the agreement is 
not exact, it should be considered good in view of 
the difficulties involved in the precise measurement 
of spectral energy distributions of this type and their 
— gig ti proper color evaluation. 
ELENGTH IN MILLIMICRONS The fluorescent sun lamp is different from the 
URE 9. Relative spectral energy distribution of a fluorescent | ordinary illuminating fluorescent lamps only in the 
Black Light lamp facts that its envelope is of a special ultraviolet 
transmitting glass (having a short-wave cut-off near 
280 mu) and that the phosphor employed is one hav- 
ing a fluorescence within the erythemal region (be- 
tween about 300 and 360 my) of the spectrum [21 to 
23). In both cases the fluorescence is excited pri- 
marily by the energy emitted at the 2537-A mercury 
line. Furthermore, in both cases any 2537-A radiant 


443 


oQ- 
v4 











energy not transformed by the phosphor is absorbed 
by the glass envelope 

In addition to the ultraviolet radiant energy 
emitted by sun lamps, the visible mercury emission 
lines are present in approximately the same relative 
intensities as found in the case of the illuminating 
lamps 

Fluorescent lamps may be constructed for any 
special use wherein radiant energy from any limited 
spectral region from about 290 to 800 mu may be use- 
ful. Such a lamp (BL-360) [21] having an appli- 
cation in black light experiments is illustrated in 
figure 8 


IV. Discussion and Conclusions 


The method for evaluation of the relationship 
between the spectroradiometer response for line rela- 
tive to continuous spectra is based entirely upon 
experimental data and permits a precise representa- 
tion of the relative magnitudes of the two components 
in a mixed spectrum as found in a mercury type 
fluorescent lamp. Any peculiar mechanical or opti- 
cal characteristics of a particular double instrument 
are thus taken into account by this method of 
calibration 

The high sensitivity of the amplifying and record- 
ing system employed in this instrument, together 
with its high signal-to-noise-ratio characteristics, is 
made possible through the use of a tuned modula- 
tion and amplifying arrangement. The choice of a 
frequency contaming no harmonics of 60 or 120 
cycles allows its use in the measurement of the 
radiant energy from a fluorescent lamp having almost 
100-percent light modulation at 120 cycles 

The employment of alternating-current amplifica- 
tion permits a stable zero on the output meter or 
recorder, which allows continuous operation through 
the spectrum with no necessity for checking or for 
“cuessing’’ where the 
of the measurements 

In the work herein displayed the 
the data, especially for wavelengths shorter 
about 400mu, depends upon the accepted value for 
the emissivity of tungsten. Above 510my_ the 
mean between a color temperature and a true tem- 
perature standard wherein published emissivity data 
were employed results in fluorescent-lamp spectra 
that are slightly higher (about 1 to 5 %) in the green 
and extreme red relative to the blue than if use had 
been made only of a temperature standard 
This difference is not great but indicates an urgent 
need for an accurate determination of the emissivity 
of tungsten. The high sensitivity and stability of 
the instrument herein described suggest its promising 
usefulness in this field 

Other investigations wherein this instrument is 
adaptable are those requiring precise measurements 
of low reflectivity or of high optical densities, for 
example the spectral transmittance of dark shade 
welding glass. Tests indicate full recorder 
traces for radiant energy through 12 glass 


“zero” was during the course 


aceuracy of 
than 


color 


scale 


shade 


having as a light source a 500-w incandescent lamp 
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about 12 in. from the spectrometer entrance 
Through the use of a condensing lens it shoul 
possible to easily determine the spectral ti 
mittance of a shade 14 welding glass. A prelimi 
study of the ultraviolet spectral distributio: 
radiant energy from the sun [32] indicates wu 
applications of this instrument for investigatior 
total atmospheric ozone. 
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Penetration and Diffusion of X-Rays. 


Calculation of 


Spatial Distributions by Polynomial Expansion’ 


Lewis V. Spencer and U. Fano 


of the spatial distribution of X-rays in suitable systems of polynomials 


rge rapidly 
medi 


im. The 


ed with a view to its possibl 


I. Introduction 


Theoretical work at the National Bureau of 
Standards on the propagation of X-rays through 
large masses of homogeneous materials has so far 
been directed primarily to the study of extremely 
deep penetrations [1, 2, 3] In the complementary 
case of shallow penetration it was thought that a 
straightforward calculation of the effect of successive 
Compton scattering would prove con- 
venient. The intention in this paper is to present 
a different method of calculation that can be carried 
out conveniently for moderate penetrations (up to 
10 to 15 mean free paths of the primary radiation 
and that relies on appropriate application of con- 
ventional phy sico-mathematical techniques 

The propagation of X-rays through an infinite 
medium is governed by the transport equation 


processes 


u-grad \Nir.u,» dn’ k(n’, d 


. 


’ ae ’ 
du’ =— 6 N 


oie 


N(r, 2’ rd’ 


- SOUTCE 


where 

Nir,u,d) is the density of photons of wavelength 
\ (in Compton units), at the point r traveling in the 
direction u per unit solid angle; u(\) Is the total 
narrow beam absorption coefficient; £(\’,A) is the 
Klein-Nishina differential coefficient for Compton 
scattering with a wavelength change from ’ to X: 
and 4(z) is Dirac’s delta-function 

It is widely known from the study of transport 
equations of type (1) that moments of 
the space distribution NV are related by a recurrence 
relationship The zero-th moment, that is, the 
simple space integral of V, has already been calcu- 
lated numerically in a number of sample cases 
[4,5]. A few more moments can be calculated 
numerically without excessive labor by manual 
operation. Automatic computers might prove use- 
ful for the same purpose 


SuCcCeSSIVEeG 


Calculations are performed for various s« 
surprisingly 
use in a variety of problems 


jturce geometries in al 


good convergence of this method of calcu 


It is also widely recognized that a knowledg 
the moments of a distribution function enables on 
in principle to reconstruct the function itself. How- 
ever, little effort seems to have been applied along 
this line in the study of X-ray or neutron diffusior 
problems, presumably because the method was not 
Actually, the m 


expected to prove convenient 
conveniently 


ment method can very 
obtain rapidly convergent expansions of the dist 
bution function in suitable polynomial systems 


be used 


II. General Discussion of Polynomial | 
Expansions ‘| 


Consider, for simplicity, a distribution fun 
f of a single space variable x, which ranges fron 
to ©; f is assumed to vanish sufficient) 

rapidly at both ends of this range To obtai 
polynomial expansion of f, a suitable weight fu 
first that vanishes at 
chooses a succession of polynom Ls 

Pr(a which has an adjoint 
also polynomials) Pos Pri A . , el ni 
such that to 


d di 


tion ws Is chosen 
One then 
Po Pils 


succession 


J pP 


One can then write 


Jn On Pa (2 


The coefficients a, are given by the standard form 


Pp, Ia dr. } (| 


In connection with each function w(zr), one n 
but need not, use the special set of polynomial 


which are identical with their adjoints p 


Knowledge of all the moments i 


{ 


distribution function determines the values o 
the coefficients a, and thereby the complet: 
pansion (3) with the preselected weight functior 


polynomials. If only the first r moments are kn - 
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first 7 1 coefficients may be determined exact- 
the (r+ 1)-th adjomt polynomial is of degree r. 
the function f(r) is the solution of a transport 
on, it may not be necessary to calculate the 
ents first and thence the coefficients a, There 
recurrence relationships between the successive 
which are analogous to the relation- 
between the moments. One may 
calculate the a,’s one after the other directly. 
significance of the expansion (3) of a distribu- 
function into a series of polynomials is illus- 
d by the following argument: The first term of 
series contains a polynomial of crade zero, that 
Therefore the weight function 
tself represents the initial approximation to the 
to within the factor agp If wa 
sen just equal to the unknown function f(x), the 
vansion would reduce to a single term. In gen- 
one can expect the expansion to converge rapidly 
weight function that 
{hie iently good approximation to f(a 
For example, one may think of embodying into 
advance information that may be available on 
behavior of f(x) for very large and very small 
ies of 2 It is a matter of convenience in any 
ticular application how complicated one should 
ke w(2) to obtaim a rapidly convergent expansion. 
the practical application to the X-ray problem, 
choice of the simple exponential function for 
vields convergence that seems quite adequate 

for a surprisingly large range of 2 


heients a 


successive 


he constant dpp 


tion fiz were 


chooses a constitutes a 


The suecessive terms ef the expansion may be 
eribed as corrective terms that are added to the 
tial approximation w(r)=f(r) in order that 
cessive moments of the approximation function 

w(x) |aop 


a Pp J t A2Poe\a 


oincide with successive moments of the actual fune- 
t f(a For example, suppose that an initial ap- 


that 
Then, it 


turn out that the actual value of the first mo- 


has been chosen so 


w(x) da r)dz (i. e. 


roximation wiz Ae 


ad pP ] 


rf(rx) dr is larger than the ‘estimated’’ 


dx. To correct for this discrepancy 


initial approximation should be modified by the 
term w(r)a,p,(2) which in- 
on the side of positive x and reduces it 
negative Actual calculation of high 
r moments should usually be required only if 
has to fit the “‘tails’”’ of the distribution function 
It should not introduce substantial corrections 

e main body of the distribution 


tion of a corrective 
wid 


side 


III. Illustrative Examples 


the La- 
otherwise 
expt 7 
for «>0 


have examined 
two 


a demonstration, we 
re polynomial 
vn distribution 


expansions of 
functions, namely 


yz) and exp(—z),F\(—3/2. 1,—2), 


(Laguerre polynomials are the self-adjoint polyno- 
mials with the weight factor exp r) over the range 
Oto « These distributions represent the solutions 
of schematized X-ray or neutron problems (1, 6]; the 
former pertains to the constant mean free path case 
(u(r const the the which 
u(r a bd 

The values obtained from successive polynomial 
approximations of the “build-up factors” J)(27 
and |F 3/2, 1,—xr) are compared with the exact 
values in tables 1 and 2. The fit is surprisingly 
The addition of one more polynomial to the 
expansion appears to extend the range of satisfactory 
fit by at least an interval 6r~3 into the tail of the 
distribution function in the direction of increasing 7 
The tables include the results of expansions ID 
Laguerre polynomials, as well as expansions in the 
related set of U, polynomials that are introduced 
in the next section 


second to ease mn 


2.3 


cood 


ABLE 1 


presentation J ce r) bya 


erms of the er pansior 


To illustrate further the workings of these expan- 
sions, one may consider the integral representation of 
the distribution functions as “inverse Laplace 
transforms’’. In the case of the Bessel function, we 
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hh 
dN’ k(d’ A) PA — AEN) Nie’) + 8(x) Fd 
The deep penetration of X-rays is controll 
the absorption of their most penetrating 
ponents.[1]. Therefore, one may describe th 
pendence of the photon density on the dept 
penetration z by means of a sum of polynomial! 
/ with the weight factor exp (—az). An ob 
: choice for the parameter a in the weight factor j 
where the L,’s are the normalized Laguerre poly-  yalue of u(X) for the most penetrating compon 
nomials. Thus, the expansion in Laguerre poly- | the radiation. but a 
nomials is obtained by expanding the Laplace | wel} prove better. 
transform of the distribution function into powers of As we mentioned earlier, the ortho-normal s 
(1—t)/t. This expansion converges rapidly near | polynomials corresponding to the weight factor 
t=1; that is, in the same general region where the | (— gr) is the set of Laguerre polynomials give 
integrand exhibits a “saddle point” forz~1. Each | 7 (ar)—y.(n!/(n—r) !p®) (—ar)®.. We accord: 
term of the expansion has a pole at the point t=0, represent the function N,(z, d) - the fesue 
where the Laplace transform is singular. An expan- = 
sion including about four terms affords a good he a — ; 
. ; Ni(2, Yn41n(A)L, (ax), 
approximation of the Laplace transform from t= « 
(saddle point for r=0) to perhaps ¢ 19 (saddle point 


for r=9). ly L,, (ax)N,(x, \)adzr. 





less obvious choice of an 


IV. Calculation of Polynomial Expansions The equation that determines @:,(X) is foun: 


The procedure for obtaining a convenient poly- multiplying eq 5 by aL, (ax) and integrating ov 
nomial expansion of the distribution of X-ray from 0 to ~. rhe integration of the left side of (8 
photons varies somewhat depending upon the | 'S carried out by parts. he values of N;,, and A 
geometry of the source. Therefore. a number of | .2t the limits 0 and © vanish as shown by the follow 
different cases will be treated separately. ing arguments: 

1. No radiation penetrates to r= @ ; 
1. Plane Monodirectional Source: High-Energy 2. The density of incoming photons at the entran 
Photons Only surface should be taken equal to zero if the sow 
term 6 (x) is included in the range of integratio1 

The simplest case, perhaps, is that of a parallel 3. There are no photons going back out of the 
beam of X-ravs of infinite extension that enters a | ©Dtrance surface, as the treatment is limited to hig! 
material perpendicularly through the plane z=0 photon energies. 

The presence of the entrance surface does not cause The following identity is used in carrying out th 
any difficulty in this case because no high energy | remaining integration, which involves the derivatiy 
secondary X-rays will be present at the boundary. of the Laguerre polynomial 

In this case the photon density N (r, u, \) depends 
only on the components z and u, of the vectors r dL,(y 
and u. The dependence of N upon the direction dy 
cosine u, may be expressed by means of an expansion 
in Legendre polynomials The integrals of all terms of (8) are thus reduced | 

the form (9a). The result is 
Ni(z,d) Pi(uz) 


a 


9/4 pao lds L)Qja1 2 la, al - B(A) Gz, (A) 


The incident beam is represented in eq 1 by the Pa. : ; 
AN kN NYP, A+A')din(A) +af(A). 


source term 
4 
Sounve These equations coincide, for n=0, with t 
whose numerical solution has been discussed 
where /(A) is the incident spectrum. previous paper [5]. The equations for n>0 ca 
Multiplication of eq 1 by P, (u,), followed by in- solved in succession by the same method, b« 
tegration over u yields: 
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aS ae 


of the Le ge ndre -Laq erre 


ee 


link ages coe ficients 


sum on the left side contains only terms with 


lt is often of interest to know the density of pho- 

ns at a point irrespective of their direction. This 
iensity is represented by the Legendre coefficient 
V,(2,d). Calculation of the coefficients a, (A) of 
the Laguerre expansion of Ny, requires, according to 
10), the calculation of a set of coefficients a,, with 

0. Figure 1 shows diagrammatically how, for 
nstance, the solution of the equation which deter- 
mines @)(A) requires a prior knowledge of the func- 
LIONS Go, Aon, Aq, ANd Gyo. In turn, the determina- 
tion of a, and a2; requires a prior knowledge of a, 
and a If one calculates the moments 


A 
llow - 


ran 
oures 
On; 

f the 
high ' ad 


0 
e 


"N,(z,A) dz, 


it the 
ative 


nstead of calculating directly the Laguerre coefli- 
lents @,,, one finds the simpler scheme of interde- 
pendence illustrated by figure 2. Here, solution of 
the equation for 6,, requires prior knowledge of 
h , and b,, 1; thus, to calculate doo, @o:, Goo, 
,one must also calculate @9, @29, Ago, @1,;, G21, 
2; that is, one must solve numerically a total 
This does not involve an excessive 


and @ 

ind a 

of 10 equations 

amount of labor. 
lt may be interesting to relate the representation 

of the photon density by means of Laguerre poly- 
mials to its representation as an inverse Laplace 

sform: 


l 
9 


2ri 


dip, X dp. 


moments of the distribution N, in depth cor- 
ond to the coefficients of the Taylor expansion 
:(p, \) in powers of p [7]. The expansion of N 
aguerre polynomials is obtained by expanding 
nto a series of terms of type a@,p"/(a—p)"*'. 
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_ a Chain linkage s of the spatial moments of the Le 
cor fhicre nis, 


sets of boxes are not interlinked 


Each of these terms has a pole at p=a. Notice 
that ¢, itself has a singularity at p=4y,, where uy, is 
the smallest value of uw (A) [2]. Thus, the expansion 
of the Laplace transform in powers of p/(a—p) may 
be expected to converge more rapidly than the ex- 
pansion in powers of p if @ is given a value near 


that of py. 


2. Plane Isotropic Source 


This in which radiation is emitted in all 
directions from all points of a plane, is familiar in the 
study of neutron penetration and diffusion [6, 7]. 

The distribution from a plane source, isotropic as 
well as monodirectional, depends only on the coordi- 
nate z of each point r and on the component u, of the 
direction of flow u. The isotropic source term to be 
entered in eq 1 takes the form: 


case, 


6(r) TA) 


source 
47 


Correspondingly, one finds, instead of eq 8 


l 
sei | (+) 


*A 
T dn’'k NV NP g T A)N2(z,d’) rT d(x) f(A d10 12 


ON, 


; ON:-1  / 
+-/ ] — u(rA)N,(x,r 


Ou J 


If the source is surrounded by an infinitely ex- 
tended medium, the radiation is distributed sym- 
metrically with respect to the source plane z=0 


N(x ,uz,A)=N(—2z,—uz,d); Ni(z.d)=(—1)' Ni 
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In this case the moments of the space distribution | A number of these polynomials are given in ta 
involve an integration over the entire range 

: I : Owing to the symmetry conditions TaBLe 3. The first four U, and V, polynomia 
(13), if 7 is even, all the odd moments of N;, are zero 
and if / is odd, all the even moments of N; are zero 

In the preceding case the nth Laguerre polynomial 
added to the expansion reflected, in essence, informa- 
tion gathered from the calculation of the n-th moment 
of the distribution function. In the present case 
every other moment vanishes; that is, does not 
contribute to our knowledge of the distribution 
function; therefore, the method applied in the 
preceding case requires some modification. The n-th 
polynominal should now reflect information provided 
by the value of the 2n-th or 2n+1-th moment. An aa : 
appropriate set of polynomials is required for this lo determine coefficients of the expansior 
one may utilize the sets of adjoint polynomial 
and V;, which fulfill the orthogonality conditi: 
together with Ul, and V,, respectively. 1 
polynomials are 


purpose 

The weight function should be selected in such a 
way that the distribution function has the correct 
symmetry. To maintain an analogy with the case of 
a monodirectional source, one may choose the weight ; 
function exp (—az') for the expansion of the Ss 
even-/ Legendre coefficients V r, XA). Correspond- : — ) 
ingly, for odd values of 7 the weight function 
rexp(—ar|) is perhapsthesimplestchoice. We write 


> a r)l , l even We multiply eq 12 by al} (az), if 1 is even 
aVy (ax), if/ is odd, and integrate over 2 
to o. The integration proceeds as the int 
> a,.(r) V. ’ l odd tion of (8). The left side integrated by parts 4 
n=O one term that vanishes and an integral that cont 
the derivative dl dx (Leven) or dVy/dz (1 « 
where LU’, and V, are polynominals of degree n still | These derivatives fulfill the identities 
to be determined. 
As indicated in section IV, 1, the first n terms of dl y 
each expansion should embody all the information 
that is provided by the values of the first n usable 
moments of N r,X In other words, the higher d\ Y 
terms of the expansion (14) should not contribute to dy 
the lower usable moments. In the case of even /, 
this implies that the component exp (—al\zr|)U, (aa Thereby, the integration of eq 12 reduce 
of the function N, (vy, \) has vanishing low-order | grals of the type 
even moments 


dy 


for n’ 





These conditions determine the sets U, and V, to 

within a normalization constant. We find the fol- 

low ing expressions 

lo ler 


vely 


2 ay"? 


7 2! en 
, . (16) | ar 


Higher associated polynon 
| nomials analogous to those o 
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is again the 
The result is 


/ odd 
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20 


These equations can be solved chainwise, much like 


the eq 10. Notice that a 0 for n<(l/2—1 

since the source term is confined to the equation 
with /=0. The equations must be solved in the 
following order: doo, @yo, @21, Agi, * * +, Apr, @y}, doy + * *, 
Gea, @ (og +++. Here again, as in the case 


of the monodirectional source, the calculation of four 


terms of the expansion of N», (x, \) requires the 


solution of 10 equations. 


0 for nS (//2—1) relates to the 


The fact that a, 


known fact that the moment }, r"N,(z,d) da 


vanishes when n</ [7]. (Only the moments corre- 


sponding to the shaded boxes in fig. 3 are not zero. 
The value of each moment 6,, depends on the val- 
s of the moments 6;,;,—-;, and 5, as shown 
fig. 2. The calculation must begin with dg, 
hose equation is inhomogeneous, that is from the 
per left corner of the scheme.) Accordingly, the 
iplace transform of N/,(z,\), namely, ¢;(p,), has 
Taylor expansion that begins with the /-th term. 
he polynomial expansions (14) of N;, are obtained 
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n’-th spatial moment of 


as inverse Laplace transforms if ¢ is represented as 
a series of terms p*"/(a*— p*) for even /, or 
pr" a’ — p*)"*' for odd 1. 


n+l 


3. Plane Monodirectional Source No Energy 
Limitation 


The method developed for the treatment of X-rays 
from a plane isotropic source can be applied to the 
study of radiation from other plane sources. To 
this end, one must generalize the preceding treat- 
ment somewhat. 

It will be remembered that because the source is 
symmetric and because of the interlocking of mo- 
ments illustrated in figure 2, the only useful moments 
are those which correspond to the black squares of 
the checkerboard pattern of figure 2. The half- 
checkerboard pattern of useful moments (fig. 3 
results from the isotropy of the source (eg 11 In 
the general symmetric case, source terms, will be 
present in all equations, not just those alone in 
which /=0; and the full checkerboard scheme of 
moments (fig. 2) will be used. 

On the other hand, if the source were antisym- 
metric in z, the moments corresponding to the 
shaded squares in figure 2 would all be zero. The 
useful moments would be interrelated among them- 
selves just as before and would correspond to the 
white squares of figure 2. 

One can take advantage of this separate inter- 
locking of different sets of moments by regarding a 
plane monodirectional source as the superposition 
of a symmetrical and an antisymmetrical source. 
The symmetrical source component injects equal 
photon densities in opposite directions perpendicular 
to the plane z=0 on either side. The antisym- 
metrical source injects a photon density equal to 
that of the symmetrical source of the positive side 
of r=0 and a “negative” photon density of equal 
size on the negative side. 

The photon distribution N (x,u,,d 
antisymmetric has symmetry 
opposite to (14), namely, 


due to the 


source properties 


N 


(2 ] 


The polynomial expansions of the Legendre coeffi- 
cients N° (x, \) involve the V, polynomials for even 
l and the U, polynomials for odd /, that is, they 
conform to (14) with opposite even-odd condition. 
The total photon density from a plane mono- 
directional source is thus given by the sum of two 
terms 
N(z, Uz, )=N™ (zx, uz, +N (az, uz, d). 22 
N™ is determined by a set of equations identical with 
(20), except that the source term does not include 
the factor do. N® is determined by a similar 
system with even and odd / interchanged. 
The total density N(x, u,, \) all but cancels on 
the negative side of the source plane, inasmuch as 








N is negative there. The residual photon density 
for z<_0 represents that has been back- 
scattered behind the should of 


comparatively low energy photons. 


radiation 


source. It consist 


4. Point Isotropic Source 


The photon density generated by a point isotropi« 
source can be deriv ed from the solut ion corresponding 
to a plane isotropic source (¢ This 
derivation is especially simple if one calculates only 
the total density of photons at each point irrespective 
of their direction For this reason we shall outline 
here only sketchily an independent calculation for 
the case of a point source 

Owing to the spherical symmetry of the problem 
the distribution function N(r,u,) depends only on 
the magnitude r of r and on the component u,=u-r/7 
of u, if the origin of the coordinates is laid at the 
The source term in eq 1 the form 
trr*? and the transport eq 1 becomes 


‘use (b 7]. 


source. takes 
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The density N can be expanded in Legendre poly- 
nomials of u,, as in The transport equation 
takes then a form corresponding to (8), whose deri- 
vation makes use of the identity 


du, P 


eq 6. 


L(l+-1 


f(r) 6; 
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In this case we do not derive directly the coeffi- 
cients of an expansion of N,(r,\) in polynomials fn r, 


but we calculate the moments 


r"® N;, rx tarr?d) b,, (A). 


/V 


To this end one multiplies (24) by r" and integrates 
| : { 


over all the space. The result is 


lL d+nd- L(l+n-+-) b, 


yi b, 
2l+-1 
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— p(d)b,,(A)4 dN k(N YN) PO 


+ 1(A) bid, 


This set of equations can be solved chainw 
exactly the same way as in the plane isotropic 
Coefficients of the U, and V, polynomials ¢; 
obtained by adding and subtracting moments 
coefficients in the expansion of the total dens 
photons in LU’, polynomials are related in the fi 
ing way with the corresponding coefficients i: 
plane isotropic source case 

on) point 2n+1)(a n plane 27 a 
The coefficients a,, for the two sources in th 
140 are related in a more complicated way. 


5. Point Monodirectional Source 


This kind of source may properly be regarded as 
the elementary source, since any source can be 
represented as an aggregate of point monodirec- 
tional sources. The case of a point monodirectio 
source is more complicated than the cases considered 
above, because it involves explicitly the diffusion of 
radiation sidewise from the initial “line of fire.’ 

Cylindrical coordinates r=(z, p, ), with th 
axis lying along the line of fire are appropriate to th 
geometry of this problem. The direction vector u 
can be described by polar coordinates (0, db+-¢4 
that the zero azimuth ¢=0 refers to a plane throug! 
the axis and the point r under consideration 
Owing to the cylindrical symmetry of the source th 
photon density depends explicitly on the variables 
in the following way 


Ni(r,u,> >, p, 6,®, X). 

To calculate the expression u-grad N explicitly 
one must that the gradient is taken for 
constant U, means 6 Accordingly 


consider 
which 60 
ON 
+sin cos ® 


Op 


oO 
cos @ 


u-grad N 


sin @sin ® ON 
p Ob 
In this problem the distribution function should be 
expanded into Laplace spherical harmonics Y,,,(@, 4 
and moments should be taken with respect to both 
space coordinates z and p. Therefore we introduce 
the set of dependent variables 


: | p? ‘dp e"dz | “sin ad @ 


NrP(X) 


im 


- db Y*,(0,b) N(z,p, 0,8,X) 


A corresponding four-fold integration applied t 
| 1 reduces this equation to the form 
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[his reduces to the equation for the 
nlane monodirectional source, as it should do, when 
p=m=0. The coefficients Ne@(d) can be calcu- 
lated accordingly. Once this is done one can pro- 
ceed to p=1, m 1, and then on again. 

In the preceding cases the calculation of a three- 
variable function such as N (r, w,, 4) was reduced to 
the calculation of a two-dimensional array of co- 
efficients N?(A). The present problem involves 
five independent variables and reduces to the solu- 
ion of an “array of problems,” each of them as 
complicated as any of the previous ones and corres- 
ponding to one pair of values of p and m. The 
genetic connections in the array of problems are 
represented in figure 4 

Various polynomial representations of the radial 
distribution might suggest themselves in the course 
of any practical application. The use of a gaussian 
weight function should presumably be convenient 
at least in the lower energy range. 


equation 


V. Numerical Applications* 
1. Penetration of 10-Mev Gamma Rays in Lead 


lhe method of section IV, 1 was applied numeri- 
cally to the case of a parallel beam of 10.2-Mev 
photons entering an infinitely thick layer of lead. A 
heavy element was chosen because it is a more diffi- 
cult situation to treat with asymptotic methods, 
wing to the minimum of the absorption coefficient 
3 Mev [2]. The same initial energy was 
hosen which is used by Karr and Hurwitz in their 
asymptotic calculation with the “straight ahead’’ 
approximation [8]. The calculation was extended 
only down to photon energies of about 1 Mev 
1.49), partly because of the limitation to the 
lity of the method indicated in section IV, 1. 


litional applications were reported by L. V. Spencer and F. A. Stinson 
New York Meeting of the American Physical Society, February 3, 1951 


This calculation was carried out not ony by the 
method of section IV, 1 but also with a correspond- 
ing ‘“‘straight ahead” approximation in order to esti- 
mate the error involved in neglecting angular 
effects. The straight ahead approximation implies 
disregarding the effects of deflections experienced by 
high-energy photons in the course of Compton scat- 
tering. In this approximation one solves only the 
eq 8 for /=0, replacing ON,/Or by ON)/dzr, so that 
No (x, ») is the only distribution function to be 
determined. 

Finally, the straight ahead calculation was car- 
ried out once more, neglecting the sin*@ terms in 
the Compton cross section, that is, omitting the 
last two terms of the expression 


“2rd 
) +——2A—A h—X}* lence 
es 


The purpose of this calculation was partly to esti- 
mate the importance of these terms and partly to 
obtain a proper comparison with the work of Karr 
and Hurwitz, who used this approximation in mak- 
ing their calculation. 

Four terms of the polynomial representation (9) 
of the function No(z, \) were worked out with a 
set equal to u,,, the minimum value of y(d). For 
comparison, the coefficients were regrouped to cor- 
respond to a polynomial expansion with a set equal 
to u(A=0.05) and also to a polynomial expansion 
with the weight factor z exp (—y,2x), the last being 
perhaps the most realistic weight factor in view of 
what is known about the initial and asymptotic 
behavior of the build-up factor. The results using 
these different weight factors were all in satisfactory 
agreement. 

Figure 5 shows the differential spectrum given by 
calculation according to section IV, 1. The energy 
density is plotted as a function of the energy for 
various penetrations. Notice the accumulation of 
energy in the range~3 Mev, corresponding to the 
minimum absorption coefficient. 
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Figure 6 shows the angular calculation compared 
with the two straight ahead calculations. The energy 
density of the component \=0.49 is plotted against 
depth of Karr and Hurwitz carried 
their asympotitic calculations down to u,2=10. It 
s possible to compare their results with ours at that 
penetration. The circle gives their value and is to 
be compared with our upper curve. The two values 
are much closer together than was expected from the 
accuracy of the calculations. Notice that the angular 
calculation differs by a factor of almost 3 from the 
Karr-Hurwitz calculation at 10 mean free paths of 
the most penetrating component 


penetration. 


2. Co” Isotropic Point Source in Water 


This problem has been investigated experimentally 


White {9}. We have made numerical calcula- 
tions by using the method outlined in section IV, 4. 
he coefficients of the first 4 (’, polynomials in 
he expansion of No(r,A) have been evaluated. 
Below 0.150 Mev, the coefficients of (, and U’, were 
imated from the trend of the lower coefficients, 
ce approximate equilibrium was attained.) 
Some results are illustrated in figures 7, 8, and 9 
rure 7 shows the differential energy spectrum at 
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mmental results of 


(Comparison of the cali 


facto pith the expe 


various depths of penetration. The same results are 
plotted semilogarithmically in figure S to illustrate 
the approach to an equilibrium spectrum at the 
lower energies. At equilibrium the curves become 
parallel, indicating that the spectrum no 
changes with increasing penetration. 

Figure 9 compares White’s experimental build- 
up factors, as measured by an tonization cham- 
ber (and a counter), with a theoretical build-up 
factor. The latter was evaluated by weighting the 
calculated spectra with the theoretical response in 
roentgen of a standard free-air ionization chambet 
and then again with the specific correction factor of 
the Keleket III pocket dosimeter 110] The integral 
response over all spectral components was divided 
by the calculated response to the primary radiation 
to give the theoretical build-up factor. The slight 
downward curvature at the end of the theoretical 
curve is an artifact that indicates that the 4-term 
approximation is beginning to fail at this distance. 


longer 


VI. Discussion 


The analytical developments and the numerical 
applications presented in this paper seem to give 
some confidence in the following conclusions 

a. A method is available for calculating the dis- 
tribution of X-rays in uniform media with a moderate 
amount of labor up to fairly large depths of penetra- 
tion. 

b. This technique can presumably be applied with 
success to the broad class of transport phenomena 
governed by Boltzmann-type equations. It can 
take advantage even of rather crude theoretical 
predictions on the behavior of a distribution function 


to formulate an initial approximation that can 


455 








thereafter be rapidly improved by straightforward | treatment which will be the object of a sep 


numerical work report 
Nevertheless, much further practical experience Continuation of this work will aim primarily i 
may be required before one can assess with some | directions: 
confidence the criteria for the most economical a. To evaluate the effect of boundaries bet 
application of this technique. For example, it is not | two media and, in general, to consider the proj 
immediately apparent how far one can conveniently | tion through inhomogeneous barriers 
push the search for an improved weight function b. To include the effect of X-ray regeneratix 
The present method should probably prove con- | the cascade shower mechanism, which be 


venient for mapping out the distribution of various | important at high energies 
radiations within thick layers of scattering and 
absorbing materials. The necessary formulae for 


calculating the spac ial moments of electron distribu- VIL. References 
tions have recently been developed by Lewis [11] 
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A Study of Fatigue in Metals by Means of X-Ray Strain 
Measurement 


John A. Bennett 


X-ray determinations of lattice strain were made on steel specimens under a series of 


static bending moments. 


normalized and after various amounts of fatigue stressing 


These series of determinations were made on the specimens as 


The slope of the line relating 


the bending moment to the lattice strain was found to decrease with increasing numbers of 


cycles of alternating stress. 


However, the decrease was of approximately the same mag- 


nitude for specimens stressed above and below the fatigue limit, showing that the change 


was not associated with fatigue damage. 


I. Introduction 


In 1941 Glocker and coworkers! in Germany, 
wublished a report on the study of fatigue damage 
»y means of X-ray stress measurement. He claimed 
that as fatigue damage increased there was a change 
in the value of the surface stress measured with a 
given static load on the specimen. This claim, if 
verified, would have required a fundamental change 
in the prevalent ideas regarding the mechanism of 
fatigue, according to which fatigue damage occurs 
in a very small volume of the metal and the rest of 
the material remains unchanged. This point of 
view would have to be abandoned, if it were possible 
to detect fatigue damage by means of X-ray diffrac- 
tion, as the diffraction occurs from a relatively large 
volume close to the surface of the specimen. 


} 
t 


Because certain features of Glocker’s work made 
the results open to question, a project was initiated 
with a view to reproducing his experiments. It was 
immediately apparent at that time that in order to 
obtain an unambiguous result the sensitivity of 
stress measurement would have to be improved over 
that which could be realized with the methods 
available. The techniques and instruments that 
have been developed in achieving this objective have 
been described previously.* The improvement in 
sensitivity and accuracy was made possible by 
increasing the number of measurements made on each 
diffraction ring and by the use of instrumental rather 
than visual methods for making the measurements. 


II. Materials and Test Methods 


The specimens used in this work were made from 
SAE 4130 steel having the following composition 
n percent): C, 0.33; Mn, 0.53; P, 0.013; S, 0.019; 
Si, 0.24; Cr, 0.91; Mo, 0.21. This steel was from 
the heat used for tests reported previously,‘ being 
identified in that report as material M357. The 


JR. Glocker, W. Lutz, and O. Schaaber, Evidence on the state of faticue of 
ls obtained from an examination of the surface stresses by means of X-rays 
r. deut. Ing. 85, 793 (Oct. 1941 
A. Bennett and H. C. Vacher, Calibration of X-ray measurement of strain 
search N BS 40, 285 (1048) RP1874 
A. Bennett, Instrument for measurement of X-ray diffraction patterns, 
Sei. Tnatr. 20, GOR (1949) 
A. Bennett, A study of the damaging effect of fatigue stressing on X4130 
Proc. Am. Soc. Testing Mat. 46, 693 (1946); J. Research N BS $7, 123 (194 


grain size was small, corresponding to No. 10 on 
the ASTM scale, which made the steel particularly 
suited to X-ray diffraction work. 

The specimens were of the commonly used R. R. 
Moore type having a minimum diameter of 0.25 in 
and were heat-treated after machining and polishing 
in order to leave the surface free from plastic de- 
formation. The polishing was done manually with 
Aloxite paper, the direction of polishing being 
parallel to the length of the specimen. After they 
were polished, the specimens were heated in vacuum 
to 1,650° F and held 1 hr. They were then cooled 
by removing the furnace from around the vacuum 
chamber, resulting in a cooling rate nearly as rapid 
as air cooling. This treatment will be referred to as 
normalizing in the remainder of this report. Finally, 
the reduced section was etched 2 min in 1-percent 
nital, this having been found necessary to obtain a 
sharp diffraction pattern. 

The procedure used in all the tests reported here 
was to measure the lattice strain under a series of 
applied static loads. These series of determinations 
were made on the specimen in the normalized condi- 
tion and after various amounts of fatigue stressing. 
The strain determinations were based on the changes 
in the diffraction angle between patterns taken with 
the X-ray beam incident at two different angles 

The jig for applying the static load to the fatigue 
specimen is shown in figure 1. By rotating the load- 
ing screw, either a tensile or compressive stress could 
be applied to the front surface of the specimen. 
The load was applied through two calibrated springs 
(only one of which can be seen in the figure) and a 
pair of yokes that distributed the load symmetrically 
to the two pins in the specimen adapters. The de- 
flections of the springs were measured with a dial 
indicator having a long projection, which was in- 
serted through a hole in the shank of the spring. 

The springs were calibrated by dead-weight load- 
ing and the bending moment on the specimen was 
calculated from this calibration and the measured 
distance between the pairs of pins in the specimen 
adapters. Near the end of the investigation reported 
here, this calibration was checked by means of a 
specimen having a rectangular cross section between 
the tapers. Wire strain gages were mounted on this 
specimen and the bending moment applied to the 
specimen in the jig was measured directly. This 
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that scans the diffraction ring ts an are having 
same radius of curvature as the ring. 
For each stress determination two patterns w 
exposed, one with the X-ray beam incident nor! 
mome : in good agreement with that | to the specimen and the other with an incident ang 


bending 
ing deflection of 45 In the latter case the incident beam 


calculated from the spi 


The | ly was mounted on the assembly de- | axis of the specimen and the normal to its surface 


it ani 
ene 


n footnote 2, which provided | are in the same plane. Twelve measurements 
ingle between the specimen | ring radius were made on each pattern; for the 90 

for traversing the specimen | incidence patterns the ring radius would be expe¢ 

Che film was not oscillated to be nearly constant around the ring, so the 

iring the diffraction angle | readings were averaged For the oblique incidet 
ail he references given in foot- | patterns, however, the radius would theoretical! 

so only a brief outline will be given | vary around the ring except for conditions of 

following svmbols and definitions are | stress. However, by plotting the measurements o1 
these patterns against strain ratio, the points woul 
be expected to fall on a straight line so the intercept 
at 8=0 as determined by least squares can be taker 
as the value representative of the 12 readings. This 
method is applicable only in the case of uniaxia 


| 


radius of the Ka, diffraction ring 
specimen-to-film distance 

diffraction ang 

tan 20 RUS 

angle around the diffraction ring 
diameter of the calibrating ring 

strain in a given direction 

strain in a direction parallel to the stress 
Stram ratio= ee 


stress, 

Since the publication of the reference cited i 
footnote 2, it has been found ieee ssary to mak 
correction for misalinement between the X-ray beam 
and the center of the calibrating ring. The magni- 
tude of this correction is determined from th 

The quantity tan 2@ was determined by measur- | deviations of individual readings from the averag 
ing R and S. The measurement of S was accom- | for patterns taken with the beam normally incident 
plished by means of a micrometer mounted on a | Figure 3 shows a plot of these deviations against th 
separate camera track, the whole specimen mount | angular position on the diffraction ring. The va 
assembly being placed on this track before and after | plotted at each angular position is the average dey 
making an exposure. The radius, R, was deter- | tion observed on six patterns. The assumptior 
mined by measurements on the diffraction patterns; | made that the diffraction ring is a cirele and that 
a calibrating ring was exposed on the film immedi- | observed deviations are due to a displacement 
ately after the diffraction ring, and the distance | tween the center of this ring and the center of 
between the two measured. This measurement, cor- | calibrating ring. The correction would then b 
rected for film shrinkage, was subtracted from the | sine curve, and such a curve has been drawn 
known inside diameter of the calibrating ring to | figure 3, the phase and amplitude having b 
obtain the radius of the Aa, ring adjusted visually to give the best apparent fit to 

A sketch of the instrument used for making these | observed points 
measurements is shown in figure 2 and is fully de- The correction mentioned above was made ra 
scribed in the reference cited in footnote 3. It is | than trying to eliminate the misalinement becaus 
essentially a densitometer having a dial indicator to | the small size of the error. It will be noted that 
measure the movement of the film. The slit in the | data shown in figure 3, which are typical, indica 
optical system is curved so that the beam of light | displacement (between the center of the diffrac 
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120 180 240 
Angle 8 , degrees 


f andimdua 


300 


readings from average 0 


faken at normal incidence, 


the center of the calibrating ring) of less 
mm. The collimator for the X-ray beam 
has pinholes with a diameter of 1 mm, so it appeared 
probable that one of the causes of misalinement was 
ieven illumination from the target, which would 
ause the center of gravity of the beam to be slightly 
splaced from the center of the pinholes. This was 
orroborated by the fact that the deviations changed 
vhen it was necessary to change the X-ray tube. As 
error was inherent in the X-ray 
ipparatus it was not thought worth while to attempt 
to reduce the other factors that contributed to the 
total misalinement 
his correction, of course, made no difference in 
values of the radius obtained from 
the patterns taken at normal 
However, there was a significant change 
the intercept obtained from the measurements on 
oblique incidence patterns (when the corrections 
biained from the normal incidence patterns were 
ipplied Figure 4 the measurements, 
ted for misalinement, from a typical pattern taken 
15° incidence. The points are plotted on a 
form strain ratio scale as shown at the top of the 
figure so that they would be expected to lie on a 
straight line, and it can be seen that this expectation 
pproximately fulfilled 


and 
than 0.1 


this source ot 


the av erage 
measurements on 


dence 


shows cor- 


III. Results 


he above prom edure, then, provides a single value 


in 20 representative of each pattern. Figure 5 
plotted against the bending 
lent applied to the specimen for a typical series 
eterminations. <A positive value was assigned to 
bending moment when the irradiated surface of 
specimen was in tension, and vice versa. The 
r of applying the loads was as follows: 0, 15, 


vs these values 








ncident ongie- 90° 
Incident angie- 45 





20, 33, —35,0in.-lb. In general, a similar schedule 
was followed in all the series of strain determinations 
reported here; the maximum load in either tension or 
compression was of such magnitude as to cause a 
nominal stress of approximately 30,000 Ib in. * in the 
outer fibers of the specimen, the intermediate loads 
being about half as large 

In order to give an idea of the relationship between 
the ordinates on figure 5 and actual dimensions on 
the diffraction patterns, the change in tan 
sponding to 0.1-mm displacement on the pattern ts 
indicated at the right of the figure. On the left side 
of the graph is a curve (G) representing the galva- 


2A corre- 


nometer deflections in the pattern-measuring instru- 
noment which were observed in traversing a typical 
Ka, line. The zero point for this curve was taken as 
the minimum density observed between the Aa, and 
Kee lines, and this minimum is represented by the 
left edge of the graph That is, the total height of 
the graph represents a change of tan 2@ considerably 
less than the width of the diffraction line 
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had shown that the 
of the lattice strain was 


Pri vious work footnote 
most convenient measure 
the difference in the values of tan 2@ between the 
patterns taken with the beam incident at 45° and 
those taken at normal incidence This difference 
tan 26, tan 20,) is shown in figure 6 plotted 
against the bending moment applied to the specimen 
From a study of a number of graphs of this type, it 
was possible to determine which changes from 
linearity in the experimental points were due to 
plastic deformation in the specimen and which were 
simply scatter. The numbers adjacent to the points 
represent the order in which the respective loads were 
applied. The location of the points and the reasons 
for drawing the lines in the particular manner shown 
can be explained as follows: for points 2 and 3 the 
applied loads caused only elastic deflection of the 
specimen, so the lattice strain was proportional to 
the bending moment in this range. However, when 
the moment was increased above a certain point, 
plastic deformation began to take place on one or 
both surfaces of the specimen, and there was little 
or no further increase in lattice strain, therefore 
point 4 is offset to the right of the line joining points 
1, 2, and 3. When the high positive moment was 
reversed, it would be expected that the locus of the 
plotted values would follow the elastic relationship 
shown by the dotted line. On reaching a negative 
lattice strain of about the same magnitude as that 
at which plastic deformation occurred in the positive 
direction, the increase in elastic strain again practi- 
cally stopped so that point 5 is offset to the opposite 
side of the original elastic line. On release of this 
moment, the values again followed an elastic line 


parallel to the original one, the final point 6, repre- 


senting a residual tension in the surface of the 
specimen under study. 
The above appears to be a logical interpretation 


of the observed data, and the practice was adopted 
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Ficure 6. Difference hetween the 90° and 44 potnts shown in 


igure 


The subscripts A and B ind of tan 20 obtained with 90° and 45 


fncidence " 


respectively 


of determining the slope of the elastic line for a 
condition by averaging the slope between po 
and 3 with that between 5 and 6. On return 
figure 5 and drawing lines determined in this w 
both sets of data, it was found that the deviat 
the points from the lines was reasonably small 
other methods of analyzing the data that 
investigated resulted in deviations that were 
than would be expected in view of the precis 
the data. 

This average slope appears to be the charact 
of the data obtained in this investigation that 
nearly corresponds to the composite slope 
Glocker measured. Accordingly, in order to 
the correctness of his hypothesis, the value of 
slope was plotted against the number of stress « 
to which the specimen had been subjected wh« 
series of strain determinations were made. F 
7 and 8 show these graphs for specimens str 
above and below the fatigue limit, respectively 
ordinate being expressed in terms of the perce: 
of the slope for the normalized material. The v: 
of the slopes for these two specimens are give! 
table 1. 
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Slope of 
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curve 


s amplitude 54,000 


s amplitude 


tween 


ts the relationship between the applied bending 
1s measured by the difference in tan 26 for 


\ subsequent test, not shown on figure 8, 
performed on the specimen that had run 51 million 
»/in.2 The stress was raised to 
4.000 Ib/in and the specimen run 25,000 cycles 
before maku ya series of strain determinations. The 
slope obtained from this series was 78 percent of 
This is above the lowest value 
ressing below the fatigue limit, 
d so close to the value after 5110 cycles that 


eveles at 35,000 Il 


the original value. 
btained during st 


change can h 


ardly be considered 


signific 


normal 


pecimen between 35,000 and 40,000 cycles 


was 


‘ant. 


fatigue crack formed in this specimen between 


5.000 and 40,000 e 


veles. 


IV. Discussion 


is believed that figures 7 and 8 provide the 
er to the original problem on which this work 


as based; it can 


be seen that there is 


a decrease 


slope for both conditions, and the magnitude of 
decrease is not much greater for the specimen 
sed above the fatigue limit than for that stressed 


The change, 


to fatigue damage. 


therefore, cannot be said to be 
This is further verified by 


ibsequent test on the specimen stressed below 


aticue limit 


The conditions that 


caused 


the 


greatest decrease for the specimen stressed originally 
above the fatigue limit caused a decrease of only 
5 percent in the slope for the specimen that had had 
prior fatigue stressing. 

While the primary problem for which the project 
was initiated has been solved, this work has raised 
several new questions and has pointed the way to 
other possibly fruitful lines of research. For example, 
it is surprising that a specimen that has been sub- 
jected to millions of cycles of stress with a maximum 
value of 35,000 lb/in (the data of figs. 5 and 6 
represent this condition) can be deformed plasti- 
cally by a static stress less than 30,000 lb/in.2 The 
further investigation of this anomaly would undoubt- 
edly provide valuable information on the mecha- 
nisms of plastic deformation and fatigue. 


V. Summary 


1. Determinations were made of the relationship 
between the bending moment applied to the specimen 
and the lattice strain in the most highly stressed 
fibers as obtained by X-ray measurements. The 
slope of the curve representing these measurements 
usually decreased with increasing amounts of fatigue 
stressing applied before the determinations were 
made 

2. The greatest decrease in slope was of the same 
order of magnitude for specimens stressed above 
and below the fatigue limit. A specimen that had 
run 51X10° cycles at a stress below the fatigue 
limit was subsequently stressed for about 70 percent 
of its life at a stress amplitude above the fatigue 
limit. Almost no change in the slope was found 
after this damaging treatment. 

3. These determinations, therefore, do not appear 
to give any correlation with fatigue damage. 

1. The data indicate that even after millions of 
cycles of fatigue stressing plastic deformation 
occurred under static stress less than the prior 
fatigue stress 


The author is indebted to Lura F. Roehl, who 
made the diffraction patterns and carried out the 
extensive measurements and calculations involved 
in this study He also wishes to thank H. C. Vacher 
for many helpful suggestions in regard to the X-ray 
techniques used 
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On the Mean Duration of Random Walks 


Wolfgang Wasow 


r} mean duration of a discrete 
probability functions 


general trat The 


wil icl t« r 
This lea t x cit inequalities for the mear 


I. Introduction 


This paper is concerned with random walks in a 
bounded open domain G@ in n-dimensional Euclidean 
space. The walk ends when the moving point leaves 
G for the first time. The nature of the random 
walk is characterized by the transition distribution 
function F (Q, P). Here P and Q are points of the 
space, and F (Q, P) is the probability that a point 
at P will, after one step, be in the rectangular sector 
of the space in which all coordinates are less than 
those of . 

By duration of a random walk we mean the total 
number of steps of the walk, including the one on 
which the point leaves G. The subject of this paper 
is the mean duration; that is, the expected value 
W (P) of the duration as a function of the starting 
point P. 

It will be shown that unde very general condi- 
W (P) is the unique solution of the integral 
equation (9). In III] assumptions are 
introduced, which, in a sense to be specified below, 
make the occurrence of all but very small single 
steps very unlikely. It will be shown that under 
those assumptions W (P) differs very little from the 
solution of a certain elliptic differential equation 
This differential equation involves only the first and 
moments of F (P,Q). The method used in 
these sections is closely related to that of Petrovsky 
[4].” 

For a particularly important of random 
walks this differential equation reduces to a special 
For such random walks 
it is possible to obtain numerical inequalities for the 
mean duration. These estimates confirm the plausi- 
ble conjecture that the mean duration does not grow 
very fast with the number of dimensions, a fact 
that is of interest in connection with the attempts 
sampling procedures for the 
numerical solution of partial differential equations 

If the moving point performs a continuous 
Brownian motion rather than a discrete random 
walk, the mean duration satisfies exactly a Poisson 
equation. This is an easy consequence of the fact 
that the fundamental solution of the differential 
equation for diffusion in the presence of absorbing 
boundaries can be interpreted as the probability 
density that a Brownian particle starting from a 


tions 


sec tion 


second 


class 


case of Poisson’s equation 


to use stochastic 


random 
mean 
mean step length, can be approximated by 


domain is studied for 
integral equation 


walk in a bounded 
duration satisfies at 
an elliptic differential equation 
duration in an important special case 


given point ?, at time 0 be at a point P, at 
without having met the boundary. See 
Smoluchovski, Drei Vorlesungen Uber Diffu 
Physikalische Zeitschrift 17, p. 568 to 
1916), P. Levi, Les processus stochast 
p. 272 to 273 (Paris, 1948 It has been point 
out to the author by M. Kae that this fact co 
used as a starting point for an alternative treatmer 
of our problem based on the theory of proba vilit 
measures in function space. 

In sections I to III, where the n-dimensional! 
differs from the two-dimensional one only by 
need for more cumbersome notation, we shall ¢ 
the arguments for two dimensions only. 


etc., 


or 


II. Preparations; Lemmas on Integral 
Equations 


The transition distribution function F (Q, P 
assumed to be bounded and Borel measurable f 
all P and Q and a distribution function with resp 
to Y. In order to guarantee that the moving poi 
leaves the domain @ with probability one we intro- 
duce a hypothesis which will be called assumptior 
(A 

Assumption (A). There exist two positive nul 
bers r and s, independent of ?, such that the prob 
bility is at least s that the abscissa of the moving 
point be increased by r or more in one step. 

We refer to [4, p. 431] for the simple proof 
the probability of leaving @ is indeed one, if assum 
tion (A) is satisfied. 

Denote by F,, (Q, P), m 
functions defined by 


1,2, ... , the iterat 
F, Q), P FQ, P 


F,.(Q, P 


F'n (QV, Ridk R, r) 


/G 


Here, and in the sequel, the Stieltjes differ 
are formed with respect to the first argument 
an integral without explicitly indicated dom: 
integration is to be extended over the whole sp: 
F,, (Q, P) is the probability of a point being 
domain r< re, Y< Ye after exactly m steps, st 
from P and having all of its intermediate pos 
in @ Here Ta, Yq) are the coordinates of ( 
consequence of assumption (A , and the bo 
ness of G, there exists a positive integer m 
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trve number ¢, less than one, such that 


dF,.(Q, P)<e, for all P in G. 


e G 


shall need the following lemma concerning 

ral equations formed with distribution func- 
For a proof, we refer again to [4]. 

there 

which 


is satisfied, 


u(P) 


mma 1: If assumption (A 
s exactly one bounded function 
es the integral equation problem 


(Su Wdk(QY, P)inG 
u(P 
ly P) outside G, 

where ¥(2) is a prescribed bounded and Borel meas- 
urable function. This solution can be found by 
means of iterations starting from any 
function that is equal to ¥(/) outside G, and bounded 
and Borel measurable in G. 

The solution u(P) of 


successive 


Corollary: (2) satisfies the 


neq ialities 
glob. WP)<u(P)< lou. b. Wd). 
Take as initial approxima- 


Proof of the corollary: 
, Which is equal to 


tion to uw (P) the funetion uw, (P 
gl.by(P) in G. Then 


u(P)= ful QdF(Q, P)>g.1.b.¥, 


and, by induction, u,>g.l.b.y. A the 
limit proves the first one of the inequalities (3). 
The proof of the second inequality is analogous. 

Lemma 2: Let $(P) be nonnegative and Borel 
measurable in G, and assume that assumption (A 
Then the integral equation 


passage to 


s satisfied. 


u(P u(@) dF(Q, P)+ oP 


e/ ( 


has a unique bounded solution in G. If 9, (P) isa 
nondecreasing sequence of nonnegative functions 


converging to @(P), and u,(P) is defined by 


un (Q)dF (Q, P)+¢,(P), 


U'(P)=lim u,(P 
u(P) is 


solution non- 


s the solution of The 
egative. 
is proved to be non- 


Proof: The sequence U, P 
induction, 


reasing by means of mathematical 
ising the relations 


u, (P)—uy(P o (YW) dF(QY, P)>0, 


( 


WO) dF (QP 
@Q P). 


In order to show also the boundedness of the se- 


quence u, (/’) we iterate (15) and find 


(dF (Q, R)dF(R,P)\+- 2 (P 


u Qdol | fk , R dk R. P)| 7 @ pP 


QV) dF, (Q,P)+ 02 (P), 


Q, 
G 


Repeating this procedure m times 
the integral equation 


Q) dF (Q, P)+¢,(P): 


assuming n> m 


Q) dF, (Q,P)+oet” (P 


is obtained, where 


o=+ (P)< +1) l.u.b. O(P). 


we hav ec 


if i the 1 .b. of u,(P) in G 


therefore, using (1 


denotes 


QD P 


Lo1:<¢eL.-a4 u. b. 
cL ‘ b @Q P). 


hence 


-b o(P). 


The sequence u P), being bounded and nondecreas- 
ing has therefore a limit, and this limit satisfies the 
integral equation (4), since a passage to the limit 
under the integral sign in (5) is legitimate. 

The uniqueness of the solution follows by a stand- 
ard argument from the uniqueness of the solution 
of the corresponding homogeneous problem. The 
uniqueness of the latter is assured by lemma | 
Corollary: Equation (4) has a unique solution, 
even if the assumption that ¢(P) be 
omitted. 

Proof: 
Lemma (2) can be applied to the integral equations 


nonnegative 1s 


Let >? bea negative lower bound of o P 


us (P u, (QdF (Q,P)—¢ 
G 


u.(P Us A, dk (QV, P - @ r 2 
ty 





u(P 


Le mma 8: 


u,(P)—u,(P) is then the solution of (4). 
If ¢ and @ are constants such that 


o~¢o P)< re) in G. 


and if W(P 


is the solution of 


WIP W(QdF(Q, P)+1, 


a ( 


then the solution of (4) satisfies the inequalities 


oW(P)<u(P)<oW'(P). 


Proof: The function 7(P oW P) is the solution of 


w(P u(Q)dF(Q,P)+¢ 
JG 

Subtracting (4) from this we see that u(P)—u(P 

is the solution of a similar integral equation with 

@—¢(P) as nonhomogeneous term. From lemma 

2 it follows that 7(P)—u(P) is nonnegative. The 

other inequality of (7) is similarly proved. 


Ill. The Integral Equatian for the Mean 
Duration 


Let p,(P) be the probability that the duration of 
a random walk starting at P be exactly n. Because 
of the remark following the assumption (A), we have 


» > P 


Pay 


if we as we shall—that assumption (A) is 
satisfied 


For a point outside G 


assume 


not in G) 


p,»(P bon (PP? 


where 6 » 38 Kronecker’s symbol. 


duration of a random walk is 


The expected 
by definition 


W(P)= do np,(P 


if this converges. Otherwise, the expected 
duration is infinite. We shall prove 
Theorem 1. If assumption (A) is satisfied, W(P 


is finite and is the unique of the integral 
equation 


series 


solution 


WP W(QdF(Q,P 


satisfies the rela- 


Proof: 
tion 


The probability p,(P) 


Pn dF (YU, P),n>0. 


Pn- i(P) 


W’,(P) >> vp(P), 


then we obtain from (10) the relation 


W.i(P)= | W.(QdF(@,P) +S ple 


Reference to lemma 2 completes the proof. 

If F (Q, P) is a discrete distribution with a 
number of discontinuities, and if the position of 
discontinuities is independent of P, then the int: 
equation (9) reduces to a difference equation 
simplest cases can be solved explicitly; see, 
3, p. 286.] 


IV. The Asymptotic Differential Equation 


In many applications the steps of the random walk 
are likely to be small. In order to describe mor 
precisely the sense in which this statement is to bh: 
understood we assume, following Petrovsky, that the 
transition function depends on a small paramete: 
uw in such a way that the following assumptions ar 
satisfied. 

Assumptions B). Denote by a, (P. ps), ba (P. zu 

i, k=1, 2) the first and second moments of F ((, P 
about the point P. Then 


a; Pip a; P)\ u+o(u 


ba(Pypy Bu(P)ut+o (pu). 
The functions o(4) may depend on P, but the 
to have the indicated order uniformly in G. TI 
a;(P) and By (P) are assumed to be in class ¢ 
G+ Cc" 

Assumption L . Let 
P0<r, with center at P. 


K.P 


Then 


denote the « 


PQOdF(Q,P 


K,(P 


here o (wu) depends on P and r. The relation 
hold uniformly for P in the closure of G, for 
fixed r>v. 

In the applications » may measure the sma 
of the mesh, if the walk takes place in a lattic 
smallness of the standard deviation, if the transi! 
distribution is normal, the sma'lness of th 
intervals between observations, if the walk is act 
a motion in time, ete. 

The definition of the By is readily seen to 
that Bir Boo—Bi2>0. We require, beyond this 

Assum ption ke). In Gand on its boundary 


B11B 22 - Bio 0. 
This assumption will enable us to make use « 


existence theorems in the theory of elliptic diffe: 
equations 
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ally, we have to require a certain degree of 
hness of the boundary C of G. In particular, 
sults remain valid for polygonal domains. 
sumption (S). C has a continuously turning 
nt. 

is restriction can be weakened. ‘The necessary 
nents can be found in [4, p. 438]. 

shall show that the function w=yuW remains 
as uO and converges to the solution of the 
ential problem 


_s 0 in G 


V4 


0 on 


btaVn ta taVartal, 


+ ay V, 


(12) 


ye 


Because of assumption (E) the differential problem 
is elliptic and possesses in G a unique solution. 

The basic idea of the proof is to show that V 
satisfies an integral equation little different from the 
one for W. In order to avoid extraneous difficulties 
near the boundary ( it is convenient to replace 
V(P) by a slightly different function V;(P) defined 
na larger domain G’’. This can be done by con- 
tructing a twice continuously differentiable mapping 


x’ = f(x, y;4), y’ = g(z, y;4), (x,y in G) 
which is continuous in 6, for 6>0, together with its 
first and second derivatives with respect to z and y, 
nd has the following properties: (a) It reduces to 
the identity for 6-0. (b) It is, for all 6, the iden- 
transformation a subdomain G’ of G that 
0) It maps G into a domain G” 
containing G in its interior for 6>0. 

For the explicit construction of such a mapping 
vith the help of assumption S) we refer to [4]. 


V; (P) in @ by 


in 
c) 


tends to G as 6 


If we define 


a V (zx z,y in G) 


,y), 
this new function is defined and twice con- 
iously differentiable in G’’ It tends to V(P), 
formly in G@ together with its first and second 
vatives. We extend the definition into the whole 
ie Eby setting 

V;(P) 


0 in E—@Q”. 


In order to estimate the integral 


J=| Vi(QdF(Q,P 


P in G@ we denote by K(P,é4) a circle about P 
| radius 7(6) depending on 6 but not on P, such 


K(P,é) is in G@’’ whenever P isin G. We define 
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J, | VidF, Jy= | 


V;dF, 
E—K(P,é) 


then J m - ds. 
Now if A is an upper bound of 
0<5<%&, 


V;|, for all 6 with 


. 


A 


r*(6), E 


J; 


< 


POdF (Q,P). 


K(P 


Hence, by assumption (L), 


J.=0 (p), for fixed 6. 13) 


For the calculation of J; we use Taylor’s formula, 
for P in G, Q in K(P, 4), setting 
V3 (Q) 


Vi(P)+M[Vs(P+Q1Vi(P+R[Vi], (14) 


where 


’ 


Q(V) 


~ 


Q(V (P\—Q(V(P), Pin K(P, 


R([ Vi] =PQ*p (P, Q, 8) 


with 


lim p ' Vv, 6 


> 


0 uniformly for P in G, Q in K(P, 4). 


16 
In these formulas z, y and &, n are the coordinates 
of P and Q, respectively. We now integrate (14) 
termwise with respect to Q over K (P, 6) and discuss 
the terms separately. We have, e. 
tion (B), 


g., by assump- 


dF, 


df 


by assumption (L), and hence 


Ma, + O\p), for fixed é. 17 


Similar arguments apply to 





y)d F, ete. dF ? in G, V?= V, in E—G 


The integral equations (27) and (28) have 
o(P, d)0(u LS solutions, by the corollary to lemma 2, a 
lemma 1, respectively. If we apply lemma 
> 7 P oF 
. 27), it follows from (9) and (26) that 
by virtue of assumption (B Here p,(P, 4) tends to 
zero uniformly in G. as 6-0, because of (1.6). 


and O(u)/u remains bounded as 6-0 Formulas u(l =z} WW - ra) « u ( 1+ 54)" 


1.8) and the formulas for dF, 
, Next we apply the corollary of lemma 1 to 
y)dF. etc.. analogous to 7). when | and use the inequality (23). This shows that 
KU 


applied to (1.4), vield 


J Vi(P pL V, p P, 5)0(u)+ o(p). 19 

The inequalities (29) imply 

Now let «>0 be assigned arbitrarily (but less than 

A), and choose a 4, (0<.5< 4), so small that the fol- 
lowing inequalities hold 


uW 


a L{Vs|+1,<¢/24A inG (20) | In fact, we have, for e< 


This is possible because of (11) and the properties ‘ 
, V as described above : SA < | + iA 


» 5)0 2 (21) — ar 
pilP, 6)0(u) <we/24 A “1) | and, therefore, from the first inequality in (2 


the remark after formula (18) is used. 


rails 


€ 
in G 
By the definition of V{" and inequality 


ee 
€ 1 Q"—G@ 9: Vi |S2A. 
[f this is combined with inequality (32), one of the 
two inequalities equivalent to (31) is at hand. Th 
second inequality is proved similarly. 

Finally, adding inequality (22) to (31) the desired 
formula 


This last inequality is the only point where the 
boundary condition in (11) is needed. 

After 6 has been fixed in this fashion we choose 
uu 0 so small that the o(u) in (19) satisfies the 
inequality 

| ) ul 


e, for |u! <po 
o(u)| < we/24A in G, for |u| < po. (24) /. 
Is proved. 


: , We summarize this result as a theorem. 
From (20), (21) and (24) it follows that (19) can be — 


written Theorem 2. If assumptions (A), (B), (L), (E 
—- , (S) are satisfied, then the mean duration W (P 
(QY) dF(Q,P ‘ P, u), 2: 
VQ) dF(Q, wg SH the random walk satisfies the asymptotic relatio 


where lim pW (P V(P), 


(P, u 7 €/SA for P in G, | pw! < po. 


i | 4 | W : | , uniformly in G, where V(P) is the olution of 
For the comparison of Vy with 4 we split | elliptic boundary value problem (11). 
into 


V. Bounds for the Mean Duration in a 
where Simple Case 


From the viewpoint of numerical applicatio 
preceding theory is of interest mostly as a bas 
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reasonable bounds for the expected duration 
exists several methods for this. But since 
vestigation was originally motivated by ques- 
onnected with the *‘Monte Carlo’’ technique 
ving differential equations experimentally by 
of random walks, we shall disregard methods 
amount themselves essentially to an approxi- 
solution of such a differential equation. 
:| procedure for obtaining bounds for W (P) is 
in [2]. Here we shall limit ourselves to the 
il case that Liu] is a constant multiple of 
we’s operator. Then it is possible to obtain 
h more precise information, even if the space is 
assumed to be n-dimensional. 
In other words, we require that 


(P)=k-i:n,a(P)=0, 1,K ,2,-+-,n(33 
The differential problem (11) now assumes the form 


AV +7 0in G,V=0 on C 


where A is Laplace’s operator in n dimensions. 
The constant & is related to the mean square of 


the step length, s?= 3) 6,,(P) by the formula 
l 


nuk tT Op) 
There is therefore a u, >0 such that for |u!<y, it is 
possible to use s* itself as a parameter u. Without 
ss of generality we can therefore set 


u=S? 


nd hence, 


AV +2n 0in G, V=0 on C 36 


For those arguments below, where we let n ap- 
roach infinity we need the further assumption that 
can be chosen independently of n. Random walks 
he type described above will be called svmmet ric 
The two most interesting special cases of the type 
msidered are 
The random walk takes place in an orthogonal 
f mesh length A with transition probability 1/2 n 
ch of the 2n possible directions. 
Here a 0, bie =h*bix/n 


mptions (B), (L), and (EK) are satisfied, if u is 
ed by 


{ 


Me h? 


35) and (34) are satisfied. 
lhe transition probability is normal with 


If we define 


formulas (34) and (35) remain again valid. 


Theorem 3 If the conditions (33) are satisfied 
the mean duration of a random walk in an_ n-di- 
mensional sphere of radius a satisfies the asymptotic 
formula 


a~ 


WP -(14 10) 


where r is the distance of P from the center, s? the 
mean square of the step length, and lim e(s)=0. 


The most interesting aspect of this formula is its 
asymptotic independence of the dimension 

Proof: If Gis a sphere in n dimensions, the 
solution of problem (36) is V=a*—r*. Theorem 3 
is therefore an immediate consequence of theorem 
2 and formula (34 

Corollary 1: Let v be the volume of the n-di- 
mensional sphere of radius a. Then the value of 
W(P) at the center 0 of the sphere is 


WO 


Corollary 2: 


lim W(0)/n= A(s), where A(s) #0 $1) 


n 


This follows by a simple application of Sterling's 
formula to the right side of the preceding formula. 

If G is not a sphere, the equality (40) can be re- 
placed by an inequality, if we make use of the 
following lemma. 

Lemma 4. Let P be a fixed point in n-dimen- 
sional space and denote by By} the set of all n-di- 
mensional domains of constant volume r¢ containing 
P in its interior and having a boundary which is 
sectionally in Class C@. Let u,(P) be the value at 
P of the solution of the differential problem 


Au- 0in G,u=0 on C. 


Then u,(P) assumes its maximum with respect to 
the Class |B}, if Bis the sphere A with center at P. 

Proof: Let r be the distance PQ and denot: by 
G,(Q) Green’s function in B, if one argument point 
is fixed at P. For the sphere A of radius a about 
P Green’s function is of the form 


Gx(Q)=e 2—m__ gin) 12 


where ¢,(n) is a constant depending on the dimension 
n. For any domain B we have 
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up(P G,(QdQ 
B 


Denote by vg(P), P>O the volume of the “level’’ 
surface Gs((/)=p For the sphere A we have, 
from (42 


v n(p) =e: + q?-nyns(—m (44 
We now make use of the inequality 
vp(p)<vx(p), for all B of {B (45) 


which is an easy consequence of the classical ex- 
tremal property of the capacity of a spherical shell. 
(ef. Szegé [6], in particular p. 588. Szegé’s proof 
for the three-dimensional case extends readily to 
higher dimensions.) From (43) we have 


Up P pdvz vp , | erate) |: Vp(p)d p. 


The first term in the right member is zero, for 
v_(0)=v and lim pvzg(p)=90, since pvx(p)—which 


vy (45) at least equals pv,(p)—tends to zero, as 

p-> @, in consequence of (44). A reference to (45) 
now completes the proof of the lemma for n>2. 
The modifications for n=2 are trivial. 

From this lemma and corollary 1 of theorem 3 we 
obtain immediately the following result. 

Theorem 4: The mean duration of a symmetric 
random walk in an n-dimensional domain of volume 
v satisfies the inequality 


) 


W(P)\<x pws (——") M »y2 "lI } e(s)]. (46) 


Corollary 2 to theorem 3 implies that the right 
member of inequality (46) tends to infinity with 
order n, as n>, provided the volume is kept 
constant. But an explicit calculation shows that 
the increase, for small integers n is slower. In fact, 
if we write the right member of inequality (46) in 
the form 


318 . 385 . 450 


VI, Random Walks in an n-dimensional 
Cube 


The n-dimensional cube is another domain for 
which problem (36) can be easily solved. It seems 
intuitively plausible that the value of W(P) at the 


center of such a cube will have the same o 
magnitude in n as the value for the sphere, ¢ 
(41). This is, however, not the case. We « 
fact, prove the following theorem. 

Theorem 5: For a symmetric random walk 
lem let Ws denote the mean duration of walk 
ing at the center of an n-dimensional sphere of y 
v, and denote by Ws similarly the mean du 
for a walk starting at the center of an n-dimer 
cube of the same volume v. Then 
lim Ws/W:=0 


n 


lim n*Ws/Ws = @, (e>0, arbitrary 
n— 2 


Proof: It suffices to consider the unit cube # 


: l : . , : 
fined by |z, - 5 (t=1,---,n). Theeigenfunct 


of Au in B are 


n 


-, 2,)= IL cos K, rz;, K, > 0, K, odd 
i=! 


and the corresponding eigenvalues are 


Ne, ---x,= 9°) K?. 
i=] 


These eigenfunctions are orthogonal. Their 


is 2~", since 


II cos? A,rr,dz,--- dz, 
I 


It follows that the multiple Fourier coeffi 
_x, of the representation 


n 


l » » KCK, ; DL ..-Lm), odd A 


aK, 


valid in B, are 


(Sf ennai 


7 


If we now substitute for uw in Au 


series 
% 
DK, K, 
(odd kK, 


comparison of coefficients shows that 


AK, , Kak, 


that is, the value of uw at the center z, 
cube is 


u(Q) T 2/ 


4 ” ‘ . , " 
) Dk, ...K, I(—1)4" 


i ches 


T 


(odd K,,...,Kzx). 
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n-fold sum can be simplified by means of the | where ¢ is an arbitrarily small positive number and 
tution c, and ¢, are constants less than 2/4. Therefore 


lim u(0)=0. 


(>) K - _ :, 
f=1 The first part of theorem 5 is now an imme- 
diate consequence of theorem 2 and formulas (36) 
wed by an interchange of summation and inte- | and (41). 
on, which is not difficult to justify (see, e. g. [1] In order to estimate the rate of convergence of 
331). Then we obtain u(0) from below, we observe that 


ie w/ J0 


x *( : ["ds, - II 1) se Kidz, u(Q)> T *( 4 ) [*srtede. 


kK, K,, odd), 


In the right member we use the inequality 


e~*dt— vr | é “dt, 


0 Jwiive 


2 


f(x)2- = + vr | 


¢ 


ior | - 
(7=0 29+1 ) 


sz. 47) 
which follows from (48), combining it with the 


In order to investigate the asymptotic dependence inequality 


0) on n we make use of the formula 
] 


). 


¢ ig" for 
T - 3 Te 13 (7 at Te 
V3 >> exp{—27?+(2T-+ rin)*/2 ; 
which is readily proved by an integration by parts. 
e, e.g. [5], p. 32), which is essentially the reciproc- Then we obtain 
formula for the @-function. If we make the sub- 8yz 
tution 7 Ly/a and integrate with respect to > 
between 0 and 7/4 we obtain from this formula 


T \ 
é 


w”/16z > *(1}~—927%)> 
)>4—22")> 4 


series » - ‘ . , ‘ , ; 
for 0<2<1; N>O, arbitrary. Inserting this in 
19), we have, finally, 


Yt Any r\dy - 
u(0)> 7 ‘| exp |[—2nz"|dz 
J0 


> const. n'/* 


This, in combination with formulas (36), (41), and 
theorem 2, proves the second part of theorem 5. 


yr in the integrals above, this be- 


VII. Random Walks in n-Dimensional 
Ellipsoids 


; vz 
e~P dt+ 2) -" edt, In this section we collect some results concerning 
a the mean duration of random walks in n-dimensional 

j (48) ellipsoids. 

trom (48) we conclude The “oblongness” of such an ellipsoid can be 
defined as the ratio of the longest to the shortest axis. 
f(2)< — : ke | e-dtan™ The mean duration of a random walk of the class 
; 2 0 4 considered in the last two sections will be less for a 
walk starting from the center of such an ellipsoid 
than it would be for a sphere of same volume. The 
ratio between these two quantities will give an idea 
. of the deterioration of the estimate (46) for oblong 
convex domains. For certain domains it may even 
4\"( - —_ be advantageous to estimate the mean duration of 
= ‘( | +| + | f(x) dx’ the random walk by its value in a circumscribed 
tT] Jo Je Ji” ) ellipsoid rather than by formula (46). We shall 
make use of a lemma concerning the arithmetic and 


<s-? e+(= a1) a—9+(- e) b, geometric means of n numbers which may be of 
us 


*T/4N *e(n+4 


%*@ 


ind from the definition of f(z), 


T(z)<e 


some interest in itself. 


51 469 





The mean duration of a symmetric 
n-dimensional ellipsoid with 


Theorem 6 
random walk in an 
semiaxes d,, 


We (P 211+ e(s)] +=) >s— (50) 


i a? rea 


Proof 
problem (36 


It is readily verified that the solution of 
in such an ellipsoid is 


y sot. 
a> rms | 


| a? 


This fact and theorem 2 prove formula (50). 

Theorem 7. If Wes(P) is the mean duration of a 
symmetric random walk for a sphere having the 
same volume as the ellipse, then 


l Vy 
€\s ad 


Proof: The volume of the ellipsoid is 


‘\ 


,) 
(x = I] a; 


If we divide formula (50) for P=0 by the formula 
in corollary 1 to therorem 3, and eliminate ¢ by 
means of formula (52), we obtain formula (51 

From now on we shall assume that the oblongness 
of the ellipsoid is prescribed. In view of the remarks 
at the beginning of this section it is of interest to 
the minimum of the right member in 
To this end we prove 


investigate 
51) under this restriction 
the following lemma 

Lemma 5: Let t; (@=1, , n) be n positive 
numbers such that the ratio of the largest to the 
smallest has a prescribed value p. Then the 
quotient of the geometric and the arithmetic mean 
of these numbers will be smallest when 


with 
log p 
Proof: Without loss of generality we can assume 


ee. ? *. a, £29. 


We have to minimize the quantity 


F (t,, (54) 


oy Ue 


1 t,/(=t,)" 


under the side conditions (53). As a function of t, 
alone, (l<s<n), F has its only stationary value 


when 


bor 
Mt lt 


This stationary value is a maximum, since F 

when ¢,=0 or tt=@. Hence, as a functior 
alone, F has no relative minimum and it as 
therefore its minimum for given values of 

either when t,=t,_, or when t,=t,,,. It follow 
at the minimum of F, the ¢ of two ¢ 
such that 


consist 


et 
the minimum will therefore be among the va 


vr ———, 
h r)p| 
We find that Yr) 


which occurs at 


has exactly one stationary 


nx(p)=n ( a 
p l log p 

This stationary value is a minimum, for y(r) is 

tinuously differentiable and positive in the int 

it is infinite at 7 

this interval 


ana 


D<_P<_Tip/\p 1), x 
Mor 


np/(p—1), and r, lies in 


lim x(p)=s> x(p) increases with p, and x 


Hence 


If ro is an integer, this is the admissible valu 
for which Y(r) assumes its minimum. Otherwis 
is one of the two neighboring integers. This con 
pletes the proof of the lemma. 

Corollary 1: For sufficiently small p>1 th. 


teger r of lemma 5 differs by less than unity 
so does rs For sufhie: 


n/2. As p increases, 

large p we have 7 n 1. Also. r is never less t 
" 2. 

For n=3 we have always 7 
One shows easily that, for p>1, 


Corollary 2: 
Proof: 
¥(1)/¥(2) 

Theorem 8: With the notations and assum; 
of theorem 7 denote by m(w) the minimum o 
ratio W,(0)/Ws(0) for all ellipsoids of given ol 


ness w. 


Then 


lim {m(w)-[1 + (s)]} = 


n-+ 


[x(-v co) wx Ved) + o?(1 — x(-¥@))] 7 
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of: Without loss of generality we set 
a,>a,> » 26a, 


ay 


a» 


We apply lemma 5 to formula (51) with 


minimum of the right side of (51) is then 


t-e(a)|nayp? "az?" —-?!" /ray?+(n—r)jaz? 


\s n> @, r/n—yx(p)=x(vo), and formula (55) is 
it hand. Theorem 8 shows that for sufficiently 
arge vi and w l 


i\@)« 


On the other hand, m(w) does not tend to zero as 
For large n, random walks starting from the 


center of any ellipsoid of prescribed oblongness are 
on the average longer than those in n-dimensional 
cube of same volume. If r is kept constant for all n 
in (56), it represents the ratio W,(0)/W,(0) for 
ellipsoids of given oblongness r of whose axes have 
maximal length while the remaining n—r have 
minimal length. For this type of ellipsoids 
W(0)/Ws(0) tends to lasn->o. The correspond- 
ing random walks are, therefore, on the average 
longer, for large n, than for ellipsoids of the minimal 
type. 
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Research | 


Effect of Chromium Plating on the Plastic Deformation of 


SAE 4130 Steel | 


Hugh L. Logan 


The effects of chromium plating on the plastic deformation of SAI 


1130 steel, hardness 


Rockwell C40, have been evaluated from the results of tensile, tensile impact, bending and 


tests 


crushing 


Chromium plating materially 


reduces the percentage of elongation, true 


stress at beginning of fracture, and the ratio of the original cross-sectional area to the area 


at beginning of fracture. 


ness 


Baking after plating generally increased the ductility of the plated steel. 


Values of these properties decreased with increased plating thick- 


The tensile 


and yield strengths in tensile tests and the elongation and reduction in area of tensile impact 
tests were generally reduced less than 10 percent by chromium plating and were changed 


only ‘slightly by baking after plating 


I. Introduction 


Electrodeposited chromium has found wide appli- 
cation in machine elements. Because of its hardness 
and the ease with which it can be applied, manufac- 
turers and maintenance shops have used it frequently 
for protecting softer metals from wear or for salvaging 
worn or undersized parts. These considerations led 
the Bureau of Aeronautics, Department of the Navy, 
to request and to give support to an investigation, at 
this Bureau, of the effects of chromium plating on 
the plastic deformation and fatigue properties of 
steels used in aircraft. The results of the fatigue 
investigation have been reported.!. The effects of 
chromium plating on the plastic deformation of 
SAE 4130 steel are reported in the present paper. 

Tensile, tensile impact, bending and crushing tests 
were used in evaluating the effects of the chromium 
plating on the plastic deformation of the steel studied. 
It is the usual commercial practice to bake chromium 
plated steel articles at a temperature of about 200° C 
for several hours following plating. Hence, in this 
investigation specimens were tested as plated and 
after they had been plated and given various baking 
treatments. 


II. Materials 


The steel used consisted of one lot of SAE 4130 rod 
% in. in diameter, and one lot of SAE 4130 tubing, 
having an outside diameter of 1% in. and a %-in. wall 
thickness. The compositions of this steel, as deter- 
mined by chemical and spectrochemical analyses, 
are given in table 1. 


TABLE l Composition of steel 


rubing 


Element 


Percent 
Carbon. . 0.31 
Manganese . . 1 52 
Phosphorus 009 
Sulfur O18 
Silicon . 27 
Chromium 1. 05 95 
Molybdenum 
Nickel 


21 
Not determined 


Effect of chromium plating on the endurance limit of steels 


tH. L. Logan 
RP211; also ASTM Proc. 50 


used in aircraft, J. Research NBS 43, 101 (1049 
(1950 


Ill. Specimens and Test Methods 


All specimens prepared from rod stock were rough 
machined, prior to heat treatment, to dimensions 
0.015 to 0.020 in. greater than those of the finished 
specimens. These specimens and lengths of tubing 
for crushing tests were quenched from approximate! 
870° C. into oil. All material was tempered for 1 hi 
at 495° C. and was then cooled in air. Following 
the heat treatment all specimens were machined | 
final dimensions. 

Specimens of all types were tested as machined 
after chromium plating to one or more thicknesses 
and after plating and subsequent baking at variou 
temperatures. Cleaning of specimens prior to plat 
ing has been described in an earlier paper. (S 
footnote 1.) Plating was carried out in an electro 
lyte containing 250 g/liter of CrO; and 2.5 g/liter o 
H.SO,. The bath was maintained at 55° C, and 
the current density was 350 amps/ft®. The chro 
mium plating thicknesses on the individual specimens 
were, in general, within +25 percent of the nomi 
thicknesses given below. 


Specimens were baked in an electrically heate 
laboratory oven, with forced air circulation, ther- 
mostatically controlled to +2 deg C at temperatures 
of 100° and 200° C. Temperatures of 300°, 400 
and 440° C, controlled to +3 deg C, were obtained 
in an electrically heated tempering furnace th 
forced air circulation. 


1, Tensile Specimens 


The tensile specimens were machined with d 
ters in the reduced section of 0.4375 +0.00 
and gage lengths of 2 in. Specimens were t 
(1) as machined, (2) after plating to nominal t! 
nesses of 0.0001, 0.001, 0.010, and 0.015 in., an 
after plating to above thicknesses and baku 
various temperatures. The yield strength, t 
strength, and percentage of elongation in 2 in 
determined on each of the specimens. True s 
true strain data were also obtained in the 
range extending from just above the yield st: 
to fracture on specimens as machined and 
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le im pac t tests made in accordance 


New York Naval Ship Yard. 


| Specimen for tens 


lrawing prepare 1 at the 


¢ thicknesses of 0.0001, 0.001, and 0.010 in., 
e conditions as plated and as baked. 
eld strengths were determined at 0.2 percent 
from the modulus line from autographically 
ded load-strain curves obtained with a Templin 
type high magnification extensometer. After the 
specimens had extended so that a yield load could 
be determined, the recording extensometer was re- 
moved. Diameters of these specimens on which 
stress-strain data are reported were measured 
manually, to failure of the specimens, using a dial 
gage caliper reading to 0.001 in. All tensile proper- 
ties were computed on the total diameters of the 
specimens, that is, diameter of the basis steel plus 
the thickness of chromium plating. True 
strain data were obtained from simultaneous read- 
ings of the applied loads and the minimum diameters 
of the specimens; values for true strains were com- 
puted as the log, A /A, where Ap and A are the 
original (including plating) and current minimum 
areas of the cross sections of the specimens. 


stress- 


2. Tensile Impact Specimens 


Tensile impact specimens were machined in ac- 


ordance with the drawing shown in figure 1. The 
specimens were tested at the New York Naval Ship 
Yard as machined, as plated to a nominal thickness of 
015 in. and after plating and various baking treat- 
ents. The tests were conducted at room tempera- 
ture, the striking velocity was 27.8 ft/sec. 


3. Bend Test Specimens 


The rod material used for bend tests was cut into 
)-in. lengths. These lengths were ground to a 
finished diameter prior to plating of 0.500 +0.001 in. 
Specimens were chromium plated nominal 

ness of 0.015 in. and were tested as plated and 
various baking treatments. Bend were 
in a hydraulic type testing machine, as shown 

gure 2, at a free (unloaded) cross-head speed 

05 to 0.10 in./min. The load was applied until 

specimen failed (usually with sudden fracture 
itil the limitations of the apparatus (angle of 
bout 145°) prevented further bending. 


to a 


tests 


FIGURE Chromium plated steel rod, in position for bend test. 


Supports were 10 cm apart 
4. Crushing Test Specimens 


Specimens for crushing tests were prepared from 
the steel tubing by machining to an outside diameter, 
of 1.480 +0.0005 in., an inside diameter of 1.273 + 
0.0005 in. and a wall thickness of 0.1033 +0.0005 in. 
Specimens were 0.902+0.001 in. long. Specimens 
were plated on both the inner and outer cylindrical 
surfaces to a nominal thickness of 0.010 in. They 
were tested to failure by compressing between 
parallel plates, the direction of the load being 
along a diameter of the tube. Tests were made in 
a hydraulic type testing machine at a free cross-head 
speed of 0.05 to 0.1 in./min. Failure of the specimen 
was indicated by a loud report and a sudden decrease 
in the load 


IV. Results and Discussion 


1. Effect of Chromium-Plating on the Tensile Prop- 
erties of SAE 4130 Steel Specimens 


The results of tensile tests on unplated specimens 
are given in table 2; data on specimens plated to 
four thicknesses, as plated and after various baking 
treatments are given in table 3. Many of these 
data, plotted in percentage of values obtained for 
the unplated steel, are shown graphically in figure 3 
True stress-true strain curves for the unplated steel 
and for specimens plated to nominal thicknesses 
of 0.0001 to 0.010 in., inclusive, are shown in figures 
5 to 7, inclusive. True stress-true strain curves for 
the unplated specimens heat treated alike but at 
different times (lotsa. b and ¢, table 2) differed slightly 
in slope and are shown with the curves for plated 
steel for which they are applicable in figures 5, 6, 
and 7 
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TABLE 2 Results of tensile tests on unplated SAE 4130 steel specimens 
. } } 


Data reported were obtained from 5 specimens of each group 


lensile strength 


Standard 
deviation 


Original 
men/area it begir 


True stress at begin 4 of spec 
* ning of fracture 
ning of fracture 


lhiin 
186, O00 
187, 300 
TSA, TOO . 
heat treated at different times under approximately the same conditions ‘ 
e reported in percentage of values for unplated steel specimens 
TABLE 3 Results of tensile lests on chromium plated SAE 4130 steel specimens 
1 4 0r 5 specimens except as indicated Total diameters (basis steel plus chromium plating) of specimens were used in comy D 
reported s 


i strength 
Nomina 
plate 
thick 
Percent Percent 
age of 
unplated 


age of Standard 
unplated deviation ' 


steel 


times under approxim 
e of the average values 


The effects of chromium plating on the tensile 
properties (of specimens tested as plated) were most 
noticeable in the reduction of values of (a) true 
(load/corresponding cross-sectional area) at 
beginning of fracture, (b) A)/A, 
tional area beginning of 
fracture), (c) percentage of elongation in 2 in., and 
(d) the flow stress at equal values of true strain 
The properties listed (a), (b), and (c¢) and the true 
stress-true strain relationships indicate the ability 
of the material to withstand plastic deformation and 
will be referred to hereafter as deformation proper- 
ties. The reduction in these deformation properties 
became greater as the plating thickness was increased 
from 0.0001 to O.OLO in 

The yield and tensile strengths of specimens de- 
with plating thickness as the thickness 
became greater than 0.001 in., figure 3, B and A, 
and were approximately 91 percent of that of the 
unplated steel for specimens plated to a thickness of 
0.015 in These data indicated that the tensile and 


st ress 
(original cross sec- 


cross-sectional area at 


creased 


vield strengths of the chromium plate were appreci- 
ably 


less than those of the basis steel 


Standard 


deviation 


Origin re 

lrue stress at be . he : 
men/are t 

ginning of fracture . = 


ing of ira 


Percent Percent Percent 
Standard 


deviation 


we of Standard ize of 
unplated deviation > unplated 


steel steel 


ize of 
unplated 
steel 


Baking of plated specimens for6toS8 hr at 200° ¢ 
prior to testing, produced substantial increases : 
the deformation properties over results obtaine 
plated but unbaked specimens, as is shown in | 
3,C,D,and E. For specimens plated to thick: 
of 0.0001 and 0.001 in., baking at 440° C pro 
greater increases in values of these deforma 
properties than were obtained by baking at 20 
For specimens plated to a thickness of 0.01 
some improvement was obtained in these prop 
by baking at 440° C for 1 hr, but this improy 
was not as great as that obtamed by bak 
200° C for 6 hr. Greater improvement in 
deformation properties might have been obtau 
baking specimens with the 0.010 in. thick plat 
a longer period at 440° C, preferably in an 
atmosphere or in a protective bath 

In general the vield and tensile strengths 
plated specimens were not appreciably chang 
baking at 200° or 440° C, figure 3, A and B 
apparent increase in the tensile strength of 
mens plated to thicknesses of 0.015 in. and ba 
200° C was probably not significant 
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Bae 


.005 
PLATE THICKNESS , INCH 





Relationship between <A, tensile strength, B, yield 
elongation, D, true stress at beginning 
train at beginning of fracture and plate 

plated SAE 4130 steel 


( plated and baked 


¢ nlage ot 


specimens. 


In making the tensile tests on the specimens plated 
hicknesses of 0.010 in. or more it was noted that a 
of the chromium spalled off the specimens 
after the maximum load had been reached 
ng this part of the tests, loads and cross-sec- 
nal areas were changing extremely rapidly, and 
her loads nor specimen diameters could be ac- 
tely determined; therefore the true stress-true 


n curves are represented by broken lines in 


se regions (fig. 5 
pecimens from which part of the chron:miimn had 
led are shown in figure 4. Measurements of the 
ced sections of these specimens at the conclusion 
he test and after the remaining chromium had 
stripped from the steel indicated that about 
» to .0035 in. of chromium remained on speci- 
s originally plated to a thickness of about 0.010 
There was no appreciable spalling of chromium, 








Specimens plated to a nominal thickness of 0.010 in 
and broken without subsequent baking 


FIGURE 4 


Note that a part of the chromium plating has spalled off from the reduced 


section 


during tests, from specimens plated to thicknesses 
of 0.001 and 0.0001 in. 

True stress-true strain relationships plotted for 
specimens chromium plated to a nominal thickness 
of 0.010 in., figure 5, were well below that obtained 
for specimens of the unplated steel. There was a 
considerable amount of scatter in data obtained on 
all specimens plated to this thickness. The scatter 
in data for specimens plated and baked at 200° C 
was less than in that for specimens tested either as 
plated or after they had been plated and subsequently 
baked for 1 hr at 440° C 

True stress-true strain relationships for specimens 
plated to thicknesses of 0.001 and 0.0001 in. (figs 
6 and 7) were closer to those for the unplated steel 
than was the case for specimens plated to a thickness 
of 0.010 in. ‘True stress-true strain relationships for 
specimens plated to nominal thicknesses of 0.001 
and 0.0001 in. and subsequently baked at 200° C 
for 6 to 8 hr, approached those for the unplated 
steel much more closely than did the relationships 
for the unbaked specimens. Relationships for speci- 
mens baked for 1 hr at 440° C approximately coin- 
cided with those for the unplated steel 

The electrodeposition of chromium is_ electro- 
chemically a very inefficient operation. Only 15 to 
20 percent of the current is effective in depositing 
chromium, and large quantities of hydrogen are re- 
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3 4 5 
TRUE STRAIN, LOG, A,/A 
true strain irves for specimens plated 


réss-l 


of 0.010 in. 


, lens tested plated and baked 6 hr 
, Specimens plated IV, unplated specimens testé 
without baking . rinning of cture for individual plated specimens: @ 
of fracture for unplate« ‘0 data were available in broker 
} > on 
N e AeA 


reCaUse 


to a nominal thicknes 
at 200° ¢ 





be ginning 
sections of curves 
stress-true strain curves for specimens 


of st ‘ 
; 7. True 

to a nominal thickness of 0.0001 in 

d baked ¢ 


1s plated; II, specimens plated and 
und baked | hr at 440° C; IV, unpl 
ng of fracture for individua 


Degini 


mens l ited Spx 
without baking; plated spe 
zy of fracture for unplated spe 


Spex 


ig 





. ain cuves jor 
nominal thickness of 0.001 


ed as plated; IT, 
ited and baked 











iGURE 8. End of fractured specimen showing central spot 
differing in appearance from remaining area; X &. 


It has been reported ? that 
much as 128 


considered 


the cathode. 
of chromium 


leased at 

volume contains as 
volumes of hydrogen. It is generally 
that hydrogen produces embrittlement in steel. In 
the present investigation atomic hydrogen was de- 
posited on unplated tensile specimens by making 


them the cathodes in a 10-percent H,SO, solution. 


Specimens on which tests were begun approximately 
min after they had been removed from the bath 
broke before yield loads (0.2-% offset from the modu- 


The elapsed time from the 
emoval of the specimens from the bath to their 
racture was 6 to 8 min. ‘Testing of one specimen 
vas delayed for 5 to 6 min after it was removed from 
solution and a true stress-true strain curve very 
similar to that of specimens tested as machined was 
tained. ‘These results indicate that the effects of 
epositing atomic hydrogen on a steel specimen are 
transient. It 1s difficult to see how hydrogen, 
osited on the steel during the chromium plating 
would be retained unless there was diffusion 
m the chromium into the steel and a state of 
librium was set up. If such a condition existed, 
ng of specimens at 200° C, which is reported * to 
xpel approximately half of the hydrogen from the 
hromium, would be expected to decrease the hydro- 
in the steel and hence increase the plastic de- 
lation prior to fracture. Baking at 440° C if 
tinued for a sufficiently long period should remove 
than 90 percent of the hydrogen from the 


Is line were reached. 


rocess 


chromium and might be expected to further increase 
the plastic deformation that could occur in the basis 
steel before fracture occurred. 

Tensile specimens plated to a nominal thickness 
of 0.010 in. and broken without subsequent baking 
showed areas near the center of the fractures that 
differed markedly in appearance from that of the 
remainder of the fractured surface. <A typical area 
of this type is shown in figure 8. No markings of 
this type were found in specimens that had been 
baked before testing nor in specimens with thinner 
chromium plating nor in specimens loaded cathodi- 
cally with hydrogen. No satisfactory explanation 
for the occurrence of these markings has as yet 
been proposed. 


2. Tensile Impact Tests 


Tensile impact tests were made at room tempera- 
ture on specimens as machined, after plating to a 
thickness of approximately 0.015 in. and after plating 
to this thickness and baking for periods of 24 hr at 
100° C, 6 hr at 200° C, and 1 and 2 hr at 300° C 
The results are given in table 4. The data indicate 
that chromium plating slightly reduced the energy 
absorbed by the specimens and also the ductility as 
measured both by percentage of elongation and 
reduction in area. Baking of plated specimens at 
temperatures up to 300° C produced no significant 
changes in these properties of the plated specimens. 
No data were obtained on specimens heated at 
temperatures above 300° C 

Seabrook, Grant, and Carney * have shown that 
the effects of hydrogen embrittlement decreased 
with increased strain rates. The results obtained 
in the tensile impact tests and tensile tests are in 
ayreement with their conclusions if the reductions 
in plastic deformation of the steels under test are 
due, to hydrogen embrittlement. 


3. Bend Tests 


Bend tests were made on specimens as machined, 
after plating to a thickness of about 0.015 in. and 
after plating and baking for periods of 24 hr at 100 
C, 16 hr at 200° C, 7% hr at 300° C, and 1 hr at 
440° C. The data are given in table 5. The ap- 
plied loads reached maximum values while the angles 
of bend (180° minus the angle included between 
the two legs of the specimen) were still small and 
thereafter decreased until the specimens failed or 
until the limitations of the apparatus prevented 
further bending. It was noted that unplated speci- 
mens were bent through an angle of about 145 
without failure. Plated specimens failed after they 
had been bent through angles of approximately 40 
Specimens plated and then baked before testing 
could be bent through angles of 70° to 80° before 
failure occurred. Specimens plated and subsequently 
baked at 200° C or some higher temperature had 





Results of tensile 


TABLE 4. 


tal energy 


impact tests on SAE 4130 steel specimens chromium plated to a nominal 


Red 


thickness of 0.015 1 





uction In area 
Numt I 
. - Standard - Standard : Standar 
remperature Average deviation Average deviation Average leviatic 
h ( t-lh ft-lt Percent Percent Percent Percent 
N Not plated 44 27 15.8 04 55.2 0.8 
Pla bake Plated baked 431 19 14.9 ‘ 53. 6 . 
24 100 43 14 15.0 s 54.2 1.6 
f 200) 42 23 14.4 1.2 51.2 2.1 
t wi Os “MO 14.1 1 7 52.4 2.1 
Lb $16 44 14 1.0 M1 2 
approximately the same moduli of rupture [(9.975) | Taste 6. Results of crushing tests on chromium plati 
: ~ aa 4130 } 80 r by 0.902 
maximum load/ (diameter of specimen)*] computed on 130 steel tubing, 1.480 in. tn diameter by baa 
‘ 0.1033 0.0005 in. wall thickness 
the total diameters of the specimens, as the unplated 
steel. 
Baking treatment Breaking load Defe 
TABI E 5 Results of bend tests on chromium plated SAE 4130 Stand 
: ; sd- 
steel rod 0.500 in. in diameter plated to a nominal thickness ard 
= Time lemperatu \ ge A verage 
of O.U015 in devia 
tion 
Baki Stents Ar t Maximum Modulus of Hours < Li Li n 
. S tive load rupture Not plated Not plated 2 900 oo 0. 202 
Plated,not Plated, not 650 330 O1LS80 
baked baked 
Star Stand Stand 8 100 4, O45 aw 0610 
I per \ rd Aver ard Aver ard 6 200) 4, 285 45 1110 
‘ g levia ge levia 1 TD 4, 28 16 111 
tion tior l 1) 4, 155 12 1m 
OO * Deformation is original dia { plated specimen minu 
Ar ( Deare De Li Li lhiir lhiinx of specimen after failure 
ly pl 1 la 14 < 4 { 42 2 
plate Plates ‘2 6.020) 80 | 400 V. Summary 
ba bake ' 
' 10 ) 6, 130 1 10) ! 
a 4 . 4B 4 - 1. The effects of chromium plating on cert 
140 82 S 6330) 60 | 42 mechanical properties of SAE 4130) steel wer 
determined from tensile, tensile impact, bend 
® Specimet inott iK 


crushing tests. Data were obtained on specimens as 
plated and after baking at various temperatures 

2. The most noticeable effect of chromium pla 
was to appreciably reduce the plastic deformatio 
prior to fracture in all types of test except the tens 
impact. The reduction of elongations in tensile tests 
was particularly noticeable. 

3. Tensile and yield strengths of plated sper 
decreased as the plate thickness was increased 
0.001 0.015 in. However, these prop 
determined on specimens plated to a thickness 
0.015 in. were 91 percent or more of the correspo! 


Spalling of chromium from the steel usually began 
before the maximum load was reached; all or nearly 
all of the chromium had spalled from some of the 
specimens, in the region of maximum deformation, 
before failure occurred 


4. Crushing Tests on Tubing 


’ to 
Crushing tests were made on tubes 1.480 in. in 


outside diametet by 0.902 in. long with a 0.1033-in. 
wall thickness. Specimens were tested as machined, a f the | steel! 
after plating on the inside and outside surfaces to a | PPOPETUes O ee 
thickness of about 0.010 in. and after plating and : Phe plastic deformation gen bee ire a 
baking for periods of 8 hrs at 100° C, 6 hrs at 200° C, tensile tests, became less as the plating thicknes 
1 hr at 300° C, and 1 hr at 400° C. The data are increased. 

given in table 6. It was noted that plating increased 
the load hecessary to crush the specimen and reduce 
the deformation at failure to about 9 percent of that 
of the unplated steel. Baking of plated specimens 
increased both the breaking load and the deformation 
at fracture. After baking the specimens for 6 hrs 
at 200° C, the average load required to break the 
specimens Was approximately 1.4 times that required 





to crush the unplated specimens; however, the A 8 C 
deformation was only about 55 percent of that of the - 7” 
—_ . FiGuRE Y Tube specimen hat had failed 1 he 
unplated steel ['ypical specimens that had failed 
: . A, Unplated specimen; B ted 
in the test are shown in figure 9 €. aaaitinis ticked cit aed ok a 
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Baking of chromium plated steel specimens at 
eratures of 100° to 440° C for various periods 
ased the plastic deformation that could be 
iced in the plated steel above that obtained on 
mens tested as plated. The tensile and yield 
gths of the plated specimens were not appreci- 
changed by the baking treatments. 

It is suggested that hydrogen deposited on the 
with the chromium may be a factor in reducing 
ymount of plastic deformation that the steel can 
stand prior to fracture. 


Specimens used in this investigation were plated by 
personnel of the Bureau’s Electrodeposition Section 
under the direction of W. Blum and Vernon Lamb. 
Tension impact tests were made by personnel of the 
New York Naval Ship Yard, Brooklyn, N. Y. The 
following members of the staff of the Bureau’s Metal- 
lurgy Division have given valuable assistance in the 
experimental work reported in this paper, Harold 
Hessing, Houston Babb, and Thomas P. Royston, Jr. 


> 


WASHINGTON, October 30, 1950. 





Researcn Par 


Stochastic Processes and Dispersion of Configuration: 
Linked Events 


Chan-Mou Tchen 


By the use of a function for the transition probability, as introduced by Kolmogoroff in 
diffusion problems with continuous motion (the so-called continuous stochastic processes), 
the dispersion of the configurations of linked events is studied. A long-chain molecule of 
irregular configurations, representing configurations of events linked in a time sequence, is 
chosen as a model, Attention is given in particular to the inversion of the dispersion law of 
stochastic processes, and to the interactions between nearest neighbors. A very elementary 
application is given to the perturbance on the distribution function for the configurations of 
a chain molecule, as a result of interactions by other neighbor chains 


we will discuss certain general properties of th 
persion function and other derived functions 
will show how such considerations of stochasti: 
esses are relevant to the description of the confi 
tions of long chain molecules. 


I. Introduction 


We can distinguish problems concerning the motion 
of particles as occurring in physics and mechanics 
into two types: (1) problems relating to completely 
determinate systems; (2) problems relating to irregu- 
lar systems 


The particles (atoms, monomer units, etc.) of 
long chain molecule are linked in a random arrange 

In the problems of the first type, the motion of the | ment. This may be considered as a spatial mo 
particles is controlled by certain known equations, so | of the configurations of events linked in a tim 
that its whole course can be derived from prescribed | sequence. Thus the variable s in equation 1, which d 
initial conditions. In the problems of the second | termines the time sequence of the occurrence of 
type, on the contrary, the motion is not exactly | position r is, in the model, equal to the number o! 
known, even if the initial conditions were given. | links separating the given particle from the orig 
Such cases are found when particles are subjected, for | of the chain. For the sake of abbreviation s w 
instance, to irregular influences of mutual impacts | be called “link-number”’. We consider a great nun 
and thermal agitations. In such cases it is onlv | ber of chains of identical overal structure, but whic! 
possible to calculate a distribution of the chances for | at a given instant give configurations independet 
the occurrence of particular positions of the particles | from each other. A chain will be called OA,, A,, t! 
at a later instant, that certain initial | particles A,,, A, being labeled by their respectiv: 
conditions have been given. Such systems can be | link-numbers s’ and s. A _ system of coordinates 
with fixed directions is attached to the origin ( of t! 
chain. The position of the particle A,, in the syste 
of coordinates is given by the vector OA,,=r’ 
components r’;, 7’, 7r’s. It is assumed that 
given instant there are many chains that hav 
A,,-particle in the positions between r’, and r’ 
Let 


assuming 


called “stochastic processes” or “irregular processes.” 


It is assumed that at a given instant s’ there are 
many particles with between rr’; and 
r’,+dr’, (in a volume element dr’; dr’, dr’;), where r’ 
is a vector with components 7’), r’s, r’3, and 7~=1,2,3. 
Let the density distribution of these particles be 


positions 


ir dr’; dr’. dr’ 
Let 
be the probability for the particle A,’ to be sit 


in the region between r’, and r’,-+-dr’ 


be the probability for a particle that at the instant 
s’ started from the region between 7’, and r’,+dr’, 
to arrive in the region between r,; and r,;+dr,; at the 
instant s. The transition probability function Pp 
will be called the dispersion function, since it con- 
trols the dispersion process of particles as due to 
irregular motion. 

The above general description of the stochastic 
processes has been studied by Kolmogoroff (1931, 
1933) in the diffusion of particles by continuous 
motion. The present author has extended the 
theory in such a way that it can be applied to the 
dispersion of the elements of volume of liquid in a 
turbulent motion (see Tchen, 1944). In this paper, 


In an analogous way we can define 
dry a? 


as the probability for the particle A, to be sit 
in the region between r, and r;+dr,. The trai 
from the state (s’,r’) to the state (s, 7) can be di 
by a transition probability 


> 8,7 dr; dro di 


It is the probability for a particle A, to be sit 
in the region between r, and r,+ dr,, being give 
a preceding particle Ay was situated in the 





en r’, and r’;+dr’, The function p has a 
lete formal analogy with (2), and in the present 
‘ation will be called the dispersion function of 
nked particles. 
s to be noticed that the introduction of a dis- 
nuous link-number s in the above picture is a 
er of convenience. In fact the picture can very 
be generalized to a continuous chain, where we 
to with “elements of chain” instead of 
ed particles” Then the position s of the ele- 
in the chain can be determined by its contour 
h from the origin, which length may be expressed 
nsionlessly, for instance, in terms of a “‘corre- 
length’’. About the correlation length, see 
n 1947 
me properties of the dispersion function and its 
rential equations will be discussed briefly for the 
dimensional case in section I], by following the 
dimensional procedure of section 5 of Tchen 
14) and section 1.3 of Tchen (1947) in the study 
fusion problems 


be 


do 


remarked that the above dispersion 
tion refers to increasing s. However, the inverse 
oblem might formulated: we could 
to find from which preceding distribution an observed 
stribution may have originated. If s refers to the 
time variable in the dispersion processes of discrete 
stems, this will mean to determine the past from 
future distribution. Such an inversion of a law 
nature presents special difficulties and will be 
in section III. (The notion of inversion 
fa law of nature has been introduced by Schrédinger 
As an application, we may consider cases 
here the distribution of the positions at the final 
nstant can be perturbed by interactions between the 
stributions of positions at the earlier instants. In 
rder to calculate such interactions, it will be neces- 
iry to know the distributions of positions at the 
tial instants, being given the distribution of posi- 
ons at the final instant. Such a transition from the 
nal state to the initial state is an inversion of the 
spersion processes 
In order to study further details in the mechanism 
nteractions, in particular when we consider two 
rticles as nearest neighbors, or when we consider 
nteraction between the far particle of a molecule 
ticle situated in the farthest perimeter) with 
rounding molecules, it will be necessary to formu- 
te the nearest neighbor and farthest neighbor 
robability functions. Those functions will be 
tudied in section IV. While the derivations of the 
egoing functions and the analysis of their proper- 
form the subject matter of this paper, elucidating 
mathematical basis of the dispersion of chain 
gurations, a very elementary application is 
en finally in section V, in order to illustrate as 
xample, the perturbance on the distribution as 
ult of interactions. This at the same time will 
nish an opportunity of pointing out what reason- 
ngs are underlying the mechanism of interactions, 
where the difficulties arise. 


to 


also be ask 


scussed 


ane 


II. Properties of the Dispersion Function and 
Its Differential Equations 


The distribution function W and the dispersion 
function p must satisfy the relations 


®t 


The notation 


di represents | drs | dr; 
for the sake of abbreviation. Formulas in (3) ex- 
press the condition of normalization, and formula 

4) expresses the condition that all transitions start- 
ing out of the state (s’,r’) must find their places in 
the totality of the elements dr, - - -, forming the avail- 
able domain of the variable r;. p is the transition 
probability from a known initial state (s’,r’) to a 
final state (s,r). So, if the initial distribution W(s’r’) 
is known, the final distribution W(s,r) can be deduced 
from the relation 


Wis,r 


In order conform to the laws of diffusion 
phenomena, the dispersion function moreover must 
satisfy the condition (see Smoluchowski, 1913 and 
Chapman, 1928 


Lo 


p(s’, r’ 


. means that the 
diffusion process is integrable, 7. ¢., the process from 

sr’) to (8, controlled by the function 
p(s’ can be obtained by calculating first the 
dispersion in the interval s”—s’, and then the dis- 
persion in the adjoining interval s —s” being an 
arbitrary intermediary variable within the interval 

; 


Ss Ss 


where s’<s”<s. This equation 


as 


; 
8,7 


Under the conditions expressed by (4) and (6), 
the dispersion function will satisfy two partial dif- 
ferential equations given originally by Kolmogoroff 


(1931 and 1933). They have been derived by 

Tchen (1944) independently of Kolmogoroff’s meth- 

ods. In the following lines, we shall first discuss 

some developments into series of the dispersion func- 

tion before returning to those partial differential 
| equations. 
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It is in the nature of the phenomena of motion, | sembles to those of the Dirac operator, the 6-func: jo, 


whether regular or irregular that we must expect which is a kind of improper function. Certain 
ae ; ations involving an improper function may in pa 
lim p(s’,r’; s,r)=0 for r#r’ the rigour. For instance, the developmen: oj 
— p(s’, r’; 8,r) into a Taylor series will become ly 
i. e., when there is a finite difference between r and r’. | Convergent or even cease to be convergent at th 
Nevertheless relation (4) must remain valid. Con- | limit s=s’. As an example let us examine the 
sequently for r=r’, the function p(s’, r’; s, r), | development in the derivation of the second p 


which will be finite for ss’, must take necessarily | differential equation of the p-function by Kol- 
an integrable infinite value, when s tends indefinitely | mogoroff. Kolmogoroff (1931 and 1933) intro: 
to s’, so that (4) is not violated. Such properties | an arbitrary function R, that satisfies  certaiy 
of the p-function at the limit when s tends to s’, re- | boundary condition and showed that 


“4 Op(s’ rs 8,7 
di Ror) 
a Os 
lim dy Ro dr” p(s’, r’;8,r”) p(s,r"; s+r,7 lim d) Rr) pls’ .7 
, np Oo. ,; ” ” 
lim di di ~ Ro -—=> ins ,fr 58,7 ) W8,? s+rir)(r—) 
, “) e . ( or ' 
oO” , ” ” / yr” / r” } 
+ — — Pp S } S.7 ) PAs rss T.7 
or” ,or”; 2 1} 
The quantity in is a development into a ie “4 i b Pita ae nl 
Taylor series, which may become badly convergent é‘ ore e a9 re oe 


for. values of + approaching zero. In view of this - 
difficulty, special care about the series development a, *4 
will be taken in our derivations. dl, ... Ul, P(s,r; 17,1 

In Gauss’ function the parameters s’, r’, s, 7 : 
figure exclusively in the form of the differences 


s—s’, r—r’. In functions of a general type, how- Following the assumptions of Kolmogoroff, moments 


ever, the parameters may occur separately. Let of higher orders are negligible and the ratios 
=. = 
s—e’=r: r—p’=l l/r, ll,/2r tend to finite constant values at t 
limit 7-0. However, we can make, as was don 
then we can write the dispersion function p(s’, r’; s,r) | Tehen (1947), a broader assumption: we assum 
as a function P of (s’, r’; r, J) or of (s—r, r—l; r, L) that they tend to finite constant values for sma 
values of r but still exceeding the correlation measur 
p(s’, r’; 8,7,)=P(s’,r’; +r, D=P(s—r,r—l; 7,D. which was studied in detail there. 
; In order to derive the partial differential equations 
The second mode of writing is useful when we want | for p, we start from the integral equation (6). W 


” 


to express that the variations of a dispersion function suppose that r=—s”—s’ be a small quantity 
are much slower with respect to s’, r’ than with | write /=r"—r’ .' Ifs—s’ is not particularly si 
respect to r, J, in particular when 7 is small: neither will s—s” be small. So we can develop 


oP /dr!<oP/dl, for r small. eee P er 

p(s” ,r”; 8, r)= pls”, r’ +1; 

Hence in the development into a Taylor series that 

will be used in the derivations of the partial differ- into a Taylor series in powers of / by keeping 

ential equations, we will use the dispersion function | dispersion function in the form of p. Substituting 

in the form of p when + is not small, and in the form | this series into (6), applying (4) and (7), and 

of P when r is small. A development of P into a | tend tos’, we obtain 

Taylor series with respect to r means a development 


of p(s’, r’; 8, r) intoa series proceeding simultaneously ol = 

with respect to r’ and to r, with equal increments of op(s’,r’; sr) ke, r’) OP lik Op 

both variables. Os’ c ' ‘ee fea & 
With the use of the dispersion function, the mo- 

ments or the mean values of the displacement and of I aE ee Oe a 

its higher powers can be definied as follows: function considered is p(s’,r’; #”,r’’), such as in (6), then r=s 
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H and in the following, a summation will be 
stood whenever an index occurs twice. 


the second place let us suppose that r—s—s 


s , small quantity. When we put /=r—r’, we 
write (6) in the form 
S. 7) 
| dl . P(s’,? s” r—l P(s”’,r—l: r,.D. (6a) 
\ "—s’ is not small, we may _ develop 


‘; 8”, r—l) P(s”, r—l; 7) into a Taylor series. 
Substituting this series into (6a), we obtain after 
gration with respect to d/ and applying (4) 


a i, 
Or, 7; 2 


obtain 


i Oo 
iS. 7) } — 
TP I Or; Or; 


and (9) are the partial differential 
quations for the changes of p with respect to s’ and 
respectively. 
\ partial differential equation W can also be de- 
rived by an analogous procedure as for p. We start 
, Which can be written as follows 


rm 
(/.(s”,r) p]+ 


»| 


] 


to) 


By decreasing indefinitely 7, we 
ee" 

OWS 15 8,7 ofl 
Os Or T 


Kquations (Ss 


from (5 
Wis.7 dh... We—s.r 


, The integrand can be de- 
veloped into series, and by using the assumptions of 
holmogoroff, we obtain at the limit r—0, 


ith 7 S s’ / r 


Wd (wi a 
O 0 , rae of 

oe or (HW : +55, | Ws ): (10) 
rhe equation in this form is called the Fokker- 


Planck equation See Burgers 1941, Fokker 1914, 
Planck 1917 and Tchen 1944.) 


Ill. Inversion of a dispersive phenomenon 


Being given relation (5), we seek a _ function 
‘(s’,r’; s,r) such that for every pair of functions 
V(s’,r’), W(s,r) connected by relation (5), we have 
he following relation: 


Wis’. r’) dr Wis ir p'(s’,r’;8,r) (11) 


If such a function is found, then p* is called the inverse 
dispersion function. Now, in general, it is probable 
that such an inverse function will not always be 
positive. As an example let us take a discrete sys- 
tem in which the integral is replaced by a sum: — 
Wi= DOW; Dis t,j=1,2,---N, 


(12) 


for all values of j, 


and 
Pis= 0. 


We can solve this system of equations for the vari- 
ables W, when W; are given. Let D be the deter- 


minant of the pi, 


P Pi2 “ee. P N 
P2 

D Ps (13) 
PN eee PNN 


Let us denote by D,, its cofactor. Then we have 


, Dus orm 
+. ( 
W, Dp) Ww, 1 4) 


so that for every arbitrary distribution of the W’s 
the relation is inverse. 
Now there exists the known relation 


1, if j=k 
0,if 74k 
As p> 0 for alli and k, it is not possible that always 


Dy; 
>0 
— 


When negative values occur among the D,,/D, these 
quantities cannot be considered as determining a 
probability. Thus we cannot use them to answer 
the question: What is the probability that a par- 
ticle that is situated at r at the phase s (the variable 
s may be the time), was situated at r’ at a previous 
phase s’? 

As an example we consider a case of travels from 
the cities C, and C, to the cities C; and ©,, For 
each city of start, we can define the probability of 
Cc; or ¢ 


a travel to » independently of the number 


of travelers leaving the cities of start, e. g., 
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from C 0.6 to OC), 0.4 to C; 


from (C; 0.4 to C,, 0.6 to C». 


Thus relation (12) can be written into the equations 


W,—0.6 W,+0.4 W,) 
W.—0.4 W,+0.6 W,) 


(15 


where W,, W, are numbers of travelers in cities of 
start ©, and C,, and W, W, are numbers of travel- 
ers in cities of destination C; and C,. 

can solve the system of (15) for the 


We have 


Inversely we em 
variables W,, W, in terms of W,, W2. 


W,—3i 
W,=—2 


ow, ) 
Ww.) 


Wy 


} 


It is seen that p* (coefficients of W, and W;) presents 
also negative values and hene e cannot be considered 
as determining a probability. 

The dispersion function p as governed by the 
differential equation (9) has the task of bringing the 
distribution at a certain initial phase s’ into a broader 
distribution at a later phase s. So, on the other 
hand, in the inverse problem, we must expect a 
narrowing of the distribution. However, such an 
inverse dispersion in the narrowing transition of 
distributions cannot exist without having certain 
negative values either in the inverse dispersion 
function or in the distribution function. From the 
above considerations we can conclude that in a 
stochastic world, we can find a law of transition 
(dispersion function) that tells us about a distribu- 
tion of the future from a distribution of the past. 
However, for someone living in the inverse world 
where all events present in the wrong way, it will not 
be possible to find an inverse law of transition (in- 
verse dispersion function) that tells about a distri- 
bution of the past from a distribution of the future, 
if the law of transition and of distribution are to be 
positive. 

In order to answer in the inverse problem the 
question: what is the preceding distribution from 
which an observed distribution has originated, we 
may write: 


for the case of travels: 


W, W, ' 


Wi, 
W,=0.6 +0.4 —' W,, 
l , W, 2 


w,=0.4 2 Wi+0.6 4 wy, 
V; V, 


for a discontinuous system in general: 


fu nection . 


WwW here 


and for a continuous system: 


%+o 


Ww 


W(s’,r’) 
W(s,r)- 


; 


7? ; ; 
P\s ,? 58,7) PS ,r; 8,7 


We may consider p’ as the probability for a par 

to have started from dr’, at the phase s’, being 

its position in dr, at the phase s._ In a certain sens 
the function p’ plays the role of a kind of invers 
function of the dispersion function Pp. Howe ver 
p’ is dependent upon the distribution really present 
and has no independent meaning such as was th 
case with p. p’ will be called retrograde disp: 

By eliminating p from (5) and (17 
seen that p’ satisfies the condition of normalizstior 


. 


dr’ 


Let us now specify the above considerations of 
retrograde dispersion for a chain molecule in th 
description of the intramolecular interaction 
sider a chain A,A, A, with origin at A,. 
given the position r’ of Ay, the probability for tl 
particle A, to be situated between r; and r 
by definition 


p’ 


Since the chain is symmetrical, we can tak: 
form 


depending upon the difference 
Hence (17) becomes 


» 7) 


This relation means that the  probabilit 
p’(s’,r’; s,r) dr’, dr’,dr’; for the particle A, to b 
situated in dr’,dr’,dr’;, being given the position 
of the particle A, from A,, is equal to the product 
of the following two probabilities: 


(a) The probability W(s’,r’) dr’, dr’, dr’; for th 
particle A, to be situated in the region between r; 


| and r’,+dr’, from Ag. 


(b) The probability W(s—s’,r—r’) dr’, dr’,dr’ ; for 
the particle Ay to be situated in the region between 
r,—r’,and r,—(r’,+dr’, from A,, with A,A,= 
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his latter probability, argued from a symmetrical 
rsion of the chain, assures the condition that the 
icle A, is situated at r; from A,. By taking 
product of (a) by (b), we have 


} dr’ 
W(s 


constant r’) W(s’,r’) dr’, 
This expression becomes identical to (20) after de- 
ination of the constant by the condition of 
nalization (18). 
, the study of intramolecular interactions of long 
n molecules, the above retrograde dispersion 
tion is in general not sufficient. We have to 
know, in general, the probability of the distance 
ween an arbitrary pair of particles, say A; and A, 
which none of them are ends of the chain. In the 
following we shall show how such a more complicated 
retrograde dispersion function can be constructed on 
the basis of the simple one obtained in (20). 
Let us now consider a chain A,A;A,A,. Being 
viven the distance r; of the particle A, from the 
particle A,, we ask the probability 


p'(A;Aj,h; s,r) dh, ... 


the particle A, to be situated in the position 
between Ah, and hy+dh, from A,;. It is the product 
f the following probabilities: 

a) Both ends 7 and / must be subjected to certain 
onditions, namely A, must be in the position between 

h, and Ayr dx h, from A, with probability 


W,A\—h) dr, 


b) A; must be in the position between r;,— A, and 


dd, from A, with the probability 


W(s—l,r—d) dd, 


The product of the above probabilities integrated 
th respect to \ gives 


A A, h; $,T 


’ 


onst. | M4, .. WG, W(s—lr—n), 


pil ~jh; lr Pp ld: sy 


(aa, 


-@ 


sing (17) and (19), 


W (l—j, h) 
W (s,r) 


(lL—47,h; s,r) » using (6). (21) 
Jy", 4 


(17), the member is 


Hence 


\ccording to right-hand 


“9, h: s.r). 
p’(AyAi, h; 8,7)=p'(—j,h; 8,7). 


This shows that the retrograde dispersion function 
for an intermediary pair of particles A; and A, is 
equal to the retrograde dispersion function for the 
end pair of particles A;, and A,. Thus in the 
following section, by studying the nearest neighbor- 
hood dispersion function, we can consider only the 
latter dispersion function for the end pair of particles, 
by putting s=/—j and h=r’. 


[It will be useful to calculate the p’-function of 
(22) from (21) for a W-function of the Gaussian type: 


(<) 7 “Pdr, ees 


with «?=(3/2b")(1/s), and b=length of one link. 
The use of (19) and of polar coordinates transforms 
(21) into the following formula after integration with 
respect to the polar angles 


W(s,r) dr; 


(23) 


p' (AjsA,,h; 8,r) 4rh?dh=p'(l—j, h; s,r) 4ah?dh 


wh, 
i 
Vr 


pr(h 


where 


3 4 

Dh? (24a) 
This formula is an extended form of a formula ob- 
tained by Kuhn, Kiinzle and Katchalsky (1949) for 
the special case of 7=0. 


IV. Nearest Neighbor Probabilities 


In the study of the mechanism of interaction, we 
have to know not only the distribution of the distance 
between two arbitrary particles, but also often such 
a distribution as they become nearest neighbors. 
For discrete particles that are totally independent 
from each other, Herz (1909) has obtained a formula 
for the nearest neighbor distribution functions in 
terms of the distribution function. Applications 
have been made by Jaffé (1940) to the study of inter- 
actions of charged particles. Applications to linked 
particles, such as is the case with a chain molecule, 
has been made by Simha (1948), in connection with 
the problem of volume effect. However, such a 
direct application offHerz’s formula to the density 
distribution of linked particles yields only an ap- 
proximate result. In fact the different particles 
on the chain behave differently as nearest neighbors, 
according to their order of linkage, in such a way 
that the pure randomness of the particles, lying in 
the indifference of order, will be distorted. It is 
the purpose of the following lines to study the 
nearest neighbor probabilities for such mutually 





dependent particles Herz’s method will be extended | 
to the calculations of the nearest neighbor function 
of the retrograde dispersion function 

Consider again a chain OA, A, (s>s’) with origin 
of coordinates at O, the case of spherical symmetry 
being considered. Being given the scalar distance r 
of the particle A, from the origin, the probability for 
the particle Ay to be situated in the region between 
r’ and r’+dr’ is 


“dr’, 


tr 


where p’ is the retrograde dispersion function. Its 
nearest neighbor function will be designated by p’, 
It is the propability for the particle Ay to be situated 
in the region between r’ and r’+dr’, nearest 
neighbor to the particle i), being given the distance 
OA,=r of the particle A It can be calculated by 
taking the product of the following two probabilities 


as a 


a) the probability that the particle Ay is situated 
and r’+-dr’ from O, being 
r for the particle A,; 

b) the probability that 
situated in the region between 0 and r’ 


at a distance between r’ 
given the distance OA 
neighbor is 


no nearest 


The probability indicated in (a) is given by (25 


but owing to the fact that A, cannot be situated in 
the region between 0 and r’, the probability in (a 
will be increased to the value 


| dr’ 4ar” p’ ): 


On the other hand, the probability of finding any 
one of the s particles as a nearest neighbor to the 
particle O in the region dr’ is 


tar? dr’ 


say denoted by o 4ar’*dr’ and the probability that no 
nearest neighbor is situated in the region between 0 
and r’ is 


27a 


Hence by taking the product of the probabilities 
in (26) and (27a) and using (27), we obtain the fol- 
lowing integral equation for o 


a4nr’ dr’ 


In order to simplify the writing, we put? 
pill: I 


4 


then (28 


An integration gives the solution under one of 
following forms 


om 4ar” 


Oa 


It can be shown that o satisfies the condition of 
normalization We write from (29a 


‘ R 


* , 1 
ao 4qnr’ 7 ed 
=) 


By integrating with respect to r’, we have 


This condition means that one of the s particles at 
least must be a nearest neighbor to the particle ( 
in the totality of the domain. 

The function p’; gives the probability of the dis- 
tances between two nearest neighbor particles that 
are linked on a chain. If the particles are not linked 
and behave like gas molecules, we can define a dis- 
tribution W(r) of the distances r_ between 
particles, and a distribution W,(r) when the part 
are nearest neighbors. For this latter case, 
nearest neighbor distribution W, can be obtai 


7 (*) means 


d 
dr’ 





(29) by changing p’, p’; into W7/s, W;/s, re- 
tively, W and W, being now independent of s’. 
or discrete particles formula (29) degenerates to 
formula of Herz (1909): 


Wi(r) Wo) [1-2 far ders w] , (31) 


re s is the number of such independent particles. 

n the study of intramolecular interactions, 7. ¢. 

interactions between the different parts of one 

nt chain, the configurations of any particle of the 

n are wholly controlled by the dispersion func- 

s p, p’, and p’, as studied above. On the other 

d, in the study of intermolecular interactions, 

the interactions between two particles belonging 

two different chains, we might ask: what is the 

‘bability for the particle A,, belonging to one 

in, say chain B, to be situated at a distance 

tween 7” and r’+dr’ from the origin O, at which 

tance it becomes the farthest one of all particles 

B, and therefore becomes the most exposed to 

ractions with surrounding chains? Such prob- 

ibilities, as will be indicated by the subscript ‘2”’, 

e. g. Wo, p's, will be called “farthest neighbor 

probabilities’, in contrast with nearest neighbor 

probabilities. The farthest neighbor probabilities 

an\be calculated in the same way as for the near- 

est neighbor probabilities, except that the limits 

0.r’) have to be changed into (7’,“). Following 
29a) the result is 


p’ 
dr’ 4ar” p’ 


V. Intermolecular Interactions 


The intramolecular interaction, ¢. ¢. between the 
different parts of one giant chain molecule, is com- 
parable to a kind of internal repulsion within the 
hain molecule. Therefore an expansion of the chain 
can be expected. On the other hand, the inter- 
molecular interaction, i. ¢. between a number of 
separate molecules, is a kind of external repulsion 
exerting on every molecule, and will cause a com- 
pression of the chain. As indicated by Simha, the 
observed anormality of the variations of viscosity 
with concentration in connection with the solvent 
fect could point to a shrinkage of two coil-mole- 
cules on approaching each other at concentrations of 
a few tenths of one percent. 

The perturbance on the distribution function 
W(r), by those interactions, can be calculated on the 
basis of the retrograde dispersion function, nearest 
neighbor and farthest neighbor functions derived 
above. The perturbance by intermolecular inter- 
actions will be discussed here for the following very 
clomentary case. 


Let us consider the one-dimensional interaction of 
two chains A and B, which are supposed to move only 
along the direction Oy. We suppose that the chain A 
is fixed at the end A,, and that the chain BP is fixed 
at the end B,, the point A, being taken as the origin 
of coordinates. The particle of B, which forms the 
volume hindrance to A, is the farthest particle of B, 
of coordinate c. Of course ¢ is a random variable 
having a certain distribution that can be calculated 
from the considerations given in section IV concern- 
ing the farthest neighbor probability function. This 
forms a difficult detail in which we shall not enter 
here, but assume only that we could find an effective 
fixed position ¢ forming a fixed obstacle to the motion 
of A. Then the distribution of any particle A,, of 
the chain A will be modified and becomes: 

W,,(8",y’) = W(s" xy’) + W(s8" 2c’ 
where W(s’,y’)dy’ is the distribution free from 
obstacle, W(s’ ,2e—y’) dy’ is the additional distribution 
reflected back by the obstacle, and W(s’,y’) dy’ is the 
total distribution. 

It is easy to verify that the total distribution sat- 
isfies the normalization condition 


| dy’ W,(s’,y’)=1 


if the free distribution is normalized to unity 


(" dy’ W 


Let us define 
dy y’ W(s',y’); 


and 
dy’y' W, s’,y'); yy dy’ y’? W,(s’,y’) 


as the mean distance and the mean square distance 

between the particles A,, A, of the chain A, in the 

free configuration and in the restricted configuration, 

respectively. With the use of (33), and for a free 
distribution 

Wis’ .y')dy’ i 

\ 7 


oy" dy’ 


with «=(2s6*)~', b=length of one link, it is found 


that: 


er? xc): 





The value of y’,>< y” indicates a compression of the VI. References 
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Heat of Combustion and Formation of Cyanogen’ 
John W. Knowlton * and Edward J. Prosen 


Calorimetric measurements were made of the heat of combustion of cyanogen in oxygen 


in a flame at constant pressure. 
combustion: 


C,N, (gas) 


AHc® at 25° ¢ 


+ 20, (gas) = 2CO, (gas) + 


t= — 1095.97 4 


The data yielded the following value for the reaction of 


Nz (gas), 
1.80 kj/mole 


261.94+0.43 kcal/mole. 


Utilizing the value for the heat of formation of carbon dioxide, the following value is 


given for the reaction of formation: 


2C (solid, graphite) ; 


>* 


AHf?® at 2! 


I. Introduction 


This investigation is a continuation of the program 
of determining the heats of formation of chemical 
substances. The heat of combustion of cyanogen 
(gas) Was measured because of the discordance in the 
values for this property given in the literature. 


II. Method and Apparatus 


The calorimetric method and apparatus used in 
the present investigation have been previously 
deseribed [4, 5].2 Platinum resistance thermometer 
No. 262, 214 and resistance bridge No. 404 were used 
[6]. White double potentiometer No. 156,585 was 
used in the electrical energy measurements [4]. The 
same reaction vessel that was used in determining 
the heat of formation of water [4] was used in this 


work. 


The unit of energy upon which values reported in 
this paper are based is the absolute joule, derived 
from mean solar seconds and absolute ohms and 
volts, in terms of which certification of standard 
resistors and cells is made by this Bureau. For con- 
version to the conventional thermochemical calorie, 
ee following relation is used [1, 2]: 1 cal=4.1840 
as ). 


All molecular weights used in this paper were cal- 
culated from the 1947 table of international! atomic 
weights [3]. 


This investigation was completed in 1941, but the preparation of the results 
ublication was delayed by the transfer of the senior author to work in the 
Bureau of Agricultural and Engineering Chemistry and the U. 8, Bureau of 

es. This work formed the basis of a thesis submitted by the senior author to 

Faculty of the Graduate School of the University of Maryland in partial ful- 

ent of the requirements for the degree of Doctor of Philosophy. 

deceased 


gures in brackets indicate the literature references at the end of this paper, 


f 


N; (gas) =C 


oNo (gas), 
C= 308.94+ 1.80 kj/mole 


73.84 + 0.43 keal/mole. 


III. Chemical Procedure 
1. Preparation and Purity of the Materials 


The cyanogen was prepared by a method similar 
to that of Cook and Robinson [7] as described by 
Ruehrwein and Giauque [8]. This method con- 
sisted of precipitating silver cyanide from a saturated 
aqueous solution of silver nitrate by addition of a 
slight excess of a 78-percent aqueous solution of 
potassium cyanide. The silver cyanide was filtered, 
washed free of excess potassium cyanide, and dried 
in an oven at 130° C for 15 hrs. It was then trans- 
ferred to a 3-liter flask, evacuated to 10-* mm of 
mercury and then slowly heated to 380°C to drive 
off cyanogen. About 350 ml of liquid cyanogen was 
collected. 

The sample was degassed and then distilled in a 
l-m vacuum-jacketed distillation column packed 
with glass helices. The condenser temperature was 
kept constant at —24.0°C. A reflux ratio of 10 to | 
was used and the take-off rate maintained at 30 ml/hr. 
The first fraction of 50 ml and residue of 45 ml were 
discarded and the two middle fractions of 120 ml 
each were combined and stored in a 1-liter brass 
bottle for use in this investigation. Three 5-ml 
samples, taken after the first fraction, between the 
middle fractions and after the middle fractions, were 
used to get an indication of the purity of the cyanogen. 
Differential vapor-pressure measurements were made 
on these 5-ml samples according to the method of 
Shepherd [8], but not using the micrometric-reading 
devices used by him. The observed differences in 
vapor pressures of these samples taken over the 
temperature range of —24.8° to —23.3° C were 0.0 
+0.5 mm, indicating that the sample was probably 
sufficiently pure for this investigation. The cyanogen 
prepared by Ruehrwein and Giauque [9], utilizing 
this same procedure except that they passed their 
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cyanogen through phosphorus pentoxide as it was 
generated, was found by them to have a purity of 
99.995 mole percent 

To test the effect of cyanogen on the materials 
with which it would be in contact, a small amount of 
cyanogen was placed in a glass bulb containing 
pieces of copper, brass, and steel as well as a smear of 
stopcock grease. No effect on the materials was 
observed after 10 months of storage. 

The oxygen (prepared commercially from liquid 
air) used for the combustion was freed of combustible 
impurities by passage through a tube containing 
copper oxide at 600°C. It was then freed of carbon 
dioxide and water by passage successively through 
tubes containing ascarite, magnesium perchlorate, 
and phosphorus pentoxide 


2. Purity of the Combustion Reaction 


The cyanogen was dried by passage through 
anhydrous calcium sulfate and was burned in an 
atmosphere of the purified oxygen in the calorimeter 
reaction The combustion was initiated by 
passing an electric spark across the ends of platinum 
wires at the tip of the burner. The flame burned 
quietly in the reaction chamber 

The main products of combustion of cyanogen 
are carbon dioxide and nitrogen. Other possible 
products are the oxides of nitrogen, carbon monoxide, 
carbon, and unburned cyanogen 

No significant deposit of carbon was observed in 
the reaction vessel after about 50 ignitions of the 
eyanogen, with combustions of varying duration, 
including those of the calorimetric observations. 

Because nitrous oxide itself supports combustion 
and since the combustion of cyanogen occurred in a 
flame in the presence of about 50 percent excess 
oxygen, the possibility of the presence of nitrous 
oxide as a combustion product appeared to be very 
remote. The other oxides of nitrogen are acidic, and, 
mixed with oxygen, are completely absorbed in 
alkaline solutions. Nitric oxide will oxidize under 
these conditions [10] Tests were made by passing 
the effluent gas from the combustion of cyanogen in 
the calorimetric reaction vessel through a solution of 
sodium hydroxide and determining (see below) the 
total nitrogen in the solution. Fixed nitrogen was 
found in amounts varying up to 2 percent of the total 
nitrogen in the cyanogen burned. Attempts to 
reduce this amount by varying the flame velocity of 
the combustion, by changing the technic of ignition, 
or by use of the reaction vessel used in the experi- 
ments on ethane [11] were unsuccessful, or 
difficulty in securing ignition 

Saturation of the oxygen used for combustion 
with water vapor at 0° C, as was done by von 
Wartenberg and Schutza [12], reduced the amount of 
fixed nitrogen to about 80 percent of its previous 
amount, and this technique was used in one of the 
series of calorimetric combustion experiments re- 
ported here. Since these amounts of nitrogen 
oxides were appreciable, it was necessary to deter- 
mine the amounts quantitatively for each calori- 
metric experiment 


vessel 


gave 


In the actual calorimetric combustion ey jer;- 
ments, the effluent gas from the calorimetric rea. | jo) 
vessel was passed through a U-tube packed .ye¢- 
cessively with ascarite, magnesium perchlorate, and 
phosphorus pentoxide, separated by layers of dry 
acid-washed fibers. This tube absorbed 
all of the carbon dioxide, water, oxides of 1 tpo- 
gen, and unburned cyanogen issuing from the ca 
orimeter. The complete absorption of the oxides 
of nitrogen was verified in a test combustion ex- 
periment of cyanogen run for 4 min under the 
conditions as in the calorimetric experiments. The 
gas issuing from the ascarite tube was bubbled into a 
10-percent sodium hydroxide solution. No fixed 
nitrogen was detected in this solution in a test that 
would have shown the presence of less than 0.01 mg 
The complete absorption of any unburned cyanogen 
was verified by passing cyanegen directly into the 
ascarite tube. No effluent gas was observed wher 
much as 0.04 g of cyanogen was passed into the tub 

The determination of the products of combustion 
was made as follows: After a combustion, the as- 
carite tube was weighed, the stopcocks were re- 
moved, and the stopcock grease carefully removed 
with ether. The layers of phosphorus pentoxid 
and magnesium perchlorate were removed, and th 
traces adhering to the glass were washed out wit! 
a jet of distilled water in a manner such as to avoid 
wetting the ascarite. The asbestos plug over th 
ascarite was removed and the ascarite dissolved i 
water. This solution was filtered and diluted to 
200 ml in a volumetric flask. From this solutior 
100 ml was pipetted out for the nitrite determinatiot 
and the remainder was used for the determinatio! 
of total nitrogen. The amount of nitrate was 
obtained as the difference between the total nitroge: 
and the nitrite. Thus, by determining both nitrit’ 
and nitrate, more information is obtained on the 
oxides of nitrogen formed than if only total nitroger 
had been determined 

For the titration of nitrite, 30 ml of 0.017 mola 
potassium permanganate was added from a burett 
to an Erlenmeyer flask. This was diluted to 75 m 
with water, and 25 ml of 18 N sulfuric acid was added 
The 100-ml sample was added to this with the flask 
cooled in tap water. This cold solution was the! 
treated with 0.04-molar sodium oxalate until de- 
colorized, a definite amount of oxalate being added 
The end point of the titration reached by 
further addition of the permanganate, with th 
solution heated to 50° to 60° C before addition of 
the last few For 
permanganate equivalent of the oxalate was 
termined in the manner. The 
used as the primary standard, with 1 mole of oxalat« 
equivalent to 1 mole of nitrite. 

For determination of total fixed nitrogen, the re- 
maining alkaline from the 
reduced for 1 hr with 2 ¢g of Devarda’s alloy in 4 
Kjeldahl still, and then the ammonia was dist: le 
into 50 ml of 0.05 N sulfuric acid, which was 
titrated with 0.05 N sodium hydroxide 


asbestos 


was 


drops each cetermination th 


same oxalate was 


solution ascarite 1s 
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ri- ontrol tests were performed in which known | 3. Determination of the Amount of Reaction 
ion a .ounts of sodium nitrite and sodium nitrate were | 
in. a ‘ded to the asearite. The results indicated that The amount of reaction was determined from the 
nd the method was sufficiently accurate for this in- | mass of carbon dioxide formed in the combustion, 
rv, vy stigation, To the increase in mass of the first ascarite absorp- 
red it was also shown in a control test that unburned | tion tube connected directly to the outlet of the 
ro- cyanogen absorbed by the ascarite had no effect on | reaction vessel (which absorbed all the carbon diox- 
al- the fixed nitrogen or nitrite determinations. ide, water, oxides of nitrogen, and unburned cyano- 
des l'o determine any carbon monoxide formed in the | gen issuing from the calorimeter), there was added 
eX- combustion, the effluent gas from the ascarite tube | the increase in mass of the second ascarite tube, 
me was passed through a tube packed with copper | which was connected to the first through the copper 
The oxide and maintained at about 600° C, then through | oxide combustion tube (and which absorbed the 
Oa a second ascarite tube packed in the same manner as | carbon dioxide formed by oxidation of the carbon 
xed the first asearite tube. In this way, any carbon | monoxide). From this were subtracted the masses 
hat monoxide formed in the combustion was oxidized | of nitrogen oxides as determined from the nitrogen 
ng to carbon dioxide by passage through the hot copper | analyses (present as nitrite and nitrate and calcu- 
gen oxide tube and was weighed as carbon dioxide in the | lated as nitrogen trioxide, N,O;, and nitrogen pent- 
the second ascarite tube. It has been previously shown | oxide, N,O;, respectively), the average mass of un- 
| as that no oxides of nitrogen nor unburned cyanogen | burned cyanogen (0.01730 g), and the average mass 
ibe passed through the first ascarite tube. Also, tests | of water introduced in the experiments using water- 
ion for nitrogen on the ascarite in the second tube gave | carrying oxygen. This result gave the total mass of 
as- negative results. In the series of calorimetric com- | carbon dioxide representing the amount of reaction. 
re- bustion experiments in which dry oxygen was used | One mole (44.010 g) of carbon dioxide was taken as 
ved for combustion, it was found that an average of | equivalent to }; mole of cyanogen. 
‘ide 3.5 percent of the total carbon in the cyanogen The average mass of water introduced in the 
the burned to carbon monoxide, In the series of cal- experiments in which the oxygen used for combus- 
ith orimetric combustion experiments in which water- | tion was saturated with water at 0° C was determined 
oid carrying oxygen was used for combustion, an average | jn 10 separate experiments. In these experiments, 
the of 0.35 percent was formed. rhe amount of carbon water-carrying Oxygen was passed through the reac- 
| in monoxide was determined in the above manner for | tion vessel and absorption tube at the same rate and 
to each calorimetric combustion experiment. for the same length of time (46 min) as was used in 
ion lt was to be expected that some cyanogen would | the average combustion experiment. The average 
on pass through the calorimetric reaction vessel un- | amount of water absorbed was 0.05234 g (standard 
ion burned during the process of ignition and extinction | deviation 0.00137 g). An uncertainty of 0.0014 g 
was of the flame. This amount of cyanogen had to be | of water would cause an uncertainty of +0.029 
gel determined since if would be absorbed In the first | percent in the determination of the amount of reac- 
rit asearite tube and the increase in mass of this tube | tion by the above method. 
the was used to determine the amount of reaction (see 
gen below). This amount of cyanogen was deduced | 
from a series of eight separate experiments carried | 
lar out calorimetrically under the same conditions as in IV. Calorimetric Procedure 
tte the main calorimetric combustion experiments, but 
ml with extinction of the flame immediately following : 
led ignition. The average amount of cyanogen burned 1. Electrical Energy Experiments 
ask in these experiments was 0.0077 g. The increase in 
hen mass of the ascarite tube was considered to be due The results of the experiments to determine the 
de- to absorption of carbon dioxide and unburned cyan- | electrical energy equivalent of the calorimeter system 
led ogen. From the temperature rise of the calorimeter, | ®ve given in table 1. There are listed: AR,., which 
hy together with its electrical energy equivalent, the | is the increase in_ temperature of the calorimeter 
the total energy evolved in each experiment was caleu- | System corrected for thermal leakage and heat of 
“ar lated. From this total energy was subtracted the | Stirring [4] and expressed as an increase In resistance 
the energy of sparking (see below) to give the energy | in ohms of the given platinum thermometer, as meas- 
due to combustion of evanogen. The mass of carbon | ured with the given resistance bridge; the total elec- 
de- dioxide, corresponding to this quantity of energy, | trical energy added to the calorimeter system ex- 
was which would have been formed in the pure combus- | pressed in absolute joules; and the electrical energy 
ate tion of cyanogen was calculated, using the approxi- | equivalent of the calorimeter system, corrected to the 
mate value of 1096 kj/mole for the heat of combus- | standard amount (3600.00 g) of water in the calo- 
re- tion of cyanogen. This mass of carbon dioxide was | rimeter can. The average temperature of the 
vas subtracted from the increase in mass of the ascarite | calorimeter was 25.00° +0.01° C in all experiments. 
1a tube to give the mass of unburned cyanogen. The | The voltage across the heater was about 40 v, the 
ed average mass of unburned cyanogen was found to be | resistance of the heater was about 60 ohms, and the 
L. 0.01730 ¢ (standard deviation of 0.00013 g). time of heating was 2280.00 see in each experiment, 
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TAaBLe 1 Results of the electrical calibration experiments 


Deviation 
from 
mean 


ohm 
0.0 
0 

i) 

l 

l 


7.8 


Mean 
Standard 


calorimeter car 


2. Correction Experiments 


The spark energy was determined calorimetrically 
in five separate experiments in which only sparking 
energy was supplied to the calorimeter with oxygen 
in the reaction The average value of the 
sparking energy was found to be 2.07 j/see (standard 
deviation 0.02 j sec 


\ esse] 


TABLI 


2 Results of the 


Spark 
nergy 


Net reactior on 
ma 
energy * " 


r deviatio from the standard calorimeter syst 
imount of nitrite determined 
amount of nitrate determined 


rage mass (0.01730 g) of unburned cyanogen 


* Corrected fi 
> Calculated from the 
* Calculated from the 


Corrected for the ave 


ibsorbed j 


TABLE 3. Results of the calor combustion experiments on ¢ 


t reaction 
nergy * 


ations from tl! dard calorimeter system 

1 he amount ¢ te determined 

m the amount termined 

the average mass of 1 ur i cyanogen and of water al 


orbed | 


The gas energy, column 3 of table 2, was 
5 em. 


lated from the known heat capacities of the 
and the difference between their temperature 
as that of the air near the tubes leading to t! 
action vessel) ‘and the average temperature, 
respect to time, of the calorimeter. For this 


pose, the heat capacities of oxygen and cya 
were taken as 29.2 and 44.6 j/degree-mole, r 
tively [13]. 

No correction was made for the heat capacit 
the effluent gases from the calorimeter, since 
average temperature was substantially 


25° C. 


3. Combustion Experiments 


The data of ten calorimetric combustion ex 
ments on cyanogen using dry oxygen for combu 
are given in table 2. The net reaction energy is th: 
energy evolved in the experiment after subtracting 
the gas energy and spark energy. The energy evolved 
was calculated from the corrected temperature ris¢ 
AR., combined with the electrical energy equivalent 
of the standard calorimeter system. Correction was 
made for deviations from the standard calorimeter 
system due to small variations in the amount of 
water in the calorimeter can in the individual ex 
periments. 


calorimetric combustion experiments on cyanogen, using dry orygen for combustion 


yanogen, 


Mass of | Mas 
NOs rot 


0. 01884 
01086 
O1161 


O2188 


0. 00602 
OT675 
CTHH6 
09560 
10390 
09928 
10204 


01802 
02536 


O1835 
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ie data of seven calorimetric combustion experi- | 
is on cyanogen, using oxygen saturated with 
r at 0° C, are given in table 3. In those ex- | 
nents for which the time of flow of the water- 
ving oxygen differed from 46 min, the amount 
ater absorbed was calculated from the average 
e of 0.05234 g for 46 min and assuming the 
unt to vary linearly with the time. 


V. Results 


lor the final calculation of the heat-of combustion 
of cyanogen, it is necessary to identify the particular 
oxides of nitrogen formed in the colorimeter during | 
the combustion in order that appropriate corrections | 
for their respective heats of formation may be ap- 
plied to the calorimetric observations. For this 
purpose there must be considered nitric oxide, nitro- | 
gen dioxide (in equilibrium with its dimer, nitrogen 
tetroxide), nitrogen trioxide, and nitrogen pentoxide, 
as well as nitric acid vapor in the case of the experi- 
ments in which water vapor was present in the | 
reaction chamber. 

Nitrogen pentoxide is known to decompose com- 
pletely at 25° C [10, p. 552]. The dissociation of | 
nitrogen trioxide into nitric oxide and nitrogen di- | 
oxide is also virtually complete [10, p. 449] under 
the conditions of the present experiments. This is 
corroborated by a calculation based on the equilibri- | 
um data of Verhoek and Daniels [14]. The fraction 
of nitrogen tetroxide in equilibrium with nitrogen 
dioxide, computed from the data of Giauque and 
Kemp [15], was found to be negligible under the 
conditions of the present experiments. 

It thus appears that the oxides of nitrogen in the 
effluent gas from the reaction vessel were nitric 
oxide and nitrogen dioxide. Concerning the absorp- 
tion by solid alkali (or ascarite) of these substances, 
Mellor reports the following: (1) nitric oxide reacts 
with solid potassium hydroxide in the presence of 
excess oxygen, to form mainly potassium nitrite, but 
that some nitrate is formed as well [10, p. 432]; 
2) nitrogen dioxide should react with solid alkali to 
give equimolal amounts of nitrite and nitrate, but 
that larger amounts of nitrate may be formed 
10, p. 5389-40]; and (3) nitric oxide and nitrogen 
lioxide in equimolal amounts are absorbed by alkali 
as pure nitrite. Burdick and Freed [16] state that 
the assumption that these absorptions of nitrogen 
xides are pure reactions is subject to some qualifica- 
cams 

A reasonable interpretation of the process occur- 

ng in the calorimeter with respect to the fixation 
{ nitrogen during the combustion of cyanogen, is 
hat both nitric oxide and nitrogen dioxide are 
rmed in the amounts z moles and y moles, re- 
spectively; 


(1) | 


(2) | 


r[% No(g) + O.(g)=NO()I, 


yl No(g)+O2(9)=NO,(g)]. 


The subsequent reaction of these oxides with ascarite 
is assumed to take place according to the following 
reactions 


2z{|NO(g) + NO,(g) + 
2NaQH(s)=2NaNO,(s) + H,O(s)], (3) 


y—xz{|NO,(g) + NaOH (s) = %NaNO,(s) + sNaNO,(s) 


Letting @ and b represent the number of moles of 
nitrite and nitrate determined, respectively, then the 
amounts of nitric oxide and nitrogen dioxide formed 
are 

W(a—b), (5) 


y=(a+3b). (6) 


In the series of experiments in which water vapor 
vas present in the reaction chamber, it was assumed 
that nitric oxide, nitrogen dioxide, and nitric acid 
vapor were formed and that the nitric oxide and 
nitrogen dioxide were present in equimolal quantities 
z’ and the nitric acid in quantity y’: 

NO(g)], (7) 


z’[’No(g) + 4O2(9) 


x’ [4%N2(g) +-O2(g) = NO,(g)], (8) 


y’(4N2(g) +%O2(g) +4H,0(g)=HNO,(g)]. (9) 


) 


The subsequent reactions with ascarite are assumed 
to be 


+-2NaQOH(s) 
H,O(s)}, 


zx’ [NOwg) + NO,(g) 2NaNO,(s) 4 


(10) 


y’ |[HNO,(g) +- NaOH (s) = NaNO, (s) +H,O(s)]. (11) 
Thus, the amounts of nitric oxide, 2’, nitrogen di- 
oxide, x’, and nitric acid vapor, y’, formed are 


Va (12) 
ty’ =, (13) 


where a and 6 represent the number of moles of 
nitrite and nitrate determined, respectively. 

The calorimetric data given in tables 2 and 3 
were reduced on the basis of these assumptions and 
are presented in tables 4 and 5, respectively. 

The heats of reactions (1), (2), and (9) were 
taken as AH/°=90.37, 33.85, and —12.85 kj/mole [2] 
and the heat of combustion of carbon monoxide 
was taken as A//c°=— 282.99 kj/mole [2]. 

It is recognized that, in this work, there remains 
an uncertainty in the actual composition of the 
effluent gas from the calorimeter. Cognizance of 
this was taken in estimating the over-all uncertainty 
in the heat of combustion. 
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Tasre 4. Calculation of the heat of combustion of cyanogen from the data of table 2 in which dry oxygen was used for coml 


Amount Cc tior Corres Corres Corrected {mount Heat of Devia 
of CO orrection tion tion reaction of CoN combustion tion from 
for NO - : 


formed for NO for CO energy burned ut 25” ¢ mean 


Moles kj/ mole ki/ mole 
0. 0542553 1093. 62 1. 38 
Zi 0542108 1095. 08 +0). O05 
OO2U9R2 OS38550 1005. 94 + O4 
OOSSIRS Ls i7 GOS v 0542215 10906. 13 +1. 1 


0040268 O545178 1005. 22 +(), 22 


i iv 


004144 % 73 ry 5 0545066 1004. 8&2 18 
0542280 1006. 5 + 5 
0038404 O47360 1095 + 11 
OO44110 5 My ; 587 7 0547630 1093. 3 As 
OO44335 7 § ; SOS. O} 46823 1094. 2 st) 


Os 2208 


1095. 00 
+O). 34 


trite and nitrate, assuming that the absorption of nitrogen oxides proceeded as indicated by eq 
I 


TABLE 5 Calculation of the heat of combustion of ¢ yanogen from the data of table 3 in which oxygen saturated with wate 


C’ was used for combustion 


Amount Amount Corree Correc Corre Correc Corrected Amount Heat of De 
of HNO of CO tion tion tion for tion reaction of CyN combustion fr = 
0 

5° ( 


formed* formed for NO for NO HNO for CO energy burned it 2 


Moles Moles Moles kji/ mole 
0. OO05488 0. 0008540 5 § 0. 0534423 1007. 85 
OOONATT 002990, O534804 10v6 
oOoo21 50 o005440 0534100 1005 
OOO40451 OO0S3876 Sug 1097 
OOO8337 ooosl22 1007. 77 
0004806 ooos479 1097 
OOOS398 OO00S 108 1097 


nitrate, assuming 1 absorption of 1 id procees s indicated t 


The heat of combustion of cyanogen, taken as The heat of combustion as determined in this 
the mean of the two series and including a correction | investigation is compared with the following existing 


of 0.00, kj to adjust the value to the standard state, | data: 
be Omes 
N ‘vanogen, gas) +: ; Oy | 
am & . , Year 


Investigation 
reporte 


Alle 2! } 1095.97 1.80 kj/mole 


261.94 0.43 keal/mole - 
bos Dulong [17] 1838 


4 Berthelot [18] 1879 
The assigned uncertainty is the estimated over-all = , J 
Thomsen [19] 1SS6 


uncertainty considering the uncertainty in the exact vam Wastenbers & Gehuten 
products of combustion, in the amount of water [12] 1933 
earried through the calorimeter by the oxygen, in MeNorris & Badger [20] 1933 
the electrical energy equivalent, and in other factors Present oan 
considered to be of lesser importance 
Combination of this value with the value for the 

heat of formation of carbon dioxide A// 393.513 
kj/mole [2] vields the heat of formation of cvanogen 


2C (ec, graphite)+N, (gas C,N, (cvanogen, gas), T 
ic 


(15 encouragement and assistance of M. 


Haring, formerly of the University of Maryland 
of F. D. Rossini, who directed this work 
73.84 +0.43 keal/mole gratefully acknowledged 


AHf? (25° C) =308.94 1.80 kj/mole 
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Researcl 


Behavior of Bromophthalein Magenta E (Tetrabromophe-. 
nolphthalein Ethyl Ester) With Organic Bases and Its 
Bearing on the Brgnsted-Lowry and Lewis Concepts 
of Acidity 


Marion Maclean Davis and Hannah B. Hetzer 


In a continuation of spectrophotometric studies of its indicator properties in benzene, 


bromophthalein magenta E 
react with 1,3-diphenylguanidine and 
correspond to the equations: l R N- 

R,N=(R NHNR xX blue 
with diphenviguanidine are 2.2x 10° and 
phenylguanidine, A, is 
phenylguanidine 


as increasing the 


tetrabromophenolphthalein ethyl ester 
with 
HX= 
The association constants A 
15.5, respectively 
525 and kK, is of the same order as the K, for the reaction with di- 
Re placing benzene by ethanol or py ridine produces the same kind of effect 
concentration of di- or triphenylguanidine. 


has been shown to 
1,2,3-triphenylguanidine in two steps that 
R,;NH+... X-(magenta); (2) R;5NH+.. .-X 
and A, for these reactions 
For the reactions with tri- 


The results indicate that the 


Brénsted-Lowry formulation of an acid-base reaction should be modified, and that hydrogen 


acids should be regarded as a special class of 


I. Introduction 


In connection with a project at this Bureau to 
develop methods for measuring acidity and basicity 
in organic solvents, the utility of indicator dyes has 
received special consideration. Indicator dyes, if 
available, would be particularly useful for solutions 
of acids and in inert solvents like benzene, 
because such solutions do not contain significant 
quantities of free hydrogen ions, and measurements 
by means of the hydrogen or glass electrodes are 
therefore not applicable. In a previous paper {1],? the 
preparation of the acidic indicator dye, bromophthal- 
ein magenta E (the ethyl ester of tetrabromo- 
phenolphthalein), was described and a comprehensive 
survey of its reactions with organic bases in benzene 
and some other organic solvents was presented. A 
quartz photoelectric spectrophotometer was used in 
following the course of the reactions. This paper 
deals with a continuation of such studies, which are 
& prerequisite to practical applications of the 
indicator 

Apart from the bearing of the results on qualitative 
and quantitative applications of bromophthalein 
magenta E, the latest measurements are of impor- 
tance in clarifying the conflict between the Brgnsted- 
Lowry and Lewis concepts of acidity. As is well 
known, Brgnsted and Lowry formulated an acid-base 
reaction as the transfer of a proton from one base to 
another {2 to 5], and Lewis formulated it as the 
addition of an acid to a base [6 to 8]. Compounds 
that can donate a proton, which are the only sub- 
stances recognized as acids bv Breénsted and Low ry, 
form a homogeneous group for which valuable 
quantitative interrelationships have been derived. 


bases 


American Chemical Society at 
wish to acknowledge a gr 
rk 
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the end 


of this paper 


Lewis acids 


| This has not yet proved true for Lewis acids, which 
include such diverse substances as silver ion, iodin: 
and sulfur dioxide. However, the Lewis theory has 
the merit of correlating a much wider range of 
chemical phenomena in a qualitative sense. Thus 
two different but useful concepts of acidity hav 
coexisted, in much the same way as the wave an 
corpuscular theories of light. 

Over 40 years Lewis said, We ar 
beginning to realize that the phenomena of aqueous 
solutions are but special instances of the widely vary- 
ing phenomena occurring in other solvents; and it 
seems unlikely that a satisfactory understanding of 
the behavior of aqueous solutions can come excep! 
through a careful study of nonaqueous solutions 

i [9]. Two decades later Bronsted wrote, “Th 
properties of acids and bases may then by no means 
as is generally believed, be investigated in aqueous 
solutions but must on the contrary be investigate: 
in solutions of benzene-like, aprotic character if wi 
would gain the most comprehensive and the deepest 
insight into the nature of these substances” [5 
Within the past half-century systematic investig 
tions have been made of acid-base systems in non- 
aqueous solvents such as liquid ammonia, li 
sulfur dioxide, and glacial acetic acid, but apart frot 
a few exploratory studies, outstanding among whic! 
are those by LaMer and Downes (10, 11], investiga- 
tions in aprotic solvents have been neglected. Ever 
Brénsted left the systematic study of acid-bas 
behavior in inert solvents to later investigators 
except for one paper that dealt briefly with rela 
acidic strengths in benzene and included no ex} 
mental data [12]. This neglect is understand 
because of the unavailability of suitable indy 
acids and bases. The lack of indicators still exists 
as pointed out by Kolthoff, who recently w 
“With the extension of acid-base titrations to org 


ago 
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mts the need for more acid-base indicators is 
ving. The titration of Lewis acids in aprotonic 
ents is still awaiting systematic research on 
librium constants between these acids and 
ious bases. Such studies include the search for 
able indicators and their equilibrium constants 
iven systems” {13]. This need will probably not 
met completely until a number of new indicators 
e been synthesized that have adequate solubility 
the solvents concerned, as well as adequate 
sitivity and stability. However, various dyes 
it are described in the literature or known to dye 
nufacturers would probably prove very useful if 
ade available to investigators. For example, the 
mmpound termed by us “bromophthalein magenta 
has been known since 1897 [14], although its po- 
tential use as an indicator of basicity in inert solvents 
was not known until a systematic search at this Bu- 
reau for suitable indicators disclosed its valuable 
properties [1]. 

Before proceeding to a discussion of our latest 
studies, it is well to review briefly some of our earlier 
observations. With regard to bromophthalein ma- 
genta E, it is a brick-red solid that dissolves readily 
in benzene and other inert solvents to give a yellow 
solution. It is a monobasic acid and is monomeric 
in solution; these properties result in greatly simpli- 
fied behavior of the indicator as compared with the 
sulfonephthalein indicators, which are dibasic acids, 
and with carboxylic acids, which are partly dimerized 
in solution. An instantaneous reaction occurs be- 
tween bromophthalein magenta E and an aliphatic 
amine in an inert solvent, as shown by a vivid change 
in color. In contrast to the behavior of an acidic 
indicator with bases in aqueous solutions, the color 
produced is not the same for all bases. Thus, with 
a primary amine such as benzylamine it is red-purple, 
with a secondary amine such as piperidine it is purple- 
blue (“‘cornflower” blue), and with a tertiary amine 
such as triethylamine it is magenta. A magenta 
color is also produced by 1,3-diphenyl-, 1,3-di-o- 
tolyl-, and 1,2,3-triphenylguanidine and by various 
other bases such as hexamethylenetetramine. <A solu- 
tion of a quaternary ammonium salt of bromophthal- 
ein magenta E in an inert solvent is blue, without a 
purple tinge. An explanation of this differentiating 
behavior of bromophthalein magenta E toward differ- 
ent classes of bases was presented in the reference 
cited {1}. 

For the moment, let us focus attention on the 
reaction that produces a magenta color. This 
reaction may be formulated quite simply as the 
combination of one molecule of the indicator (A) 
with one molecule of the base (B) to give a molecule 
of a salt that will be referred to in this paper as the 
primary salt (S,). In the examples studied, the 
reaction did not proceed to completion, but produced 
an equilibrium mixture of A, B, and S,. It therefore 
offered a means of determining the relative strengths 
of bases in terms of association constants. That is, 


[S,] 


“[A} (BT (1) 


} ne 





Provisional values of the association constant were 
given previously for triethylamine, diphenylguani- 
dine, and di-o-tolylguanidine [1]. It was desirable 
to check some of these measurements with better 
control of the temperature than was feasible formerly 
and at additional concentrations made possible by 
the acquisition of shorter absorption cells. This 
paper is concerned with measurements of the reaction 
of bromophthalein magenta E with triphenylguani- 
dine in benzene and an extension of the earlier 
measurements with diphenylguanidine, as well as 
effects produced by adding ethanol or pyridine. As 
will be shown, these measurements have demon- 
strated that salt formation between bromophthalein 
magenta and a base such as diphenylguanidine occurs 
in two steps, yielding products that we shall term 
the primary salt (S,) and the secondary salt (S,). 
The first step would be obscured in an active solvent 
such as water. The theoretical significance of these 
results will be discussed. 

A few words should be said about the character of 
benzene as a solvent. Although classed as an inert 
or “aprotic’’ solvent, it is known to exhibit basic 
properties in some instances. A full discygsion of 
the subject is outside the scope of this paper, but a 
few examples may be cited here. For instance, 
O’Brien, Kenny, and Zuercher obtained evidence of 
hydrogen bonding between benzene and hydrogen 
chloride and gave a provisional value of 0.16 (in 
terms of molalities) for the association constant; by 
comparison, a value of 1.23>10° was given for the 
association constant when hydrogen chloride com- 
bines with water [15]. Staveley, Jeffes and Moy 
discussed evidence from solubility measurements that 
benzene can sometimes associate with other molecules 
that contain a hydrogen atom capable of forming a 
bridge (for example, water), and suggested that the 
hydrogen atom of the donor molecule may be at- 
tracted toward the center of the benzene ring in a 
direction normal to the plane of the ring [16]. 
Recently, Benesi and Hildebrand have obtained 
interesting evidence of 1:1 combination of benzene 
and iodine in carbon tetrachloride, with the asso- 
ciation constant 1.72 [17]. Mulliken has discussed 
various possible structures of the complex of iodine 
and benzene, the most probable of which he considers 
to be a structure in which the iodine molecule lies 
above the plane of the benzene ring, with its axis 
parallel to the plane of the ring [18]. As to possible 
acidic properties of benzene, Rao and Jatkar con- 
cluded that one or two of the hydrogens in benzene 
may participate in hydrogen bridge formation, for 
example, to chlorine atoms of chloroform, methylene 
chloride, or carbon tetrachloride [19]. The available 
evidence may be summarized as indicating associa- 
tive forces between benzene and many other types 
of molecules that account for its excellent solvent 
properties; in some cases, the associations may be 
likened to acid-base interactions. The associative 
forces are relatively very small as compared with 
those involved in the interaction of compounds 
familiarly known as acids and bases, such as the 
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compounds with which this paper is concerned, 
although in some instances it may prove necessary 
to take them into account 


II. Equipment and Procedure 


Spectrophotometric data were obtained with a 
Beckman Model DU quartz photoelectric spectro- 
photometer The absorption cells, which were all 
provided with glass stoppers, were assembled from 
endplates and spacers constructed in the optical 
shop of this Bureau. For 1-cm cells, both endplates 
were of crystalline quartz, 2 mm in thickness, and 
the spacer was made of fused quartz. For shorter 
cells, one of the 2-mm endplates was replaced by 
the tubulated fused-quartz endplate previously de- 
scribed [1]. The spacers used for 1-, 2.5-, and 5-mm 
cells were made of stainless steel. The shortest 
cells were made by combining a tubulated endplate 
with an endplate of crystalline quartz, 3 mm in 
thickness, having a central indented area 0.1 mm 
or 0.5 mm in depth. The cells were placed in a 
thermostated box maintained at 25.0°+0.1° C [20] 
half an hour or more before measurements were 
begun. The temperature of the room was auto- 
matically controlled to stay between 24° and 27° C 
but was usually between 25° and 26° C 

The procedure followed in the preparation and 
dilution of solutions has already been described [1] 
When the molar concentration of the base was twice 
that of the indicator or greater, the same amount of 
base was added to both cells. Otherwise, the 
reference cell contained only the pure solvent. 

The following symbols and terminology are used. 
C., €,, and C,=the initial concentration of the salt, 
acid, and base, respectively, in moles per liter; [S], 
[A], and [B]=the molar concentration of the salt, 
acid, and base, respectively, in an equilibrium mix- 
ture; S,=the primary salt and S,=the secondary salt 
T. (transmittaney of the solute Wcatiell Baotie® 
absorbancy of the solute) logy7.: @y (molar 
absorbancy index)=A,/(6.\4), where b=length in 
centimeters of the absorption cell, and \/=concen- 
tration of the solution in moles per liter 


III. Materials 


Benzene, bromophthalein magenta E, 1,3-diphenyl- 
guanidine, and 1,2,3-triphenylguanidine were of the 
same quality and were purified in the same way as 
the materials described previously [1]. Commercial 
“absolute’’ ethanol was purified by the method of 
Lund and Bjerrum [21] immediately before use. A 
high grade of ps ridine was dried over sodium hy- 
droxide and distilled from fresh sodium hydroxide 
immediately before use. Separate fractions of the 
distillate were collected in a desiccator protected by 
a soda-lime tube; the middle fraction was used. 


IV. Data and Discussion 
1. The Primary Reaction 


As previously reported [1], the spectral absor! 
curve of BPM-—E (bromophthalein magenta | 
benzene shows a maximum near 405 mug, w 
molar absorbancy index of 26,200. 

When diphenylquanidine is added in succe 
small increments to a solution of bromophth 
magenta E in benzene, the original yellow col 
the indicator changes progressively through tor 
orange and red to the magenta color that indi 
formation of the primary salt. The primary 
shows maximum absorption near 540 mu. The 


mediate tones correspond to mixtures of unchar og 
indicator with the primary salt. Curves obtain 


for some of the solutions are shown in figur 
Curve A is for bromophthalein magenta E (2.5) 


M, 5-mm cell Curves 1 to 5 are for its mixtur 
with 0.2, 0.4, 0.6, 0.8, and 1.0 molar equivalent 


respectively, of diphenylguanidine. For the « 
plete conversion of BPM-—E to its magenta salt 
excess of diphenylguanidine was necessary. | 
experiments under discussion, the absorbancy « 
remained practically constant when the molar « 
centration of diphenylguanidine was from 25 t 
times that of the indicator. At much higher 

centrations of diphenylguanidine, the magenta 


! 


shifted toward blue, indicating a secondary reactior 
The best “limiting curve” for the primary reacti 


was considered to be the one that had the lowes 


value of absorbancy at 405 mu and the highest val 
at 540 mu. On this basis, the curve for a solutio 


that contained 2.510 VU of BPM-—E with 


molar equivalents of diphenylguanidine was adopté 
as the curve for the primary salt, S,; the mol 


absorbaney indexes for 405 mu and 540 mu are 2,2 


and 40,800, respectively. The curves in figur 


show sharp isosbestic points near 355 my and 455 mz 

As bromophthalein magenta does not show measu! 
able absorbancy at 540 my under the conditions 
used, this is a very suitable wavelength for quant 


tative studies of the formation of the primary salt 


The sensitivity of the absorbancy at 540 mu to sn 
increments of diphenylguanidine is shown in fig 


The molar absorbancy index was given earlier as 24,500: the newer 
calculated from the average of numerous Measurements of transmittar 
differed by 0.1 or 0.2 percent 

‘ To permit comparisons, values of the molar absorbancy index (wl 
into account the leneth of the licht path as well as the molar concent 
ire plotted along the vertical axis in most of the ficures. Strictly snea 
course, the term ‘‘molar absorhancy index” does not apnoly to mixture 
pounds, such as those in which bromophthalein magenta E has been | 
transformed into a salt 

5 Diphenylguanidine has a strong absorption band in the ultraviolet 
but an approximately 0.01-M solution in benzene does not show me 
absorption at wavelengths longer than 400 my with the cells used. Slig 
s sometimes evident in more concentrated solutions, which may ind 
trace of some oxidation product of aniline The spectral absorptior 
such substance should be nullified by the addition of equal amounts of diy 
guanidine to the “solution’’ cell and the reference cell. Absorption 
diphenylguanidine is probably shifted to shorter wavelengths when it ! 
formed into a salt The effect of such a shift, which would be very small ! 
of the low concentration of salt formed in any of the solutions, would 
detectable in the spectral region investigated. 
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WAVE NUMBER, CM X 10 . 
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700 


Changes wn molar absorbancy accompanying the 
eaction of bromophthalein magenta E, 2.5 10 V 


with diphenylquanidine 


RE | 
mary 
be nzene, 


and 1.0 molar 
equivalents 


5. with 0.2. 0.4, 0.6, 0.8 


1 to 2 
with 45 molar 


respectively; S 


BPM-E without added bas« 
ilent of diphenylguanidine, 


nylguanidir 


my is 
cX- 


ive 1), in which the absorbaney at 540 
plotted against the concentration of the base, 
pressed as molar equivalents of the indicator. The 
successive changes in color are most evident for 
small increments of diphenylguanidine; a concentra- 
tion at least as small as 5 107° M can be detected 
readily, either visually or spectrophotometrically. 

Experimental data and calculated values of the 
association constant for the primary reaction of 
bromophthalein magenta E with diphenylguanidine 
in benzene are given in table 1. The initial concen- 
tration of bromophthalein magenta E, C,, was 2.5 
10° M. The molar concentration of diphenylguani- 
dine at the moment of mixing is given in column | 
and the measured absorbancy after mixing, in column 
2. The equilibrium concentration of the primary 
salt, [S,], is C,(A,/0.5100). [A]=C,—[S,], and [B] 

( [Si]. Values for Aoassn. 1) are given in 
column 6; the average value is 2.210°. Values 
calculated similarly from data at 405 my, with the 
average value 2.3 10°, are given in column 7; we 
consider these values to be less accurate, because 
the changes in absorbancy are smaller at this wave- 
length, and the calculations involve fwo limiting 
values that may be slightly in error (one for the 
indicator itself and one for the primary salt). As 
indicated in the table, the average value of Ags... Was 
computed from the data for increments of diphenyl- 
guanidine ranging from 1.0 107° M to 3.75 107° M, 
inclusive. With very small increments of the base 
or increments close to the limiting concentration, the 
effect of any experimental error is magnified; also, 
there is evidence that diphenylguanidine is adsorbed 
to a small extent on glass, and errors from this 
source would be greater for low concentrations. 
According to studies by Hunter and Marriott [22], 
diphenylguanidine is somewhat associated in naph- 


fsec eq 


thalene, through hydrogen bonding; any such associ- 
ation in benzene has no evident effect upon its reac- 
tion with bromophthalein magenta E at the concen- 
trations under consideration. 


° 
b 


° 
wu 


ee ee ee ae 


° 
Ld) 
Sowa 


ABSORBANCY AT 540 Mp 


°o 
re ere ee 





4 i 4 


10 20 30 40 
MOLAR EQUIVALENTS OF BASE 





Fiat RE - Change wn absorbane yY act ompan jing the primary 
reaction in benzene of 2.5 %10~°-M bromophthalein magenta 
E with I, and II, triphenylquanidine; 


J-m m cells 


diphen ylguanidi ne 





-2.0 





4.0 6.0 


—LOG [8) 


Ficure 3. Graphical evaluation of the association constant 
for the primary reaction of bromophthalein magenta E, 
2.5X10-§-M in benzene, with I, diphenylguanidine and 11, 
triphenylguanidine. 

5-mm cells. Numbers adjacent to experimental points indicate the number of 
moles of base per mole of the indicator. Double circles, values tor 540 my; solid 
circles, values for 405 mp. Values indicated by a cross inscribed in a circle are 
for the primary salt of diphenylguanidine and bromophthalein magenta E 
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TABLE 1 


Experimental data and association constant for the primary reaction of bromophthalein magenta E with diphenylgu 


in benzene 


ution 
average « 


ter 


) 


f 10 readings 


1 asterisk were no 


A graphical means of evaluating the association 
constant is shown in figure 3 (line 1), in which log 
[S,]/LA] is plotted against log [1, 20}. As 
required by the reaction the experi- 


1) 
(Bi 
postulated, 
mental values fall on or very close to a line of slope 

1. As the line has been drawn, the intercept 


corresponds to the association constant 2.2 10°. 

From a solution that contained approximately 
3.2 10~°-M bromophthalein magenta E and 0.015-M 
diphenylguanidine in benzene, rosettes of short, 
lustrous green needles deposited after several weeks. 
The supernatant liquid was decanted off, and the 
crystals were washed four times with benzene by 
decantation, then dried in a vacuum oven at room 
temperature. This solid was used in the prepara- 
tion of three different stock solutions in benzene, 
on the assumption that it was the primary salt of 
bromophthalein magenta E and diphenylguanidine 
(molecular weight, 873.26). Absorbancy curves 
were obtained for these three solutions and for 10 
others made from them by dilution. Within ex- 
perimental error, the curve for a 2.5107°-M 
solution of the salt was the same as that for a 
solution made by mixing one molar equivalent of 
diphenylguanidine with 2.510-°-M BPM-E. Be- 
cause of the variation in the concentration of the 
solutions and the different cell lengths employed, 
the absorbancies were converted to the molar basis 
in computing the constant for the dissociation of 
S, into A and B. The data obtained for one of the 


From measurements at 


540 mw 


B)s 


hohe her 
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’ 


tancy 


he omputation of the 


three series of solutions are presented in table 2 
Values for Kain. are given in column 5 and values 
for its reciprocal, Kosen., in column 6. The averag 
value thus found for K,,,,. is 2.3410°. 
average based on all 13 solutions was 2.35 10 
This is nearly the same as the value obtained b 
mixing the indicator with different increments 0 
diphenylguanidine; we regard it as less accurat 
because the measurements were made at a perio 
when the temperature of the room was subject 
more than the usual variation. Molar absorbance 
curves for most of this series of solutions are show: 
in figure 4. They show fairly sharp 
points near 355 my and 455 my (compare fig. | 
A third less well-marked intersection occurs nea 
295 my. In figure 3, nine values for the salt ar 
indicated by crosses inscribed in circles; four ad 
ditional values are not shown because they coincide 
with ones obtained for mixtures of the indicato 
with the base. The values show a little scattering 
but are close to line I. 

In the study of the reaction of triphenylguanidin 
with bromophthalein magenta E, it was immediatel) 
apparent that triphenylguanidine is much less reac- 
tive than diphenylguanidine. (See fig. 2, curve I! 
It was also evident that a secondary reaction overlaps 
the primary reaction. Consequently, in evaluating 
the primary association constant it was assumed that 


isoshe st 


* It will be observed that the vertical scale of fig. 4 is twice as long as t 
fig. 1. 
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Fieurge 4. Effect of dilution on the molar absorbancy of the 
primary salt of bromophthalein magenta E and diphenyl- 
guanidine in benzene. 
Concentration of the salt in moles per literx105: (1) 15, (2) 12, (3) 9, (4) 7.5, 


6, (6) 5, 


(7) 4, (8) 3, (9) 2.4, (10) 1.6. 
1ivalents of diphe ny lguanidine 


(S:) 2.5<X10--M BPM-E with 45 molar 


TaBLeE 2. Dissociation constant of the primary salt of bromo- 
phthalein magenta E and diphenylguanidine in benzene, and 
the corresponding association contant 


: : . ae 
Co at 540 my [s'}* [A]}Pe Kedicen Kaeo 

15.0X10-8 34. 00108 | 12. 50X10-5| 2.500K10-5 | 5.0010-* | 2.00105 | 

12.0 33.77 9.932 2.068 4.30 2.32 } 

‘9.0 33 10 7.301 1,609 3.95 2.53 | 
{7.5 32.29 5.936 1.564 | 4.12 2.43 
‘6.0 31.49 4.631 1.369 4.05 2.47 
{3.0 27.96 2.056 0.9441 4.34 2.31 

124 26.97 1.586 8135 | 4.17 2.40 } 

Limiting | | 

value_. 40.80 | 

; | 

Average 4.28X10- | 2.3410 


* Initial concentration of the primary salt in moles per liter of solution. 
In moles per liter of solution. 


B}=[A]. 
Reciprocal of Kaisse. 


* Stock solution. | 
Prepared by dilution of stock solution 
| 


ihe limiting curve for the reaction is the same as the 
ne found for the reaction with diphenylguanidine 
the one designated S, in all of the figures). This is 

very reasonable assumption, as trie ‘thylamine, 
iphenylguanidine, and di-o-tolylguanidine have all 
een found to have the same limiting curve in reacting | 
ith bromophthalein magenta E. (See [1], fig. 8.) | 
he experimental data and calculated values of the | 
sociation constant for the reaction of bromo- | 
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Ficure 5. Effect of triphenylguanidine on the molar absorb- 
ancy of 2.5 X10-*-M bromophthalein magenta E in benzene. 
5-mm cells. A, BPM-E without added base; 


and 1,600 molar equivalents of triphenylguanidine, respectively; 
BPM-E with 45 molar equivalents of diphenylguanidine. 


1 to 6, with 25, 50, 100, 200, 400, 
Si, 2.5X10--M 


phthalein magenta E with triphenylguanidine in 
benzene are given in table 3. The average value for 
Koesn. obtained from data at 540 my is 525; from con- 
siderations similar to those discussed for diphenyl- 
guanidine, this is probably more nearly correct than 
the value 547, which was obtained from data at 405 
my. 

Graphical evaluation of the association constant 
for the primary reaction is shown in figure 3, line II. 
The postulate that the primary reaction corresponds 
to the simple pattern, A+ B=—S,, and the assumption 
of the same limiting curve as for the reaction of 
BPM-E with diphenylguanidine appear _ well- 
grounded, as the values fall on or very close to a line 
of slope —1; the intercept of this line with the z-axis 
for log [S,]/[A]=0 corresponds to the association con- 
stant 525. 


TABLE 3. Experimental data and association constant for the 
primary reaction of bromophthalein magenta E with triphenyl- 
guanidine in benzene 


From measurements at— 





Cy 540 my 
euceadieaes — seen 405 my 
Kame. | 
A,> [Si] * [A] * [B}» K seen 

510-5 0.0592 | 0. 2002K10~-5 | © 2.210K10-5, 24.71K10-5 531 
37 5 . 0819 e 2.009 37.10 516 
5O . 1055 ¢ 1. 983 49. 48 527 
ft) 1080 41.991 49. 47 537 
62. 5_. . 1246 e 1, 889 61. 89 522 
75 . 1445 41.812 74. 29 526 
87.5 . 1612 41.730 86. 71 527 
100 . 1746 e 1.644 99. 14 525 
112.5 1894 41. 592 111.6 522 
125 . 2026 41.527 124.0 524 
175 . 2447 | 1 41. 321 173.8 522 
Average. .... 525 547 


* In moles per liter of solution. 

> Cale ulated from the average of 10 readings of the transmittancy. 
¢ C,=2.5X10— mole/liter; limiting value, 0.5100 

4 C,=2.52X 10 mole/liter; limiting v: alue, 0.5143. 








Representative absorbancy curves for the primary With the inclusion of the two associatioy 
reaction of bromophthalein magenta E with tri- | stants reported in this paper, a total of 11 asso: 
phenylguanidine (see fig. 5) closely resemble those | constants have been obtained at this Burea 
for the primary reaction of BPM-—E with diphenyl- | spectrophotometric data. These values are su 
guanidine (see fig. 1). The curves in figure 5 show the | rized in table 4. Some of these are provisional \ 
same isosbestic points near 355 and 455 my. It will | having been determined under less favorable « 
be noticed, however, that curves 2 through 6 show | mental conditions However, the redeter 
greater absorbancy than that of curve S,; at wave- | Value for the association of bromophthalein m:; 
lengths longer than about 560 my; the corresponding | E with diphenylguanidine agrees with the 
solutions had a purple tinge. These effects are indic. | PFo* isional constant [1], and it is probable th: 
other provisional values are at most only sligh 
error. 

The value given in table 4 for the associat 
N,N-dimethylaniline and hydrogen chloride in ¢| 
roform was calculated at this Bureau [23] from 
solutions of the salt of bromophthalein magenta E | trophotometric data obtained in another labor 
and diphenylguanidine [24] 


ative of the overlapping of the primary with a second- 
ary reaction. Similar effects shown in figure 4 (no- 
tably in curves | through 3) are very likely caused by 
a partial secondary reaction in the most concentrated 


1) 


sociation conat ifs primary acid base reactions in inert solvents 


lrinitro-m-cre ol 
Bromophthalein u 
Pieric acid 
Bromophthalein m 
do 
Bromophthalein 
Bromophthalein m 
do ce 
Bromophthalein m 
Bromophthalein magent 
Bromophthalein m 
Hydrogen chloride 


It is of particular interest that the association | this point, a number of examples of dipole mome 
constant for the reaction of tribenzylamine with | of such compounds, taken from the work of ot! 
picric acid in benzene is the same as the value ob- | investigators, have been assembled in table 5. T 
tained from measurements by Maryott of the | dipole moment in Debve units (uX 10" esu) is ¢ 
dielectric constant of benzene solutions of tribenzyl- | in the third column. The first five examples 
ammonium picrate [25]. The suggestion has been | Lewis salts, and the next group are Bronsted sa 


made (26, 27] that tribenzylammonium picrate may | The reader will note that the dipole moments of | 

exist in isomeric forms, R;sN.HPi and R;NH*Pi, | tributylammonium halides are of about the sai 

which are slowly interconvertible, but the arguments | magnitude as the moments of the trialkylami 

used in support of this hypothesis are not entirely | sulfur trioxide addition compounds and of pyridit 

convincing. No anomalies that might be explained | sulfur dioxide. The dipole moments of trialkyla 

on this basis were observed in spectrophotometric | monium picrates are about one and a half times 

studies of the reaction of tribenzylamine with picric | great. In contrast, the dipole moments of 

acid in benzene [20], and Maryott, in measuring the | constituent acids and bases are considerably small 

dipole moment of tribenzylammonium picrate, could | showing that in the formation of the salts ther 

not detect such an equilibrium [25]. Witschonke | appreciable separation of charge. 

and Kraus [28] were also unable to furnish experi- The similarity in the reactions that lead to 

mental support for such an isomerism. formation of the two types of salts is illustrated 
One of the authors has already pointed out [23 the following examples: 

that in an inert solvent there is no essential difference 

in the product of the reaction of a base with a Et,N:+S0,.—Et,N+*t:SO, 

“Bronsted” acid and the product of the reaction of a 

base with a ‘‘Lewis’”’ acid. Both are in effect Et;N:+HCI—Et,N:H?t...Cl- 

highly polar addition compounds.’ To illustrate 

primes aces In the first reaction, the separation of charg 

there may be some disagreement with this terminology. Kraus and Fuoss nave | DFOught about by displacement of the “lone 

spplied the term “‘ion-pair’’ to addition compounds of the type listed as “Brénsted | of electrons on the nitrogen toward the sulfur 


salts”’ in table 5 [29 For other points of view in regard to the term “‘salt’’, : } 
see [30], p. 159 and [31], p. 88 In the second reaction, the proton of hyd 
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Dipole moments of salts and thei 
and bases 


Compound Solvent | ~X10!'§ esu 
Lewis salts 











methylamine-sulfur dioxide Benzene °4.9 
ethylamine-sulfur trioxide do 6.84 
Do Dioxane one 
butylamine-sulfur trioxide Benzene r 
vridine-sulfur dioxide do °6.2 
Bre@nsted salts 
butylammonium chloride Benzene 
Do roluent e7.3 
ributylammonium bromide Benzene 67.61 
butylammonjum iodide do 8. 09 
ethylammonium picrate do $11.7 
ripropylammonium picrate Dioxane 12.7 
butylammonium picrat Benzene 111.9 
Do Carbon disulfide © 10.7 
Do i Dioxane 12.2 
Do Toluene 11.7 
soamylammonium picrats Benzene 111.9 
Dy Carbon disulfide © 10.7 
Do Carbon tetrachloride il. 1 
Do Dioxane © 13.6 
liethylbenzylammonium picrate Benzene 11.8 
rribenzylammonium picrat do £12.0 
A mine 
Crimethylamine None (gaseous state & 0.65 
lriethylamine Benzene Os 
+-butylamine do 0.77 
Pyridine ck «2.3 
Dimethylaniline do «1¢ 
Acids 
Sulfur dioxide None (gaseous state | 1. 61 
Sulfur trioxide do «0 
ivdrogen chloride lo £1.06 
Hydrogen bromide do «0.79 
H vdrogen iodide do «0.38 
Piecrie acid Benzene £1.75 
Reference [24]; t=25° C 
J. A. Geddes and C. A. Kraus, Trans. Faraday Soc. 32, 585 (1936): ¢=25° C 
V. Deitz and R. M. Fuoss, J. Am, Chem. Soc. 60, 2394 (1938); ¢=35° C 


4. A. Maryott, J 
P. Mortier, W 
25 we ¢ 
« L. G. Wesson, Tables of electric dipole moments 
Research, Massachusetts Inst rechnology, ¢ 


Research NBS 41, 1 
s-Natuurkund. Tijdschr 


1948) RP 1806; 1 =30° ¢ 
10, 114 (1941); f=25° C 


for Insulatior 
Mass., April 


Laboratory 


itute of imbridge, 


chloride is attracted toward the lone pair of electrons 
on the nitrogen; it should not be regarded as trans- 
ferred completely to the nitrogen, however, but as 
forming a bridge between nitrogen and chlorine. 
This means that in each of the examples given, the 
reaction is to be thought of as yielding only one 
product. Representation of  triethylammonium 
chloride, as well as other Bronsted salts, as polar, 
hydrogen-bonded addition compounds is consistent 
with the conclusion of Walden [30] and of Kraus [32] 
from conductance measurements, that in inert 
solvents such salts are extremely weak electrolytes. 
According to the Brénsted-Lowry definition 
scheme, the reaction between triethylamine and 
hydrogen chloride would be formulated as follows: 


HCl—H* +4 


(acid; ) 


Cl 


(base, ) 


(1) 


constituent acids 
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Et;N:+H*t—Et,NH* 


(acid) 


Et;,N:+ HCI—Et,;NH*+ Cl 


(bases) (acid, ) 


( base, ) 


(3) 


(acide) ( base, ) 
















That is, the reaction would be depicted as the simul- 
taneous release of a proton by hydrogen chloride 
(eq 1) and addition of a proton by triethylamine 
(eq 2), the net reaction (eq 3) being the complete 
transfer of a proton from hydrogen chloride to 
triethylamine. In such a reaction, the product 
would consist of two discrete particles. In con- 
sidering the question of acidity and basicity ip 
benzene in an exploratory way [12], Bronsted clearly 
conceived of acid-base reactions as conforming to 
eq a. However, he recognized that the ions of a 
salt in benzene and other inert solvents are almost 
completely associated electrically and spoke of a 
transition from the electrovalent to the covalent 
state. His concluding statement was “The provisory 
nature of these measurements must be expressly 
pointed out. An exhaustive study of the molecular 
state in benzene and similar solvents must precede 
the exact quantitative investigation of acidic and 
basic strengths in these solvents.” 

There is much corroborative evidence in the 
literature that hydrogen acids form addition com- 
pounds with bases in the absence of water and other 
bases, and recognition of this fact is implicit, if not 
actually expressed, in many articles. For instance, 
in eryoscopic studies of the dissociation of pyridine 
o-chlorophenoxide or quinoline o-chlorophenoxide 
in benzene and in p-dichlorobenzene, Glass and 


Madgin formulated the equilibrium in terms of 
the association BN+HOR@BNHOR [33]. The 


acid-base dissociation of phenylammonium, phenyl- 
dimethylammonium, and pyridonium picrates in 
nitrobenzene was studied by Witschonke and Kraus 
by the conductance method [28], and they recognized 
that the dissociation constant ‘“‘provides a quantita- 
tive measure of the relative strength of acids or 
bases in non-aqueous neutral solvents.’’ However, 
it appears to us more logical to express such acid- 
base equilibria in terms of an association constant, 
and we believe that if the acid-base equilibria studied 
by various other investigators had been formulated 
in this way instead of as a dissociation, the analogy 
in the behavior of Lewis acids and Brgnsted acids 
with bases would not have gone unrecognized for 
so long. To think of acid-base equilibria in terms 
of dissociation was natural, of course, because of the 
customary way of expressing the relative strengths 
of acids and bases in terms of ionic dissociation 
constants in water, especially when the starting 
material was so frequently a salt (which is usually 
the most convenient substance to study). 
Although, as just indicated, a number of investi- 
gators have arrived independently at the conclusion 
that acid-base reactions in inert media are addition 
reactions (and thus conform to the Lewis concept), 
the problem remains of harmonizing this view with 










the behavior of acids and bases in solvents such as 
water. In their well-known book [31], Luder and 
Zuffanti discussed the reaction between hydrogen 
chloride and water and, as a means of unifying the 
Lewis and Bronsted-Lowry concepts of acidity, 
suggested that formation of the coordinate bond and 
ionization may occur simultaneously, as indicated 
by the following equation 


H:0:<7 +H:Cl:—#H:0:H*!+:Ck:-! 
H H 

base acid 
However, it seems to us that a clear understanding 
of the mechanism of such reactions is not possible 
without consideration of the effect of any excess of 
base. In aqueous solutions, the solvent itself plays 
the role of a base, as has long been recognized, but 
in order to render the formulation of a reaction as 
simple as possible, the role of the solvent is usually 
ignored. This concession to convenience has had 
additional justification in the lack of concrete 
knowledge of the extent to which the solvent becomes 
bound to dissolved particles, but it has led to incor- 
rect interpretations. A partial step toward acknowl- 
edging the role of water has been made in the now 
quite general practice of indicating the formula of 
the hydrogen ion in aqueous solutions as H,O* or 
H,OH*. Some time ago, however, Huggins [34, 
35] concluded that hydrogen ions are best con- 
sidered as H,OHOH,* ions. If this view is accepted, 
the stepwise reaction of hydrogen chloride with 

water can be formulated as follows: 


H,O+ HCIA@H,OH*...Cl-, (1) 
base acid 
H,OH*...Cl-+H,O@H,OH* ..OH,+CI-. (2) 


We believe there is strong experimental support for 
this formulation of the reaction. The product of 
the first step is an addition compound with a hydro- 
gen bridge between oxygen and chlorine. The second 
step is the displacement of the chloride ion by water, 
not necessarily because water is the stronger base, 


but because it is present in such large excess. The 
proton now forms a bridge between oxygen and 


oxygen. (The chloride ion may also be solvated, 
but such solvation is probably of secondary im- 
portance.) The first step would naturally be ob- 
scured when hydrogen chloride is added to a large 
volume of water, and the apparent reaction would be 


2H,0+HCl2@H,OHOH,+Cr-. (3) 


In a similar way, the reaction of an ammonia 
derivative, say trimethylamine, with water can be 
formulated as a two-step reaction: 


(CH,),.N +HOH2(CH,),NH*.. OH (4) 
base acid 
(CH,),NH*.. OH-+-H,0<(CH,),NHOH,+ 0H 
(5) 
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In this case also, the first step is addition and th, 


second step displacement. The over-all rea: tion 
then becomes 
(CH;),N +2H,O=(CH,),NHOH,+OH-. f 


The first step is essentially the same as in Moor 
and Winmill’s formulation of the reaction of a 
amine with water (36, 37], but the second step diffe 
in being a displacement rather than a dissociation 
In our view, the formulation currently in vogue fo; 
such a reaction, namely, 


sf 


(CH,),;N +H,O* =(CH,),NH*+H,0, 


obscures the real nature of the reaction.* 

We have presented quantitative evidence that an 
addition product is formed when a base and 4 
Bronsted acid are brought together in an_ inert 
solvent like benzene. This evidence provides support 
for the first step in our proposed two-step mechan- 
ism for the reaction of an amine with water. What 
experimental evidence is there to justify the second 
step? A few examples will illustrate how a variet, 
of approaches have led to the conclusion that when- 
ever the proton is attached to an ammonium nitro- 
gen, it has a marked tendency to form a bridge to a 
second electronegative atom. 

Walden and Birr measured the mobility of thre: 
isomeric ions—tetramethylammonium, diethylam- 
monium, and mono-n-butylammonium—in water 
and in methyl cyanide. In methyl cyanide, all 
three ions had nearly the same mobility. In wate: 
the mobility of tetramethylammonium ion was th 
greatest, that of butylammonium ion was the 
smallest, and that of diethylammonium ion was 
intermediate. The authors concluded that the thre: 
ions are not solvated to any extent in methyl cya- 
nide, and tetramethylammonium ion is not solvate: 
in water, but incompletely alkylated ammoniun 
ions are solvated in water. The number of wate: 
molecules surrounding butylammonium ion is greate! 
than the number surrounding diethylammoniun 
ion [38]. We may picture butylammonium ion as 
having the formula BuN(HOH,);* and diethylam- 
monium ion as having the formula Et,N(HOH,), 
Huggins expressed the view that “‘the larger th 
number of [hydrogen] bridges connecting a given 
nitrogen to neighboring oxygens, the weaker eac! 
bridge is, hence it may be that an ammonium io! 
in water solution is a terminal for only one, two or 
three bridges, on the average, at ordinary temper- 
tures, instead of four .’ [35]. However, Trotma 
Dickenson, on the basis of later evidence, conclude 
that solvation by water is greater for secondat 
ammonium ions than for tertiary ammonium ions 
and still greater for primary ammonium ions [39 

Measurements of the ionic dissociation of various 
completely substituted and partially substitute: 


+ + 
* The reaction (CHy)s3N+ H:OHOHs2(C Hy);NHOH?+H20 may al 
but probably plays a minor role because of the minute concentration of 
present in water rhis reaction, of course, is a displacement 
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ammonium salts in different organic solvents, per- 
formed by Kraus and Fuoss and their collaborators, 
led to the following conclusions with respect to the 


latter group of salts [40]: “In solvents that have no | 


affinity for the proton, there is a specific interaction, 
due to the proton, between the positive and the 
negative ion. In a solvent such as pyridine, which 
has a marked affinity for the proton, interaction 
takes place between the proton and the solvent 
molecules, rather than between the proton and the 
negative ion. Further investigation is required to 
fully elucidate the interaction that occurs in systems 
of this kind.” 

Earlier in this paper, mention was made of the 
secondary reaction that overlaps the primary 
reaction of bromophthalein magenta E with tri- 
phenylguanidine and that occurs when bromophthal- 
ein ;magenta is mixed with a large excess of diphenyl- 
guanidine. As will be shown in the next section, 
spectrophotometric studies of this secondary re- 
action have furnished new and more concrete evi- 
dence of proton-bridge formation between ammo- 
nium nitrogen and other electronegative atoms. 


2. The Secondary Reaction 


The curves in figure 6, obtained for mixtures of 
bromophthalein magenta E and diphenylquanidine 
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Figure 6. Change in molar abserbancy of benzene solutions 
th variation in the concentration of bromophthalein magenta 
I and diphen ylguanidine 
‘:, 2.5X10-"-M BPM-E with 45 molar equivalents of diphenylguanidine; 1, 
r-M BPM-E with 0.04-M diphenylguanidine; 2, 2.5x10+-M BPM-E 
vith 0.0675-M diphenylguanidine; 3, 1.0X10--M BPM-E with 0.045-M di- 
guanidine; 4, 3.2x10°°-M BPM-E with 0.015-M diphenylguanidine 





in benzene, show that at high concentrations of the 
indicator and base, the maximum absorbancy shifted 
from 540 to 600 mg. The accompanying change 
in color was from magenta to a deep purple. The 
curves in figure 7, obtained for mixtures of BPM-—E 
and triphenylguanidine, exhibit a similar shift at high 
concentrations. If, instead of adding an excess of 
diphenylguanidine to the benzene solution containing 
S,, the benzene was largely replaced by ethanol 
(figs. 8 and 9) or by pyridine (fig. 10), the maximum 
absorbancy shifted in a similar way, and the color 
changed to blue. In figure 9, the curve obtained for 
a mixture of 2.5X10~°-M bromophthalein magenta 
E with 50 molar equivalents of diphenylguanidine in 
96 percent ethanol-4 percent benzene by volume 
(curve 2) is practically identical with the molar 
absorbancy curve for the tetraethylammonium salt 
of BPM-E in ethanol (curve 3). It was not possible 
to prepare sufficiently concentrated solutions of 
bromophthalein magenta E with either di- or 
triphenylguanidine to shift the color completely to 
blue, but it is evident from comparison of these five 
figures that for all four bases (diphenylguanidine, 
triphenylguanidine, ethanol, and pyridine), the shift 
in the curves was of the same general character, and 
that the limiting curve would be very similar to that 
of the tetraethylammonium salt. 
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Figure 7. Change in molar absorbancy of benzene solutions 
with variation in the concentration of bromophthalein magenta 
E and triphenylguanidine 


A, 2.5X10°-M BPM-E; S;, 2.5X10-'-M BPM-E with 45 molar equivalents of 
diphenylguanidine; 1, 2.5x10™-M BPM-E with 0.04-M triphenyleuanidine; 2, 
2.5<X10-%-M BPM-E with 0.2-M triphenylguanidine; 3, 5.0x10--M BPM-E 


with 0.2-M triphenylguanidine; 4, 2.0X10--M BPM-E with 0.05 M triphenyl- 


guanidine. 


505 





-| _ 

WAVE NUMBER, CM X 19° 
250 166.7 
~% 7 


200 


.4 


+6 
ah 5 


INDEX xX iO” 
o . 
eyevpouticevesssell 


uw 
oO 


b 
re 2) _ 


MOLAR ABSORBANCY 
w 
oO 








WAVELENGTH, MU 
FiIGuRE 8 Effect o henzene by ethanol on the mola 
absorbancy of a tion of the 
bromophthalein magenta E and diphenylguanidine. 
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The equilibrium constant for the secondary reac- 
tion of bromophthalein magenta E 10~°-.M) 
with from 0.01-M to 0.04-\/ diphenylguanidine in 
benzene was calculated on the assumption of the 


(2.5> 


following reaction 


S, (primary salt)+B (base)—S, (secondary salt). 

If diphenylguanidine (which behaves like a tertiary 

amine) is represented by the formula R;N and 

bromophthalein magenta E is represented by HX, 

an alternative way of writing this equation would be 

R,;NH X° +R ,N@(R;NHNR;)X 

The equilibrium constant corresponding to the first 

equation 1s 

S,] 

K a7 
* {S,) [B 

At 


concerned, 


the high concentrations of diphenylguanidine 
the concentration of A (unreacted 


WAVE NUMBER, CM'x 10° 
250 200 166.7 





iO> 


MOLAR ABSORBANCY INDEX X 








400 
WAVELENGTH, MU 
Ficure 9. 
Vf bromop! 


Effect of ethanol on molar absorbancy. 
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thalein magenta E in benzene 
2.5<x10°-M BPM-E with 3) molar equiv 
»,, 2.5K10°-M BPM-E with 
”) Dere t ethanol-4 percent 


BPM in ¢ 
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bromophthalein magenta E) was negligible, and |b 
was essentially the same as C,. The equilibr 
constant for the secondary reaction can therefor 
formulated as the ratio of the concentrations of § 
and S,, divided by C,, or 
a. 
Soe 


» 


kK, 


Computation of A, would be a relatively simp 
matter if the limiting curve for the secondary salt 
were the same for all compounds in all solvents 
From data at two wavelengths it would be possibl 
to set up simultaneous equations and calc 
S,/S, for each solution [41]. However, it woul 
surprising if the anion, which is responsible fo 
absorption of light in the visible region, remain 
unaffected by changes in environment, such as t! 
effective diameter of the cation and the type an 
amount of polar molecules present.'° For this 


| 


See, for example, [40 luding the remarks of K, Fajar 
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culation of K,. Our acce pted value for the absorb- 
ancy _ S, at 600 my is 0.0857. 
» at 600 mu was found by t 


The limiting ays 
trial and error 


WAVELENGTH, MJ 
IGURE 10, Effect of replacing benzene by pyridine on the 
BP) molar absorbancy of a solution of bromophthalein magenta E 
iy and diphenylguanidine. 
~ S 2.5X10-$ M BPM-E and 45 molar equivalents of diphenylguanidine in 
. ene; 1, 2.5<X10-°-M BPM-E and 50 molar equivalents of diphenylguanidine, 
reent by volume of pyredint 2, same as 1 but 30 percent of pyridine; 3, same 
2 but 500 molar equivalents of diphenylguanidine; 4, 5.0x10--M BPM-E, 
2 olar equivalents of « tiphenylguanidine, 80) percent of pyridine; 5, tetraethy! 
id B nium salt of BPM-E in ethanol 
brium 
ore b son, it is desirable to avoid high concentrations of 
ol polar molecules in the spectrophotometric determi- | 
nation of A,. On the other hand, the concentra- 
tions must be great enough for precise measurement 
of [S,], [S.], and [2]. One may therefore expect that 
some ranges of concentration will prove more 
favorable than others for the determination of A, 
impli The complete molar absorbancy curves for two of 
vy salt the solutions used in estimating A, for bromo- 
vents phthalein magenta E plus diphe ny iguanidine in 
yssibl benzene are shown in figure 11; curve 1 is for the | 
culat mixture with 0.02-\/ diphenylguanidine, pet curve | 
ld 2 is for the mixture with 0.04-\/ diphenylguanidine. | 
or th The other eight curves, which showed ‘the same | 
raine isosbestic points near 370, 420, and 550 my, were 
= th omitted from the figure to avoid crowding of the 
e an curves. | 
is rea The data obtained at 600 mu were used in the cal- | 
| 
» | 


for S 


ioe 


1667 


WAVE NUMBER, CM | X 
Ea oh 


r—So 


S,] 
°o 


\ 
' 
' 


Ld w b 
°o °o °o 


MOLAR ABSORBANCY INDEX X IO” 
z 








600 700 

WAVELENGTH, MJ 

FIGURE Effect of an excess of diphenylquanidine on the 
molar of the primary salt of bromophthalein 
magenta E and diphenylguanidine in benzene. 

M BPM-E with 45 molar equivalents of diphenylguanidine; 1 
MBPM-E with 0.02-M and 0.04-M diphenylguanidine, re- 

calculated curve for the secondary salt 


1] 


absor banc u 


S,, 2.5K10 


and 2, 2.510 
spectively; So, 


be 0.510. Then S,/S,=—(A,—0.0857)/(0.510—A,). 
The experimental absorbancies at 600 my and the 
calculated values are given in table 6 The values 
derived for A, ranged from 15.2 to 15.8, with the 
average value, 15.5. (The value obtained for 0.03-M 
diphenylguanidine seemed to be out of line and was 
omitted in the computation of the average.) 


TABLE 6. 
secondar Y 
diphe nylquan idine 


Experimental data and association constant for the 
of bromophthale in mage nta E with 
hen zene 


reaction 


in 


0.010 


per liter of solution. C, 


* Initial concentration of diphenylguanidine in moles 
M/liter. Cell 


(initial concentration of bromophthalein magenta I 2.5X10 
length=5 mm 

t ee me rental v ilues at 600 my 

1,—0.0857; 0.0857 is the limiting value 

600 _ when 2.5 10 M in 5-mm cel! 

4 z, the assumed limiting value found by trial and 
(S:) at 600 my, =0.510; AA,/(z 1,)=S2/S 

¢Ky=S SS, 


" 
tThe value marked with an asterisk was excluded in the 
average 


found for the primary salt at 


error for the secondary salt 


computation of the 


Graphical a of the equilibrium constant 


is presented in figure 12, line I. The experimental 
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Figure 12. Graphical evaluation of the association constant 
for the secondary reaction of bromophthalein magenta E, 
25 %10--M in benzene, with I, diphenylguanidine and I], 
triphenylguanidine; 5-mm cells 


points (open circles) fall on or very close to a line of 
slope — 1; this line was drawn so that its intercept on 
the z-axis for log S,/S,=0 is 1.19 (log 15.5). The 
validity of the equation assumed for the secondary 
reaction is very strongly supported by the close 
agreement of the values for A, and by the almost 
perfect coincidence of the plotted points in figure 12 
with the theoretical curve, line I. 

From the data in table 6, the curve for S,, and 
the estimated value (0.510) of S, at 600 muy, the limit- 
ing curve for the secondary reaction of bromophthal- 
ein magenta E with diphenylguanidine was calcu- 
lated (fig. 11, curve S, The general shape of this 
curve is similar to curve 4 in figures 6 and 7, curve 
7 in figure 8, curves 2 and 3 in figure 9, and curve 4 
in figure 10. The method of computation was as 
follows: The increment (or decrement) in absorbancy 
due to the complete conversion of 2.5 107°-M S, to 
2.5 107°-M S, is 0.5100—0.0857 at 600 mu. At other 
wavelengths, A,(S,)-A,(S,) is unknown, but the 
ratio [A, (S,)-A, (S,)]/LA, (mixture)— A, (S,)] should 
be the same for all wavelengths, and A, (mixture)- 
A, (S,) is known for each of the solutions concerned 
over the wavelength region in question. The incre- 
ments or decrements were calculated for each of 
the solutions (except the one containing 0.03—M 
diphenylguanidine), and the average increment or 
decrement was added to curve S, to give curve S). 

Curve S; shows a maximum near 580 mg, not near 
600 my, although in figure 6 two of the curves for 
mixtures of BPM-E and diphenylguanidine (curves 
3 and 4) show a maximum near 600 mu. The two 
mixtures that gave these curves obviously contained 
a much larger concentration of S, than the solutions 
that gave curves 1 and 2 in figure 11, and it seems 
possible that at high concentrations the salt is ap- 
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Figure 13. Structures of bromophthalein magenta E and 


various of its salts. 


I, Bromophthalein magenta E; X=COOC:Hs; WU, Anion of bromophtha 
magenta E, showing the two structures that contribute chiefly to the :esonar 
III, Primary salt of bromophthalein magenta E with a tertiary amine; IV, Se 
dary salt of bromophthalein magenta E with a tertiary amine (mor 
structure); V, Salt of bromophthalein magenta E with a secondary amine 
structure); VI, Secondary salt of bromophthalein magenta E with a tertia 
amine (dimeric structure). 


preciably dimerized. A _ possible explanation for : 
shift of the absorption toward longer wavelengths 
when dimerization occurs can be made with the aid 
of figure 13. Formula [| in figure 13 is for bromo- 
phthalein magenta E itself; in bromophthal 
magenta E, X¥=COOC,H;. The benzene ring marked 
2 has a quinoid structure, which accounts for th 
vellow color of the solutions. Formula II represents 
the anion and shows the two structures that con- 
tribute chiefly to the resonance that causes absorp- 
tion at lower frequencies (longer wavelengths). I! 
the anion were completely detached from th 
phenolic hydrogen of structure I, the two structures 
within the brackets would be exactly equivalent 

Formula III is for a tertiary ammonium salt 0 
bromophthalein magenta E; the phenolic oxygen !s 
partly released from the restraining effect of th 
proton, but benzene rings 1 and 2 are not fully 
equivalent and the absorption band is intermediat: 
in position between those produced by structures | 
and II. This is the primary salt, S,, which gives 
a magenta solution in benzene. In all of the 
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mary salts so far investigated, the effect of the cation 
on the electron system of the anion appears to be 
the same. Formula IV is for the secondary salt, S,; 
the positive charge on the cation is further removed 
from the phenolic oxygen, and this explains why 
absorption occurs at lower frequencies. The loca- 
on of the positive charge will vary with the structure 
of the cation, and this will affect the resonance within 
the anion. It is also possible that a displacement 
of the electron cloud toward the oxygen of ring 2 


will occur if the oxygen is near the positive end of | 


some dipole, resulting in a further shift in the absorp- 
tion toward lower frequencies. 

Formula V was previously suggested [1] to account 
for the cornflower blue color of secondary ammonium 
salts of bromophthalein magenta E in benzene 
(absorption maximum near 575 my)."' The deeper 
color, as compared with the tertiary ammonium salt, 
S,, is understandable because rings 1 and 2 are 
more nearly equivalent in formula V than in formula 
Ill. Formula VI suggests how a similar shift in 
absorbancy to lower frequencies may occur at high 
concentrations of the secondary salt (formula IV); 
two adjacent molecules of the salt may become 
oriented in such a way that rings 1 and 2 are both 
adjacent to a cation and thus become more nearly 
equivalent than they are in formula IV. Such 
association is not purely hypothetical but has a 
strong experimental foundation in the work of 
Kraus and coworkers [42]. It is more likely to occur 
when the ions are symmetrical. 

The calculation of A, for the reaction of bromo- 
phthalein magenta E with triphenylguanidine in 
benzene was calculated from data at 540 my for 
nine solutions with the initial composition 2.5> 
10-°-./ bromophthalein magenta E and from 0.002-\4 
to 0.01-M triphenylguanidine. The calculation was 
subject to greater error than in the case of diphenyl- 
guanidine, because the changes in absorbancy were 
smaller and the primary and secondary reactions 
overlapped. The value 525 given for A, in the 
preceding section was assumed to be correct and, 
as in the calculations for the secondary reaction of 
BPM with diphenylguanidine, [B] was assumed to 
be the same as C,. Since A,=525, S,/A=525.C,. 
Then f=(S,/A)/[1+(S,/A)] gives the fraction of the 
indicator that would have been present as the 
primary salt in the equilibrium mixture if none 
of the primary salt had been transformed to the 
secondary salt; multiplying f by 0.5100 gives the 
absorbancy that would have been observed under 
these conditions. Similarly, if y is the absorbancy 
at 540 mu for 2.5 10~°-M secondary salt, f.y gives 
the absorbancy that would have been observed if 
all of the primary salt formed in the solution had 
been converted to the secondary salt. Then S,/S, 

Note the distinction in the terms “secondary salt" (formula IV) and “‘second- 
iry ammonium salt” (formula V). A very striking phenomenon occurs when a 
sic solvent such as ethanol is added in minute increments to the cornflower 
blue secondary ammonium salt in benzene. The color first changes to magenta 
and then, upon the addition of a large excess of ethanol, to the blue of tetra- 


ikylammonium salts of bromophthalein magenta E. These effects were de- 
scribed and explained in [1] 








(0.5100f—A,)/(A,—yf); dividing this ratio by C, 
gives A,. For y=0.28, the values of K, range from 
18.0 to 28.2, with the average value 22.5. The 
data and calculations are presented in table 7. In 
table 8 are listed the values for K, that were obtained 
for other values of y; the most reasonable value of 
y appears to be in the range 0.2 to 0.4. The value 
for K, thus seems to be of about the same magnitude 
as K, for the reaction of bromophthalein magenta 
E with diphenylguanidine, but somewhat larger 
rather than smaller. An attempt was made to 
compute A, from experimental values of the absorb- 
ancy at 600 my, but the values of AK, thus obtained 
showed a trend toward lower values with increasing 
C, whatever the value assumed for the limiting 
absorbancy at 600 my, although they were of the 
same order as those shown in table 8. The most 
likely explanation of this effect appears to be a 
relative increase in the concentration of a dimer 
of the secondary salt as well as an increase in polar- 
izing effects caused by the free base, as the con- 
centration of the base is made greater. Effects due 
to these causes would probably be more marked at 
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Ficure 14. Absorbancy curves of salts of bromophthalein 


magenta E. 


1, Tetraethylammonium salt in ethanol; 2, calculated curve for secondary 


salt with triphenylguanidine in benzene 
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K 520, 530, or 540 gave less good agreement 
of the values for A, at different concentrations 
of the base. but had little effect on the order of 


magnitude 


TABLE 7 


secon la 


data and association constant for the 
E with tri- 


kx pe 


reaction of ty 


tmenta 
omophth alein magenta 


phen jlquanidine in benzene 


NMwhwmwhMNhhN 


s very large 


constant for the 
with tri- 


TABLE &. Calculated tlues of association 


secondar j reaction 


bromophthalein magenta BE 
* 


phenylquanidine henzene 


on cor lary reaction [B] 
' C,) , )} is the assumed value 
e assumed v fc t ima hen 2.5 * 10-5 Min benzene in 5 
cell: f . ; s): the correct value for K, is assumed to 
e measured ney I t 540 my 
las probably the best value y (see the text) 


The points located by the values given for log 
S,/S,) and log (’, in table 7 are shown as the 
solid circles in figure 12. They are fairly well dis- 
tributed on either side of a line of slope 1 (line II 

The data for the secondary reaction of bromo- 
phthalein magenta E with triphenylguanidine were not 
good enough to warrant extensive calculations to 
determine the absorbancy curve for the secondary 
salt. A few trial calculations were made for the 
solution (0.004-M diphenylguanidine) that gave the 
point which falls closest to line II. For y=0.28, 
the curve obtained is the one shown as curve 2 in 
figure 14. The resemblance of curve 2 to the absorb- 
ancy curve for the tetraethylammonium salt of bro- 
mophthalein magenta E in ethanol at an equivalent 
concentration and cell length (curve 1) i 
as to lend further support to the assumptions in 


Is so great 


regard to the secondary reaction, although we d 
believe that it is an accurate representation o 
absorbancy curve for the secondary salt of BP 
with triphenylguanidine. An absorption band 
500 mu was indicated whatever value was used | 
within the range considered acceptable. Pos 
the complex cation contributes to the absor Or 
at this wavelength 


3. Further Discussion 


The important theoretical points made in th 
previous discussion may be summarized as follows 

1. When equivalent amounts of an organic base of 
the ammonia class and an acid of the Brénsted- 
Lowry type react in an “inert” solvent like benzen 
a highly polar addition compound is formed. This 
may be regarded as a salt composed of a pair of ions 
held together by Coulombic attraction as well as by 
interaction between the negative ion and the proto 
of the hydrogen ion. This addition reaction, termed 
by the authors “the primary reaction’’, is analogous 
to the combination of a base with a Lewis acid 
The Brénsted-Lowry acids should thus be regarded 
as a special class of Lewis acids. The principal 
arguments on which these statements rest are (a 
evidence that ionic dissociation is insignificant unde: 
the conditions referred to, (b) quantitative measure- 
ments of equilibrium constants for a reaction postu- 
lated as A+ Bs58), (c) spectrophotometric evidenc 
that the anion is not completely detached from th 
proton of the cation, and (d) the high 
moments of the addition products. 

2. When the reaction occurs in the presence of 
an excess of base, the addition reaction may be sw 
ceeded by a “secondary PP displacement reaction ll 
excess of base is essential for appreciable ionizatior 
(The anion displaced by any neutral bas: 
The primary 


dipol 


to occur. 
may, of course, be regarded as a base.) 
reaction is generally obscured in a polar solvent su 
as water. The existence of such reactions has bee: 
suspected by a number of investigators, but studies 
of acid-base reactions in inert solvents were require: 
to establish the reality of the primary reaction an 
the secondary reaction 

The studies reported in this paper, while helping to 
clarify the effect of environment on the course o! 
acid-base reactions, suggest other topics for st 
Investigations of primary and secondary acid-bas 
reactions should, of course, be extended to systems 
other than those dealt with in this paper. It woul 
be interesting, for example, to devermine the cor 
tions under which the primary and 
reactions are apt to overlap. One might logical! 
predict that overlapping will not be likely in al 
inert solvent when the value of A, is large. Stet 
factors will undoubtedly prove important; in qua! 
tive tests with some aliphatic amines, there was n 
evidence of the secondary reaction when the an 
contained methyl groups than when the sma 
substituent was an ethyl group. Such = stu 
obviously have a bearing on the interpretation 


second 


510 





the 
IWS 
~(° of 
ted- 
ene 
This 
ons, 
s by 
Lon 
med 
FOuUSs 
eid 
rded 
‘ipal 
a 
nder 
ire- 
Stu- 
pnce 
the 
pole 


suec- 
An 
tion 
base 
rary 
such 
ween 
dies 
aired 


and 
a to 


dy 

base 
ems 
ould 
ndi- 
lary 


liction of the behavior of catalysts as well 
bases of pharmacological importance. 

\ccording to the aqueous ionic dissociation con- 

uts, triethylamine is a_ stronger base than 

shenylguanidine (see [1], table 2), although the 

erse order of basicity is indicated by their relative 

dencies to combine with bromophthalein magenta 

or B in benzene. This would suggest that in 

ter the secondary reaction occurs to a greater 

tent for triethylamine. In benzene, diphenyl- 
guanidine shows the behavior of a very strong base; 
as shown in the preceding discussion, A, for its 
association with bromophthalein magenta E or B is 
quite large, and with picric acid (which is presumably 
a stronger acid than the bromophthalein magentas), 
it reacts ¢ uantitatively at concentrations less than 
5x10 M [43]. Its capacity for neutralizing acids 
in inert media and its ability to form stable salts 
are likely to be greater than would be predicted 
from the aqueous ionic dissociation constant. 

The low reactivity of triphenylguanidine with 
bromophthalein magenta E, as compared with its 
moderately great basicity in water ({1], table 
suggests possible steric hindrance due to the bromine 
atoms in bromophthalein magenta E. The sul- 
fonephthalein indicators analogous to bromophthal- 
ein magenta E—bromophenol blue, chlorophenol 
blue, and iodophenol blue—show the same pH range 
in water (see [44], table 1), but differentiating 
effects might well be found in the reaction of bromo- 
phthalein magenta E, chlorophthalein magenta E, 
and iodophthalein magenta E with bases in solvents 
like benzene. So far, the latter two indicators are 
not available for such studies. 

In recent studies of the effect of the basic ity of the 
accelerator in the vulcanization of rubber, Kratz, 
Young, and Katz found that “within a closely re- 
lated series of accelerators, the sulfur combined in a 
given time at a given te mperature is a linear function 
of the logarithm of the basic dissociation constant 
of the accelerator employed” [45]. Thus, mono- 
phenyl-, diphenyl-, and triphenylguanidine show 
decreasing accelerating power in the order named, 
Our studies indicate that an extension of this gen- 
eralization to bases not closely related would be 
liable to prove fallacious. 
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